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Howard E. 


Beginning with the June, 1959, issue, IRE TRANSACTIONS 
ON ELECTRONIC CompuTeERS will have a new editor. He is 
Howard E. Tompkins, assistant professor of electrical en- 
gineering in The Moore School of Electrical Engineering at 
the University of Pennsylvania, Philadelphia, Pa. 

Dr. Tompkins is especially well qualified for the post. He 
was one of the early workers in the digital computer field, and 
continues to be active in it. His education has been directed 
toward computer applications of electrical engineering. He re- 
ceived the B.A. degree in physics from Swarthmore College, 
Swarthmore, Pa., in 1942. 

From 1942 to 1947, he was employed by the Philco Cor- 
poration in their Research Division, working first on micro- 
wave resonators, receivers, power meters, and measuring tech- 
niques, then on magnetic recording techniques and phono- 
graph pickups. Several patents resulted. 

In 1947, he received the M.S. degree in physics from the 
University of Pennsylvania, and that fall, joined the staff of 
the University as Instructor in The Moore School of Electical 
Engineering. Until 1951, he taught undergraduate electrical 
engineering, took graduate courses, and worked, first on elec- 
tromedical instrument amplifiers, then on digital computer 
components. 

In 1951, he joined the Burroughs Corporation as head of 
a transistor circuits research group, working on transistor 
switching circuits for computers, including in particular tran- 
sistor-resistor logic. In 1955, he became Educational Services 
Manager at the Burroughs Research Center, returning to The 
Moore School in the fall of 1956. 

During 1956-57, he supervised graduate student research 
on transistor high-frequency measurements, and wrote his dis- 
sertation on Unit-distance Codes. He received the Ph.D. de- 


Tompkins 


gree in electrical engineering in 1957 from the University of 
Pennsylvania. 

Over the years he has been active in local IRE affairs. He 
was local Arrangements Chairman for the first IRE-AIEE- 
ACM Computer Conference, held in Philadelphia in Decem- 
ber, 1951. For many years he has been associated with the 
IRE-AIEE Transistor and Solid-State Circuits Conference, 
either as a member or chairman of the Program Committee, as 
Publicity Chairman, or as IRE Philadelphia Section represent- 
ative. He was chairman of the Philadelphia Chapter of IRE- 
PGCT for 1956-57. He was IRE-PGEC representative to the 
Program Committee of the 1955 IRE National Convention, at 
which he was organizer and moderator of the Symposium, 
“Design of Machines to Simulate the Behavior of the Human 
Brain.” (See IRE TrRANs. oN ELECTRONIC CoMPUTERS, vol. 
EC-5, pp. 240-255; December, 1956.) Dr. Tompkins is also a 
member of AIEE, ACM and ASEE. 

The two and one-half years during which I have been editor 
have been pleasant and stimulating ones. The interest in com- 
puters has caused a growth of our group that has brought the 
number of copies of these TRANSACTIONS published each quar- 
ter from 6300 to 8700, and there are positive indications that 
this is only the beginning. 

Under Dr. Tompkins’ editorship, we shall see these TRANS- 
ACTIONS continue its growth, and he will need the assistance 
of many people in carrying out his responsibilities. I know that 
those who have been helping me, especially in the arduous task 
of reviewing papers, will continue to give their support to Dr. 
Tompkins. I want to express my great appreciation to those 
who have given so generously of their time, and who must do 
so without public acknowledgment in order that our reviewing 
procedures may have the anonymity required.—J. P. NAsu. 
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A Magnetic Core Parallel Adder* 


MAO-CHAO CHEN ft 


Summary—A logical design using magnetic core elements which 
does not have the usual carry time limitations is described. The 
synthesis uses a truth-table technique. 


INTRODUCTION 


INARY computer units are usually analyzed by 
B Boolean algebra, then built from “And,” “Or,” 

and “Not” switches by various methods. In the 
present paper, a truth-table type synthesis is employed. 
The addition function to be described is performed by 
the voltage (flux) addition, subtraction and switching 
characteristics of magnetic memory cores. Previous 
work in this field has been done by Stifler [1], Mealy 
[2], more recently by Andrew [3], and others. An ex- 
tension of these early methods has led to an arrange- 
ment for the addition function resulting in a parallel 
adder that does not have the usual carry-time limita- 
tion. 

The voltage addition and subtraction operation em- 
ployed makes admissible both negative numbers and 
numbers greater than 1. For example in Fig. 1 the value 
at point S, (to be described later) is 2 when c=x=1, and 
y=0; the value c’ in (4) is —1, when x, y, and ¢ are 
zero. Numbers in Boolean algebra are restricted to 1 
and 0, and the logical operations “And,” “Or,” and 
“Not” are employed to mechanize the algebra in the 
binary system. Furthermore, Boolean functions do not 
generally describe voltage addition and subtraction. In 
consequence, Boolean algebra cannot completely de- 
scribe the behavior of a multilevel system, even though 
the resulting numbers (1 and 0) are the same. Therefore, 
in the following description, ordinary algebra and the 
usual meaning of addition and subtraction are employed. 


THE ADDER 


One stage of the adder is shown in Fig. 1. Let x and 
y denote the addend and augend respectively, c the in- 
put carry, and s and c’ the sum and output carry. All 
magnetic cores are identical, and rectangular magneti- 
zation curves are assumed. Each driving winding of the 
cores A, B, C, and D, has the same number of JN turns. 
Each output or secondary winding also has same num- 
ber of turns. A “1” on x or y is represented by a current 
pulse of magnitude J amps (flowing in the direction 
shown) providing a net amp-turn (JVJ) is large enough 
to drive each of the cores into saturation. An “0” for 
x or y provides no current at pulse time. The x- and 


* Manuscript received by the PGEC, November 27, 1956; re- 
vised manuscript received, March 4, 1958. 
+ Physics Dept., Stanford University, Stanford, Calif. 


Fig. 1—One stage of adder. The arrows on input windings x, y and 
auxiliary windings a, b show current directions. Arrows on output 
windings show directions of induced emf’s when the cores change 
flux. The concentric arc arrow ~ inside the two circles shows flux 
direction before operation. All x windings and y windings are 
in series; they are shown separately for clearness. 


V- pulse 


driving pulse 


t 


Fig. 2—Pulse obtained when a core is switched by NJ amp-turns. 


y-pulses are applied at the same instant of time. 

In these circuits, we are dealing mainly with pulses 
as shown in Fig. 2. For convenience we define a V-pulse 
as one which is identical with that shown in Fig. 2 not 
only in shape and magnitude but in time as well. We 
initially assume that the Thevenin’s equivalent circuit 
at the input carry c, looking towards lower bits, is a 
pulse generator and a diode d; in series as shown in Fig. 
1. The pulse is a V-pulse when ¢ is “1,” and zero when 
c is “0.” Acurrent pulse J is applied to the windings 6 of 
core D and a of core C at the same instant of time when 
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x- and y-pulses are applied. A driving mmf is considered 
negative when it tends to maintain flux in its original 
direction, and positive when in the opposite direction. 
Thus a core changes state when the net driving mmf is 
greater than NJ, and remains unchanged when it is 
zero or negative. Table I gives net mmf’s of the cores 
A, B, C, and D for various combinations of x and y. 


TABLE I 
NET MAGNETOMOTIVE FORCE As A FUNCTION OF x AND y 


core A core B core C core D 
x=0, y=0 0 0 NI —NI 
s=0;7=1 —NI NI NI 0 
s=1,y=0 NI —NI NI 0 
yy 1 0 0 NI NI 


Let us define “1” as the condition when an output 
winding on a core is a V-pulse and “0” when there is no 
pulse. From Table I, we see that core A changes flux 
only when x=1 and y=0. Hence, the output of core A 
is x(1—y). Similarly, the output of core B is y(1—x). 
Since the output windings of A and B are in series, the 
value at S; is 


Si= a1 — y) +91 — x2) = x+y — 2xy. (1) 


S; and input carry ¢ are applied to the grids of two 
vacuum tubes, the plate circuits of which are connected 
through two cores £ and F as shown in Fig. 1. Normally 
the dc plate currents flow through either of the cores in 
opposite directions, and their mmf’s cancel each other. 
When the S; pulse is applied alone, plate current of tube 
1 increases causing the net mmf’s in core E to become 
negative and core F positive. The amplifiers and driving 
winding turns are so adjusted that this change in mmf 
is large enough to switch core F. Similarly when input 
c-pulse is applied alone, core E changes flux. But when 
both S, and ¢ are applied, their mmf’s cancel each other 
-again and neither core is switched. Hence core E 
changes flux only when ¢ is 1 and 5S; is zero, thus it 
mechanizes the function, c(1—.S,). Similarly core F 
mechanizes S; (1—c). Since their output windings are 
in series, voltage at s is 


$= 9911 — c)'+ 661. — Si) 96 S81 — 268). (2) 
Substituting (1) for S,, we have, 
(2a) 


This is the equation for the sum when ¢, x, and y are 
added, which can be easily checked by a truth table and 
therefore is not shown here. 

From Table I, it can be seen that core C always pro- 
vides an output at “pulse time” and is therefore a “1” 
generator. The value at S; is the same as S| except ¢ is 
in series, thus S:=S,+c. The opposing direction of the 


S=¢ +444 — 2e% — 2cey — 2xy + 4exy. 
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output winding on core C, is selected so that the value 
avezis 


ec =S,-1=S,;+c¢—-1. (3) 
Substituting (1) for S, results in, 
c=x+y—2y+c-1. (4) 


When x, y, and ¢ are zero, c’ becomes —1. Due to diode 
dz, no negative voltage can appear at c’, so 


Omar y = Daya 1), (5) 


where the function f takes the value of its argument 
when it is positive, and zero when it is negative. From 
Table I, it is seen that core D mechanizes the logical 
SATE A (5).z 


C= ay, (6) 


Thus c’ consists of two branches in parallel, one from 
core C to A, B, carry c, and diode d to ground, the other 
from diode d;, core D to ground. These two branches 
could form a closed loop. Hence, diode d, is inserted to 
prevent circulating current within this loop. Hence c’ 
then provides the voltage of either of the branches and 
the greater is then the output. From (5) and (6) the 
complete expression for the output carry is 


bt cheat iu fea y—24ey—1) Soy 
xy, ey fC a eel 


This expression satisfies all of the cases when x, y, and 
c are added, and can be checked easily by a truth table. 

Several important points are to be noted in the above 
circuit: 

1) None of the logic cores or the equivalent input 
carry generator circuit is loaded. 

2) All of the cores (A, B, C, and D) are switched by 
the same number of amp-turns (VJ) and at the same 
instant of time. 

3) The Thevenin equivalent circuit from c’ toward 
lower bits is the same as that at ¢, 2.e., a pulse generator 
and a diode in series. 

The results 1) and 2) lead to the fact that all voltage 
pulses induced have the same shape, magnitude, and 
time instant as shown in Fig. 2, and are all V-pulses. 
This is essential so that the quantities in (7) can be 
manipulated arithmetically. Therefore the output carry 
c’ is also a V-pulse. Since the input carry ¢ and output 
carry c’ are identical in their pulse representation as 
well as their equivalent circuit, all stages can be con- 
nected in cascade, and a complete adder is formed. The 
two lowest stages are drawn in Fig. 3, the extension to 
many stages being obvious. Since there is no input 
carry in the first stage, the triode 1 and core F are pri- 
marily unnecessary. They are inserted there so that the 
sum pulse can appear at the same time as other sum 
pulses. 
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B+ 


C2° C 


emo 


Fig. 3—First two stages of adder. 


DISCUSSION 


As shown in Fig. 2, there is always a time delay be- 
tween the driving pulse and the output V-pulse. A re- 
view of the operation of this adder shows that this kind 
of delay does not effect its operation. The only effect is 
in the operation speed. In performing the addition, first, 
cores A, B, C, and D, are switched, then the output 
pulses of these cores are used to switch cores £ and F. 
Let T denote the time delay in switching a single core. 
The speed of operation is 2T. 

It isa well-known switching characteristicthat second- 
ary output pulses will be changed in shape, magnitude, 
and pulse width if secondary windings supply current. 
As stated earlier, no core in this adder is loaded. This is 
essential in order to provide uniform pulse shapes. 

Vacuum tube amplifiers are used to amplify and con- 
vert voltage pulses into current pulses to drive cores E 
and F. Transistor amplifiers draw current as well as power 
from the signal; thus the V-pulses would be distorted if 
they were used instead of vacuum tube amplifiers. For 
this reason, vaccuum tube amplifiers are preferable. 

Under no-load condition, the output voltage obtained 


in switching a core might be quite large depending upon - 


driving amp-turns. For instance, 2.5 volts can be ob- 
tained across a single output turn when 2 amp-turns is 
used to drive a General Ceramics F-483 ferrite core. 
When a high transconductance tube is used, for in- 
stance, a 6AG7, g=7700 umho, this voltage can be 
transformed into a current pulse of 2.5xX7700X10-* 
=0.019 amp. Thus 100 turns winding is sufficient to drive 
cores # or F. It is obvious that high primary impedances 
do not effect the operation of these cores. If in appli- 
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cation a large pulse across s is not necessary, cores of 
small sizes can be used for cores £ and F, and driving 
amp-turns can be greatly reduced. 

In Fig. 3, we see that a carry pulse propagates through 
windings of cores B, A, and C, of successive stages until . 
it arrives at a stage where both x and y are zero. Then 
it transforms into a sum pulse and disappears from the 
circuit. But, in general, it can propagate as far as to the 
highest bit. Due to impedances of output windings of 
these cores, the carry pulse is delayed and distorted 
during propagation, and the initial assumption that a 
carry pulse is a V-pulse can not be justified. This puts 
on the restriction that the number of bits to be added is 
limited when the effect of pulse delay and distortion is 
beyond tolerance. However, it is possible to reduce the 
delay effect by making use of what has already been 
shown; one turn secondary winding is enough to drive 
the amplifier. This leads one to arrange the A, B, C cores 
of successive stages coaxially, a center conductor serving 
the secondary windings of cores A, B, C. This looks like © 
a concentric transmission system with cores as dielectric 
and carry pulse propagating along the center conductor. 
It is obvious that time delay can be further reduced by 
using cores of less thickness. In order to keep a definite 
amount of flux to give sufficient voltage, core diameter 
must be increased accordingly. Hence the maximum 
allowable number of bits to be added depends upon par- 
ticular design and is best determined experimentally. 

Reset windings must be provided for all cores to set 
the cores to their initial conditions before the operation 
of the next addition. 
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Summary—The usual floating point arithmetic makes error 
alysis difficult. This paper describes an alternative system which 
offers a means of analyzing floating point calculations more effec- 
tively and which also possesses certain advantages from an equip- 
ment standpoint. 


} 
| 


INTRODUCTION 


| ODERN computers which make use of so-called 
M floating point arithmetic, where every number 
| is represented as a fractional part with expo- 
nent, customarily deal with such numbers in normalized 
form; 7.e., nonzero fractional parts are required to lie 
between B-! and 1 in magnitude, where B is the number 
system radix. The use of such a convention, however, 
leads to the retention of nonsignificant digits in com- 
puted results, which therefore gives them a highly de- 
ceptive appearance of accuracy with no indication of 
what their true precision might be. This paper describes 
some aspects of an alternative system, which may be 
termed significant digit arithmetic. First the rules for 
this arithmetic are described, and then methods of com- 
puter implementation of the rules are discussed. Al- 
though the basic philosophy underlying the significant 
digit system is indicated, a detailed justification is not 
attempted here. An analytic study of its error propaga- 
¢ion properties is in preparation. 

The following discussion applies quite generally to 

oating point arithmetic, irrespective of the number 
oe and specific representation conventions used. For 
concreteness, however, it will be assumed that numbers 
are represented in binary form; the fractional part f 
consists of an n-digit fractional magnitude with sign 
appended, and the integral exponent e is also repre- 
sented as a signed magnitude. Obvious modifications are 
necessary to make the discussion appropriate to the case 
of complement representation, biased exponents, etc. 


ARITHMETIC RULES 


A basic feature of floating point arithmetic is its prop- 
erty of automatic adjustment for overflow at the left. 
This feature is independent of the use of normalized 
form, however. The rules to be given apply to unnormal- 
ized operands, and have the property that the results 
give some indication of their expected precision. Besides 
this, the rules have certain advantages from the point 
of view of operation speed. 


* Manuscript received by the PGEC, April 17, 1957; revised 
manuscript received, August 15, 1958. This work was supported by 
the U. S. Atomic Energy Commission. > } : 

+ Institute for Computer Research, University of Chicago, 
Chicago, Ill. 
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: Significant Digit Computer Arithmetic: 


N. METROPOLIS} anp R. L. ASHENHURST} 


Addition and subtraction 


In performing floating point addition or subtraction, 
it is necessary that one or the other operand be shifted, 
and its exponent adjusted accordingly, so that the two 
exponents agree and the numbers are, so to speak, 
“lined up.” In a normalized system, the number with 
the smaller exponent must be shifted right; otherwise 
digits would be lost at the left. It is conceivable that 
another procedure could be used for unnormalized 
operands for which some left shifting is allowed. How- 
ever, it turns out that the desired significance properties 
are obtained if numbers are treated much as in a nor- 
malized system, except that no further adjustment is 
made at the completion of an addition or subtraction 
where there is no overflow. It should be noted that 
overflow may be expected to occur much less often when 
operands are not required to be in normalized form. 


Multiplication 


There is no fundamental necessity to shift operands in 
the case of multiplication; the product can be formed 
directly and the exponents added. If only the leftmost 
n digits of the resulting 2-digit product are retained, 
however, significant digits may be lost. Hence a read- 
justed product is defined, which has the property that 
the number of significant digits in the -digit result is 
approximately the same as that for the factor with fewer 
significant digits. Conceptually, if that factor were 
multiplied by a value near 1, the desired effect would be 
automatically obtained. Hence the product can tenta- 
tively be defined as the result of multiplying the factor 
with the smaller fractional part by the other factor nor- 
malized, and adjusting the exponent accordingly. In 
this case, however, the high order part of the 2n-digit 
fractional part product has either the same number of 
significant digits as the unnormalized factor or else one 
less significant digit. Thus, continued use of the rule in 
this form would introduce a tendency toward loss of 
significant digits. To counteract this, it is advisable to 
incorporate an additional adjustment of the product, 
consisting of a corrective left shift of one place, with a 
corresponding change of exponent, in certain circum- 
stances as specified by an auxiliary rule. Perhaps the 
most obvious way of formulating this auxiliary rule is to 
have a corrective shift called for when and only when 
the number of leading zeros would otherwise increase 
by 1. There are other possibilities, however; one is to 
make the discrimination on the basis of the normalized 
factor instead of the product. In this case the corrective 
shift is ordered whenever the fractional part magnitude 
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of the normalized factor is less than some specified value, 
say 3. (Although this rule can give a product with one 
less or one more leading zero than the smaller factor, its 
average effect is approximately that desired. Of course, 
no corrective shift is applied if overflow would result.) 
The comparative advantages of these two criteria for 
corrective shift involve some subtlety, and a detailed 
analysis is needed to choose between them. Both possi- 
bilities are considered below, in the discussion of imple- 
mentation. 

If one factor has fractional part 0, then the product of 
course has fractional part 0, but the exponent should 
reflect the fact that this value may be only an approxi- 
mation. It is natural, and turns out to be justifiable, to 
subtract from the product exponent just that amount 
necessary to normalize the nonzero factor, with possibly 
an additional unit decrement for correction. If the other 
factor has fractional part 0 also, an attempt to normal- 
ize it would result in its exponent being decreased by a 
before the fact became evident. If this is taken to define 
normalization for zero fractional parts, this case re- 
quires no special treatment. 

Example (with 2=8): let two numbers to be multi- 
plied be represented (fractional part and exponent) by 


fi = — .01100010 é: = — 100 

fo = + .00010010 eg = + 001 
Since |f:| >|f2|, multiplication is defined in terms of f; 
normalized : 

fi’ = — .11000100 ey) = 


I 


— 101. 
The 16-digit product fi’fo, with exponent e;’+é, is then: 
fifa = — .00001101 11001000 e+ e, = — 100. 


If corrective shift is ordered only for|fi'| <#, then this is 
the final form (without round-off). If corrective shift is 
ordered whenever the product (upper 8 digits) has one 
more leading zero than the smaller operand fz, then the 
shift would be called for in this case. The final prod- 
uct as delivered would be: 


fs = — .00011011 é3 = — 101. 


Division 

A criterion analogous to that described for multiplica- 
tion is used in defining the division operation; the first 
step is to cause the quotient to be adjusted so that it has 
approximately the same number of significant digits as 
the operand with the smaller fractional part. If the 
numerator has the smaller fractional part, the desired 
effect could be achieved by dividing by the denominator 
normalized, which is close to 1. If the denominator has 
the smaller fractional part, it is most simple to conceive 
of the numerator being first adjusted to the right so that 
it has the same number of leading zeros as the denomi- 
nator, and then being divided by the normalized de- 
nominator as before. In either case, the quotient has the 
same number of significant digits as the operand with 
the fewer such digits, or else one more. Which alterna- 
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tive obtains depends on the relative magnitudes of 
numerator and denominator, if both are considered 
normalized. In conventional normalized arithmetic, the 
extra digit represents overflow, and adjustment must 
be made. Here it may simply be considered a legitimate 
part of the quotient unless the operands happen to be 
in normalized form to begin with. However, as for 
multiplication, it is advisable to correct the basic result 
for the tendency to change (here, increase) in number 
of significant digits. This correction most naturally 
takes the form of an extra right shift and exponent ad- 
justment whenever the extra digit occurs. The correction 
then also effects adjustment for actual overflow. 

The quotient exponent is basically defined as the dif- 
ference of the exponents of numerator and denominator. 
This result must be further adjusted to be consistent 
with the effective quotient shift. The case of one or both 
operands 0 is given no special treatment. The result has 
fractional part 0 and the exponent which is naturally 
obtained when the above described procedure is carried 
out. 

Example (with z =8): let numerator and denominator 
be represented by: 


fi: = — .01100010 1 = 
f2 = + .00010010 2 = 


Since lf = | fel , consider f; to be adjusted to the right 
(but without discarding any digits): 


fi = — 00011000 10 ey’ = 


— 100 
+ 001. 


I 


— 010 


The 8-digit quotient (remainder ignored), with exponent 
é1’—é:’ (for normalized denominator, f’), is then 
— = — 00101011 ex’ — es’ = + 000. 
Since the quotient has one less leading zero than the 


smaller operand f2, corrective shift is called for, and the 
quotient is delivered as: 


f; = — .00010101 é3 = + 001. 


IMPLEMENTATION 


In this section the realization of the significant digit 
operations in a parallel binary arithmetic unit is dis- 
cussed. It is proposed that addition and subtraction be 
carried out more or less as outlined, by comparing ex- 
ponents, shifting the fractional ‘part with smaller ex- 
ponent to the right until the exponents match, then 
performing the addition or subtraction. If overflow oc- 
curs, the customary adjustment is made. Since no 
further shifting is performed, even when leading zeros 
appear in the result, the average addition should be less 
time-consuming than the conventional normalized one. 
Furthermore, it seems likely that with use of the pro- 
posed system there will be a tendency toward less dis- 
crepancy among exponents than under the normaliza- 
tion constraint, and this should also decrease the average 
addition time. 
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_ Multiplication and division are both conventionally 
performed in a series of 7 steps, counted in the so-called 
operations counter. At each step a shift is effected, pos- 
sibly accompanied by an addition (or subtraction) of 
the multiplicand (or divisor). The pattern of such addi- 
tions or subtractions is determined by the successive 
igits of the multiplier (in multiplication) or the suc- 
cessive signs of the partial remainders (in division). 

In the multiplication process, it is necessary in some 
way to obtain the 2n-digit product shifted to the left 
by an amount equal to the number of leading zeros in 
the factor with the larger fractional part; then a cor- 
rection, consisting of a left shift of one additional place, 
is applied if circumstances warrant. The basic shift can 
be accomplished in advance, by normalizing the relevant 
factor, and then proceeding with a standard n-place 
multiplication. This has the advantage of procedural 
convenience, if the correction rule which requires in- 
spection of this factor in normalized form is used. The 
digits to be examined are then in the leftmost stages of 
the register in which the factor resides just before the 
multiplication proper begins. The time taken, however, 
is the u-step multiplication time plus the shifting time 
necessary to effect the normalization. The latter time 
can be saved by doing the multiplication process in a 
slightly unorthodox manner: let the factor to be nor- 
malized be the multiplicand. Then, instead of starting 
by shifting the partial product right at each step, shift 
the multiplicand Jeft until it is normalized. At this point 
the situation is just as it would have been had the same 
number of steps been completed with the multiplicand 
normalized in advance, using the conventional scheme. 
‘The multiplication can then proceed in standard fashion 
from this point on, and the leading digits of the multi- 
plicand can be used to determine whether a correction 
shift should be performed at the finish (which actually 
amounts to omitting the last right shift of the multipli- 
cation, unless overflow adjustment is required). 

. Unfortunately, the short-cut described above cannot 
be employed if the other correction rule is used (the rule 
which calls for the corrective shift each time the product 
has one more leading zero than the smaller factor). Un- 
less the multiplicand is normalized in advance, the digit 
to be inspected may never appear in the leftmost stage 
of the partial product register, and thus remains vir- 
tually inaccessible. 

- The factor selected as multiplier in the above de- 
scribed processes is that with the greater number of 
leading zeros, and hence it may be expected that the 
average number of steps requiring an addition operation 


(as opposed to a shift only) is minimized. Although cus- 
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tomarily this time-saving advantage has been foregone 
when complement representation for negative numbers 
is used, because leading zeros become instead leading 
ones, a multiplication scheme which is symmetric in 
this respect can be developed, where an addition is 
ordered whenever the sequence of digits has a change 
from 1 to 0, a subtraction is ordered when the change is 
from 0 to 1, and neither is required when there is no 
change. 

Division according to the proposed rules can be ac- 
complished quite naturally in 2 steps. Here the most 
convenient way to proceed is to start by shifting left 
both numerator and denominator, adding 1 to the 
quotient exponent and counting in the operations 
counter (equivalent to adding a leading zero to the 
quotient). When one of the two operands is normalized, 
the procedure changes; if the denominator is normalized, 
shifting of the numerator continues, with counting in 
the operations counter but no change in the quotient 
exponent. This is continued until the numerator is nor- 
malized. Alternatively, if the numerator is normalized 
first, shifting of the denominator continues, with the 
exponent increased by 2 at each step, and counting in 
the operations counter, until the denominator is nor- 
malized. When both operands are normalized, the di- 
vision proper begins; it continues for only that number 
of steps necessary for the operations counter to reach 
the full count of z. Note that if the numerator is larger 
than the denominator after both operands are nor- 
malized, the first quotient digit is not an overflow as it 
would be in normalized arithmetic, unless no pre- 
liminary shifting has taken place. However, as men- 
tioned previously, an extra shift right should be per- 
formed in this case anyway as a correction. A scheme for 
speeding up the process, by ordering an arithmetic oper- 
ation (subtraction or addition) only when the sequence 
of digits in the partial remainder changes from 0 to 1 
or from 1 to 0, is available for division as it is for multi- 
plication. 

The case of zero operands is included automatically 
in these implementation schemes if the attempt to nor- 
malize is always halted after m steps. 


CONCLUSION 


Significant digit arithmetic seems to offer a means of 
analyzing floating point calculations more effectively, 
and also to possess certain advantages from the equip- 
ment standpoint. It is planned to use this system on the 
Maniac III computer, now under construction at the 
Institute for Computer Research of the University of 
Chicago. 
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Minimal “Sum of Products of Sums” Expressions of | 


Boolean Functions” 
SHREERAM ABHYANKART 


Summary—The problem of economical synthesis of circuits for 
digital computers leads to the problem of finding Boolean expressions 
of minimal length equivalent to a Boolean expression f. Previous 
authors restricted themselves to ‘‘sum of products” expressions; 
dualizing this gives ‘‘products of sums.’? The next more efficient 
step is to find minimal ‘‘sums of products of sums’ expressions. 
In this paper, the basic concepts are formulated in Part I, and general 
theorems are given in Part II. In Part III, all the distinct minimal 
“sum of products of sums’ expressions are obtained in case the 
cell complex of f consists of two isolated points. For the case of 
three isolated points partial results have been obtained which will 
be published in a later communication. 


I. INTRODUCTION 
Section 1 


“YE START by fixing our notations. By a Boolean 

expression f(x1, %2,°°:*, Xn) in the letters 

X1, Xe, °* +, xX, will be meant an expression in 

1, X2, °° *, %, and their complements x1’, x2’, - - + , Xn’ 

joined by a certain number of product and sum opera- 

tions. We can render this more precise as foilows: a 

literal, z.e., an expression either of type x; or of type x,’, 

will be called a Boolean expression of zero iteration; we 

shall follow the usual custom of denoting x; and x,’ by 
x;1 and x;°, respectively. Expressions of type 


IV 


mT ik ji Se jm 
Il Vi, = Wiig? © * Vimy m 
k=1 

and of type 


™ . ° ° * 

Ik I 72 Im 
DF ti, =H, tet te, mel, 
k=1 


will be called Boolean expressions of one iteration, re- 


spectively, of type product (in short p) and of type sum 
(in short s); here z and 7, take values in the sets 
(1, 2,---,) and (0, 1), respectively. By induction we 
define a Boolean expression of (¢+-1) iterations of type 
product-sum-product- - ++ (in short psp ---+ ) to be an 
expression of type 


IT ge(, - - 
k=1 


= gi(%1, es 


5 An) 


"5 Wn) So(w1, +++ Bn) + * + Bm(1 +++, Kn), 


* Manuscript received by the PGEC, December 17, 1957; re- 
vised manuscript received, May 20, 1958. First part of a report en- 
titled “Investigation of Mathematical Methods for the Analysis 
and Synthesis of Computer Circuits,” submitted to AF Cambridge 
Res. Center, under Contract No. AF 19(604)-1818, Columbia Uni- 
versity, New York, N. Y., September, 1956. 

{ Princeton University, Princeton, N. J. 


and a Boolean expression of (q¢+1) iterations of type — 


sum-product-sum- ++ - (in short sps +--+) to be an ex- 


pression of type 


Dy, gK(X1, Lp Sade Xn) 
k=1 

= 91(%1,° °° 5 Sn) - 22+ gm(%1, °° > Xn)y 
where g)(%1, °° +>, %n))° °°» Zm(%1,°°°, ¥n) With m21 


are Boolean expressions of qg-iterations of types sps :: - 
and psp-+:++, respectively. A Boolean expression 
f(x1, °° +, %n) is then a Boolean expression of some 
iterations. Let f(x1, - - - , Xn) be a Boolean expression of 
g-iterations. Observe, that if g is even then f is either of 


type 


sp sp--+-sp(“sp” repeated g/2 times)---, (1) 
or of type 
ps ps +++ ps (“ps” repeated g/2 times)--+3; (2) 
and if g is odd then f is either of type 
sp sp-+-++sps(“sp” repeated (¢g — 1)/2 times), (3) 
or of type 
ps ps +++ ps p(“ps” repeated (q — 1)/2 times). (4) 


If h is any positive even integer then f can be considered 
to be a Boolean expression of g+h iterations, respec- 
tively, of type 


sp sp +++ sp (“sp” repeated (q + h)/2 times), (1*) 
ps ps +++ ps (“ps” repeated (¢/h)/2 times), (2*) 
sp sp +++ sp s(“sp” repeated (¢ + kh — 1)/2 times), (3*) 
or 


ps ps +++ ps p (“ps” repeated (¢ + h — 1)/2 times), (4*) 


according as f is of type (1), (2), (3), or (4); similarly, 
if is any positive odd integer then f can be considered 
to be of g+h iterations of a corresponding type. 

Note that in Boolean expressions, addition and multi- 
plication are to be commutative, thus, for instance, 


XyX2 + x34’ = x4'x3 + HK. 


By a general Boolean expression of 0 or 1 iteration 
will be meant a Boolean expression of 0 or 1 iteration, 
respectively; by induction we define a general Boolean 
expression of g+1 iterations to be an expression of type 


fi(a1, rs Ln) pe ee + fin (21, 5 


ae Sn); 


3 
: 
j 


ie te a 
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¢ of type 
file, “ee Xn) fo(x1, cade a) ac PACE rey als 
yr of type 
: [f(o, as, 02"), Ma), 
vhere fi, fo, > ++, fm,f are general Boolean expressions of 


iterations. By the length of a general Boolean expres- 
ion f(x1, + + +, %n), to be denoted by L(f(m, +--+, xn)), 
will be meant the number of sum and product operations 
i fei, > + +, a). 

Let f(x, +++, %») be a general Boolean expression. 
[t naturally gives rise to a Boolean function which we 
shall denote by f(m,--+, xn) of the binary variables 
K1, X2, ° * * , X, in the usual sense; when there is no con- 
‘usion, we shall replace Migs 2.) iY f(Gj0> op an) 
where (v1, - ++, ¢,) is a vector with components 0 or 1. 
Let g(x, -- +, Xn) be another general Boolean expres- 
sion. We shall say that the Boolean expressions 


F(x1, ++ +, Xn) and g(x, +++, Xn) are equivalent, in sym- 
Dols f(x1, + + +» Xn) ~g(x1, ++ + , Xn), if the corresponding 
Boolean functions f(x:,--+, *n) and 2(x,---, Xp) 
are equal. 

Informally speaking, f(m1, +--+, *n)~g(x1,-+ +, Xn) 


if, and only if, f can be converted into g by a finite 
number of applications of the various Boolean laws. 

_ Given a general Boolean expression f(“1, + - + , Xn), by 
applications of the laws of complementation we can at 
once find a Boolean expression F(x, x2, +++, Xn), such 
eat f(%1, -* °; Xn) ~F(%, x2, °°: , Xn), and we shall 
have L(f) =L(F). Hence, without any loss of generality 
we may, and henceforth we shall, confine our attention 
to what we have called Boolean expressions (and not to 
the more general general Boolean expressions). 

_ Given a Boolean expression f(x1, %2,- ++, Xn), the 
Boolean expression obtained from f(x1, %2, ---, Xn) by 
interchanging x; and x,’ for 1=1,---, and inter- 
changing sum and product will be called the Boolean ex- 
pression complementary to f(x1,--+, xn) and will be 
denoted by 7f(x1, - + + , Xn). 

_ Now we shall illustrate the various concepts intro- 
duced in this section by taking examples in five vari- 
ables. 


(A): 

1) x, x3’ are Boolean expressions of zero iteration. 

2) xi+x4, x3+x1+%5 are Boolean expressions of 1 
iteration and of type s. 

3) x3%4(x%1-+%2’) is a Boolean expression of 2 iterations 
and of type ps. It can also be considered as a 
) Boolean expression of 3 iterations and type sps or 

of type psp; also as a Boolean expression of 4 
iterations and of type psps; etc. 

A) (21% +302" 25) (201%4%5+%3") (2 +-%3-+24) (X3'K5H1’) Is a 

Boolean expression of 3 iterations and of type psp. 

5) [ (ocr + x3 + x4) (%2 + TOOTS os eae es x4) x51" | 

[xst+-x,] is a Boolean expression of 4 iterations 
and type psps. 
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(B): 


(x3 +-21)' +-X4%5+ (x1%2)’ is a general Boolean expression 
of length 5 and it is equivalent to the Boolean expression 
x3/ Hy" -%4%5-+21'+-x2’ of the same length which is ob- 
tained from the original one by the use of the law of 
complementation, 7.e., setting: (x3-+%1)'’=x3/x,' and 
(x 1%)’ = x4’ +x’. 


CC) s 
Let 
f(41, %2, %3, 4, 5) = %1(%2 + 43’) 
and 
(x1, H2, Ha, Xa, Xs) = XyX2 + x43", 
Then 


F(1).000).0).0) SAO: 0) 100 pete 
and 
3(1,.0,.05050) = 20 ot S00 de 1 Oe 


Similarly for all the 25=32 vectors it can be verified 
that the value of f equals the value of g. Hence f=Z so 
that f~g. This can be seen by applying the distributive 
law to f and thus obtaining g. Note that L(f) =2 and 
L(g) =3. 


(D): 
Let 
f(%1, %2, %3, L4, Ly) = Lixo + x3(%1 + a2’) + w5/(x3 + 22’). 
Then 


th (a1, X2, %3, V4, X5) = (co + 42) (as! + 4X1’ 2) (x5 Se 43’ %2) . 


Section 2 


For a given Boolean expression f(x, %2, +--+, %n) we 
shall denote by Y(f) the set of all the Boolean expres- 
sions in 1, X%2, °°, %, which are equivalent to f; by 
Ysp(f), Yps(f) and Ysps(f) we shall denote the sets of 
all the Boolean expressions in x1, x2, + + + , %», which are 
equivalent to f, and which are of type sum-product, 
product-sum and sum-product-sum, respectively; by 
LZ(f), LZsp(f), LZps(f) and LZsps(f) we shall denote 
the minimum of L(g) for gin Y(f), Ysp(f), Yps(f), and 
Ysps(f), respectively; finally by Z(f), Zsp(f), Zps(f) and 
Zsps(f) we shall denote the set of members of Y(f), 
Ysp(f), Yps(f) and Ysps(f), which are of length LZ(f), 
LZsp(f), LZps(f) and LZsps(f), respectively. Observe 
that 


Vif) Vsps(f), Vsps(f) D Ysp), Ysps) Ses); 
and 
LZ(f) S LZsps(f) S min [LZsp(f), LZps(f)|. 


Also note that since there are only a finite number of 
distinct Boolean expressions of a given length, Z(f), 
Zsps(f), Zsp(f) and Zsp(f) are finite sets. 
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The importance of the problem of finding Boolean 
expressions g(1, : + * , X,) equivalent to a given Boolean 
expression f(x), %2, * + -, Xn) such that L(g) is minimal; 
i.e., the problem of finding Z(f) and calculating LZ(f) 
is due to its interpretations as the problem of simplifying 
truth functions and as the problem of economical synthesis 
of circuits for digital computers. But concerning this ab- 
solute minimization problem, due to its very difficult 
nature, up until now nothing has been done. Various 
authors have exclusively restricted their attention to 
the much simpler problem of finding minimal normal 
expressions, 7.e., in our terminology that of finding 
Zsp(f); among these authors prominent have been 
Quine,!’? Urbano and Mueller,’ and Petrick.4 The 
methods developed by them for finding Zsp(f) are quite 
good, and we at once observe that by applying their 
methods to the complement rf of the given expression 
f, one can find Zps(f). Zps(f) coincides with the set of 
complements of members of Zsp(rf). As we shall pres- 
ently illustrate, even with expressions in as few as four 
letters, we may have LZ(f) <min [LZsp(f), Lps(f)]. We 
take the example considered by Mueller and Urbano? on 
page 129 but we label our letters as x1, x2, x3, x4 instead 
Onsen, Coy oe Oe 


Pr enmity bel <p eeel fe nian 
f(%1, %2, %3, %4) = Hy'He'H9! a4 Eye! 3) Ha" PF Les 4 
ey, / f ies / ! , 
ate tM Moiese ao beta 1 Loved 
+ 1X23 v4 + 1 XoxeX, + LiMewets. 


The elements of Zsp(f) given on page 130 of Urbano 
and Mueller have length 10, z.e., we have 


LZsp(f) = 10. 
To find Zps(f) we proceed by observing: 
Th (41, %2, 3, 4) 
oy f' (a1, %2, %3, 4) = Hy'H2/xgx4 + 41/X0!x3/H4 
a9 44' X9x3/ 4" ae Hy’ Lo%3" Ha am 4199" 90324 


qe X1%9/ 3/4 = XX 9H 304". 


Let us plot f’ on a four cube. (See Fig. 1.) Now each of . 


the basic cells I, II, III covers a vertex not covered by 
any other basic cell, and together they cover the entire 
complex.' Hence, there is only one irredundant normal 
form: g(%1, X2, %3, Xa) = X1%2%304' +2 x2%3'-+x9'%4 which is 
then the only member of Zsp(f’). Hence, 


1W. V. Quine, “The problems of simplifying truth functions,” 
Amer. Math. Monthly, vol. 59, pp. 521-531; October, 1952. 

2W. V. Quine, “A way to simplify truth functions,” Amer. 
Math. Monthly, vol. 62, pp. 627-631; November, 1955. 

3. R. H, Urbano and R. K. Mueller, “A topological method for 
the determination of the minimal forms of a Boolean function,” 
IRE Trans. ON ELECTRONIC CompuTERS, vol. EC-5, pp. 126-132; 
September, 1956. 

4S. R. Petrick, “A direct determination of the irredundant forms 
of a Boolean function from the set of prime implicants,” Computer 
Lab. A.F.C.R.C. Bedford, Mass.; April, 1956. 

5 We shall use, without explicit reference, the geometric notions 
and terminology introduced by Urbano and Mueller, footnote 3. 
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T2(%1, X25 V3, 44) 


= 41" + ae! + x3! + 94) (41 + 2! + 4) (42 + 24). 


is the only member of .Zps(f) and LZps(f) =L(rg) =8. 
Thus, 


LZsp(f) = LZps(tf) = 10 and LZps(f) = LZsp(zf) = 8. 
Now let 
k( a1, %2, %3, %a4) 
= [ae! + (ai! + 2x3! + x4) (41 + x3) |(a2 + 44’). 

Then k~rg~f and L(k) =7. Therefore, 

LZ(f) S LZpsps(f) $7 < min [LZsp(f), LZps(f)] = 8, 
and 
LZsps(f) S$ min [LZsp(f), LZps(f)| = 8 < LZsp(f) = 10. 


It is clear that when we are dealing with Boolean ex- 
pressions in a large number of letters, LZ(f) will in 
general be much smaller than Zsp(f) as well as Zps(f). 


nua ey xf Xo x3 as xy X X3 x} 


(0, 0,0, 0) 


Lett les] 


Section 3 


Thus, the minimizing efficiency yielded by Zsps(f) is 
at least as much as the minimizing efficiencies of Zsp(f) 
and Zps(f) put together. Of course the most efficient 
minimizing will consist in developing methods for find- 
ing the absolute minimal expressions, z.e. Z(f). As we 
have already stated, nothing has yet been done in either 
of these two directions; hence any sort of beginning, 
howsoever modest, is a worthwhile goal. To contribute 
to the study of Zsps(f) is the purpose of the present 
paper, while a contribution to the much more difficult 
study of Z(f) will be made in a paper which will follow 
this one.® 


Section 4 


As in Urbano and Mueller’ with a Boolean expression 
in 2 letters we may associate its cell complex in the u- 
cube.’ It is reasonable to start by studying those cell 


_ ©S. Abhyankar, “Absolute minimal expressions of Boolean func- 
tions,” in preparation. 
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complexes which possess nice geometrical properties, 
e.g. purity of dimension (7.e. all the basic cells® being of a 
same dimension), connectedness, etc.; first because 
natural problems are to be expected to yield configura- 
tions with nice properties, and second, by decomposing 
a given cell complex into nice components, one can use 


the theory of special configurations towards the minimi- ’ 


zation of the general case. We propose to start with pure 
zero dimensional complexes, 7.e., complexes which con- 
sist of a certain number of isolated points. Another justi- 
fication for studying low dimensional (dimension of a 
complex=the maximum dimension of a basic cell oc- 
curring in it)—in particular zero dimensional-complexes 
is that in a random given set of vertices the probability 
of the existence of an s-dimensional cell drops rapidly 
as S increases. 

In this paper we shall study complexes which consist 
of one or two isolated points in an arbitrary dimensional 
cube (7.e., no restriction on the number of letters of the 
Boolean expressions under consideration). We remark 
that the importance of our present investigations per- 
haps lies more in the methods used in the proofs than 
in the results themselves. Throughout the paper, we 


illustrate the results by examples in a small number of 


variables. 
II. GENERALITIES 
Section 5 
Let f(%, - - +, %n) be a Boolean expression. We shall 


denote the cell complex (on the z-cube S,) of 


| f(x, are 
_ by C(f) and the point set consisting of the points of C(f) 


_ equivalent to f(x, - 
the 2” vectors (i, %2, °°: 
Consequently, a statement about the point-complex 


+, Xn) (as considered by Urbano and Mueller®) 


by D(f), and we shall call D(f) the point of complex of f. 
Recall that the points of D(f) are in one to one cor- 
respondence with the terms of the canonical Boolean 
expression 


F(a1, °°, Xn) = >, fis, ta, + + +> 5 in) aria’? + > + ania 


- +, Xn), (the sum is taken over all 
, In) With coordinates 0 or 1). 


_ D(f) is simply the geometric equivalent of an arithmetic 
_ statement about F, and hence we can and will use these 


_ two modes of studying f(m, - - - 
_ For a point set E of S, , we shall denote by 7(S,)E the 


, Xn) interchangeably. 


_ complement of E in S,, when the reference to S, is clear 


be a Boolean expression where gi, - - - 
_ expressions of type ps: - 


from the context r(S,)E will be replaced by 7E£. 


Let 


(x1, < gees) Xn) = DS ge (oa, aE) Xn) 
k=1 


, £m are Boolean 
-,. Then 


D(g) = U D(g:) and g~f 
k=1 
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if, and only if, D(g) =D(f). 

Now assume that gi, - + - , gm are of type product-sum 
so that g is of type sum-product-sum; let E,=D(g;). 
Then we have that g € Ysps(f) if, and only if, 


mm 


U FE, = Dif); and g © Zsps(f) 


k=1 


if, and only if, 
UE, = D(f) and L(g) =m—1+ D7 L(g) 
k=l k=l 


is a minimum, 2.e., if, and only if, for k=1,---+, m, gx 
is in Zps (Ex), 


UE = 


k=1 


D(f) and m—1+ 3 LZ ps( Ex) 


is a minimum; furthermore when the latter is a mini- 
mum and g;’ is any member of Zps(F;x), then 


(Ze) 


Sip ha) 


k=1 


=m —1-+ D)LZps(Ex) = L(g) = LZsps(f), 

k= 
and hence >0%, g:’ is also in Zsps(f). Now LZps(Ex) 
= LZsp(rEx). Therefore, the problem of finding Zsps(f) 
consists in finding partitions &,,---, En» of D(f), for 


which 
m—1-+ >> LZsp(rEy) 
k=1 
is minimum, and finding Zsp(rE)), -- +, Zsp(rEm); let 
Ey, +--+, Em be any such partition of D(f) and let be 


any member of Zsp(rE;), then 07%, rh is in Zsps(f), 
and all the members of Zsps(f) are obtained in this 
manner. 


Section 6 


We introduce the following notation. Consider a 
Boolean expression f(x1, ++ >, Xn) in letters x1, °° +, Xn 
which is of the product type, 7.e., f(%1, °°, Xn) is the 
product of a certain number of x; and a certain number 
of x;’. Let S, be the m-cube of the letters x1, - ++, Xn. 
Then D(f) consists of the points of a certain subcube 
Sm (4.€., a cell) of S, (the subscripts of S,, and S, denote 
their dimensions). By I(f) or I(Sm) we shall denote 
the set of subscripts 7 for which x; occurs in f and 
x,’ does not; by J(f) or J(Sm) we shall denote the 
set of subscripts 7 for which x,’ occurs in f and x; does 
not. Finally by (1, 2,--- , 2) Gf) or KU, +33. ass) 
we shall denote the complement of J(f)UJ(f) in 
(1, 2,- +--+, 2); when the reference to (1, 2,---, ™) is 
clear from the context we shall use the shorter notations 
KG )ror KS.) 
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Section 7 


Finally we are ready to prove a theorem. Let 
f(x1, +++, %n) be a Boolean expression in letters 
1, °° * , *, and consider the point-complex D(f) of f in 
the n-cube S,. Let Sm be an m-cell of S, which contains 
D(f); in particular S, could be the subcube of S, 
spanned by the points of D(f), z.e., the smallest cell of 


S, containing D(f). Let 
(ha, Dice ak Rm) pai KY Deen, 1) (Sim). 


Then as a subset of S, we can think of D(f) as giving 
a canonical Boolean expression F(%i,, Xk, °° * » km) in 


thems: letters ve,, °°, vx. Let 
A (x1, Ce te %n) Sa iT Xt II X;'5 
1EI (Sm) JET (Sm) 


SOrenata (ives os, Xy,) GENNES Sill Sin 
Theorem 1: Let G(xx,, - - * , Xtm) be a Boolean expres- 
Siomin ZPs( Hi Xpse icy Xk.) 2) Hen 


g(%1, See aS, an) = G(x, Re ene MH (Xi; os 


is in Zps(f(x%, - - + , x,)). Furthermore Gg is a one to 
one correspondence between the members of Zps(/) and 
members of Zps(f). 

Proof: We may relabel the letters so that 


5 Wn) 


TAS eee Ay 2 ck St) 5 
T(Sm) = (m+ 1,m+2,---,m+s) 
and 
J(Sm) = (m+sti1i,m+s+1,---,m+s+i=n). 


Let A=7(S,)(S,) and B=7(S,,)(D(f)) =7(Smn)(D(P)). 
Then A and B are disjoint and 


D(zf) = 7(Sn)(DGf)) = AUB. 


Let X; and X;,’ be the (n—1) cells of S, with equations 
*,=1 and x;=0, respectively. Then X41’, Xmio', + --, 
NG pe On Pals Gris oo sia care. COntaimed. in oa. 
and hence they are all basic cells of C(zf). We are 


assuming that f is not inconsistent, 7.e., D(f) is not 


empty so that m>0, and hence S, is not in C(rf). Now 
the point of S, given by 


ty = H2 = + = hm = 1, emp = 


= Xmpe = 1, Xmtep1 = Lmtep2 = 


(z.e., the point which corresponds to x - + + XmXm41’, 
Nn taarg cir = ume h ae ae les My Anis DUS GhAn 
any of the remaining above mentioned basic cells of 
C(7f), and hence Xm41’ forms an essential star® of C(zf). 
Similarly, each of the above mentioned basic cells of 
C(rf) is an essential star of C(7f). Hence, in view of the 
results of Urbano and Mueller,’ and Petrick,* the mini- 
mal coverings of D(rf) consist of the minimal coverings 
Of Do AdjOmed ab ye t0e, COllse\ ny tai miaig aie Y ated 
Xm+st1) °° *, 4n. Arithmetizing this geometrical state- 
ment we conclude the following. If G(x, +++, %m) isa 
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member of Zsp(7rF), then 
see) oe) = G(x1, CR" ee ) Xm) + mis! -F tmpe! ig ce 
as Mes: “F Gobet ttc ete tsie eet ee 


is a member of Zsp(7f), and Gog is a one to one cor- 
respondence between the members of Zsp(rF) and the 
members of Zsp(rf). Complementing this result we at 
once obtain the required one to one correspondence be- 
tween Zps(F) and Zps(f). 

Example: Take n=6 and let 


g(x, " 


f(*1, C2, X3, La, V5, x6) 
= 25/%1%3'%6 + ax2%3/X%— + x3 (X4x6 + X3Xs). 


Then applying the distributive and cancellation laws we 
have 


fer, 5, oy XA, Loy Xe) 
~ 41%3/X5'%— + xovs’X— + X3/Xaxg + Xs XgX5 


/ / / i: 
~ 103 X5 Xe + XoX3 X— + Xz sXe. 


In each term of the last expression, x3’ and xs are com- 
mon. This means that the point-complex D(f) of f in the 
6-cube Ss is contained in the 4-cube S,; represented by 
the expression H(x1, x2, %3, %4, %5, Xs) =%X3'xe. We have: 
I(S:) = (6), J(Ss) = (3), KC, 2, 3, 4, 5, 6)CSs) = Ci, 2, 4, 5) 
Factoring out x3’x— from f we set F(x1, x2, X4, Xs) 
X1X5’-+X2-+x—. Then the point-complex of F considered 
as a subset of S, coincides with the point-complex of f 
considered as a subset of Ss. Theorem 1 now says that if 
G(x1, %2, %4, Xs) is a minimal ps expression equivalent to 
F (41, 2, %4, Xs), then g(*i, %2, %3, %4, %5, Xe) =GUri, te 
X4, X5)X3'%—s is a minimal ps expression equivalent to 
f(%1, %2, %3, X4, Xs, X6), and each such expression is ob- 
tained in this manner. 

Theorem 2: For the Boolean expression f(x, - + + , Xn) 
=Xi, °° * Xi, %;,! + + + x;,', where the subscripts 4, ---, 
1s, ly °° * Je are all distinct [m=n—s—t=0, so that f 
represents an m dimensional cell of S,; in particular if 
m=O then it is a point |, we have that f is the only mem- 
ber of Zps(f) and LZps(f) =n—m—1. 

Proof: In Theorem 1 take S,, to be the cell defined 
by f. 


Theorem 3; For the Boolean expression f(x1, + + + , Xn) 
=Xi, + + + Key - %j,' 9 + > 23, where the subscripts 
41, °:*, ts) Jus ° °°, Je are all distinct, and st+t=n, so 


that D(f) is a single point; we have that f is the only 
member of Zsps(f), so that Zsps(f)=Zps(f) and 
LZsps(f) =LZps(f) =n—1. 

Proof: Let g=gitge+ +++ +gq be a member of 
Zsps(f) where gi, g2, +++, gq are of type ps. Then 


qd 
U D(gi) = D(f). 
k=l 
Since D(f) is a single point we must have D(g,) = D(f) 
for some k, say D(g:) = D(f). The minimality of g implies 
that g=g, so that g=1 and g is of type ps. Now Theorem 
3 follows from Theorem 2. 


[ieee 


1958 


III. CAsE oF Two IsoLaTED Points 
Section 8 


Consider the Boolean expression 


#(%1, °° +, Xn) = (xe + > + tn + x1'xe! +++ a')’, n> 1. 


In this section we propose to find the minimal sum- 
product equivalents of f, 7.e., the members of Zsp(f). 
As usual, let C(f) be the cell complex of f in the n-cube S, 

For 747 let A;; denote the (~—2)-cell of S, defined 
by «:x;’, t.e., by: x;=1 and x;=0. Now xix - - 
the product of all the m unprimed letters and 
x1'%2" + + + x,’ is the product of all the ” primed letters, 
and hence neither is contained in the A,;;=x,x,’ which 
has one unprimed letter and one primed letter. There- 
fore, Ai; is a cell of C(f). Now the only (n—1)-cells 


- Xn 1S 


which contain A;; are x; and x;’; x; contains the point 
XX + 
x!%0" 


-%, of tD(f) and x;’ contains the point 
- - x,’ of rD(f) so that neither x; nor x,’ is in C(f). 
Hence, A;; is a basic cell of C(f). 

Now let £ be any cell of C(f) . If I(£) were empty 
then the point x;’x2’ - + - x,’ of rD(f) would be in E£, and 
if J(Z) were empty then the point xx. - - -x, of 7D(f) 
would be in £, each of these implications is a contradic- 
tion, and hence neither J(£) nor J(£) is empty. Fix 7 
in I(£) and j in J(&), then ECA,;. Thus we have 


| proved 


Lemma 1: The n(n—1) cells Ai; with 147 (each of 
which is 7—2 dimensional) are exactly all the basic cells 


| of C(f). 


One consequence of Lemma 1 is that any irreducible 


_ covering of D(f) consists of a certain number of the cells 


} 


_Aj;. Let H be a set of ordered pairs (27) such that the 


cells A;; with (zj) in H form a covering of D(f). 
Lemma 2: There exists 7 with (7) in H. 
Proof: Let H; be the set of members of H of type (im), 


_ and let Hy be the set of the remaining members of H of 
type (zj) with An. 


Since f = (x1%2 + + + Xn+41'%_’ > + + xn)’, in the canoni- 


“cal expression equivalent to f occur all the possible 


products xx," - - 


- x,** such that not all the & are 


_ zero and not all the &; are one, 7.e., at least one of the 


kz is 0 and at least one of the kis 1; let us denote by T 


the set of all these products [T is essentially the same 


Was D(f) |. 


For each (27) in H we write 


50; ~mw x;/ (>, xP lake Sao Xs = Verne as Cece XR) 


DS aerttacgtt «+ ah + + a9 + + + abn (14), 


where A denotes omission and where the summations 
are extended over all the 2*-? combinations of k:=0 or 


| t-with #=1,2,---,4, 


2459, =, wm, Since 


f~ ey eke 


(4j7)€H 


7 This is simply a numbering of these summations. 
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if we collect all the terms in summations (1,;) for the 
various (77) in H, we must end up with exactly the 
terms in T. 

Assume, if possible, that there does not exist any j 
with (”j) in H. Then any term in Hy is of type (27) with 
in and j#¥n. For a term of this type, the summation 
(1:3) contains x)" - - - x, x, if, and only if, it con- 
tains x1 +--+ x,y» x,’. Since x1'%_’ - - - x,’ does not 
belong to 7, it cannot be in (1,;), and hence x;/x2’ - - - 
Xn—1’ X, cannot belong to (1;;). 

Thus x1'x%2’ + + + X,1’ X, does not come from any term 
in He, but it does belong to 7, and hence it must come 
from some term in fj, 2.e., a term of type (7m). But each 
product in (lin) contains x, as a primed letter, and 
hence %1/%2’ + + + X,-1’ X, cannot come from any term in 
Hi. This is a contradiction which proves our assertion. 

Lemma 3: For any 7 with 1S7S7, there exists j with 
(aj) in H. 

Proof: Due to the symmetry in the subscripts 
1, 2,--+-+,m it is enough to treat the case of any one 
particular value of 2, say «=n. But for 1=n Lemma 3 
reduces to Lemma 2. 

Lemma 4: For any j with 1577, there exists 7 with 
(77) in HZ. 

Proof: In view of the symmetry between the primed 
and the unprimed letters, Lemma 4 follows from 
Lemma 3. 

Thus Lemmas 3 and 4 imply that each of the 2n 
literals x1, > - +, Xn, %1', °° +, Xn’ occurs at least once in 
the expression capex x:x;’. Let h be the number of 
terms in H. Then 2h=2n, 1.e., 


Lemma 5: 
eo 
Lemma 6: Let 
g(a1, °° + Xn) = Kite" + xoxe’ + + + > + Xn atn! Xn. 
Then g~f. 


Proof: Since g is the sum of a certain number of cells 
Aij, we have D(g)CD(f), and hence it is enough to 
prove that D(f)CD(g). Let xittx*? - + + x," be a point 
of D(f). Then k is zero for some ¢ and it is 1 for 
some t. Let u be the minimum value of ¢ such that 
k,=1, and let v be the maximum value of ¢ such that 
t2uand kyu=kuyi= +--+: =k: =1. If v= than k:=0 so 
that xx? - - - x,*" isin x,x,’, and if y<m then ky 41=0 


so that xx tt--- x," ig im %X X41’. Therefore, 
xyx'? - - - x, is in D(g). 
Lemma 7: x 1X_" +x0%3/ + + HXniiXn’ tnx’ is in 


Zsp(f) and LZsp(f) =2n—1. 

Proof: Follows from Lemmas 5 and 6. 

Now assume H is such that )oGjex xix,’ is in Zsp(f). 
Then Lemmas 5 and 7 tell us that H contains n pairs. 
Lemmas 3 and 4 tell us that given 7 there exists j with 
(17) in H. Hence if we label elements of H as (jx) with 
k=1,2,--++,mthen (u, %,-+°-+,%,) and (ji, jo, * - * » Jn) 
are both rearrangements of (1, 2,+--, ). Obviously 
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H is completely characterized by the permutation 
ea Bee) 
Ginga geste 
We may relabel the subscripts k so that 


c eee, 
H* = 
vB 72 Par “In 
1.e. such that 4,.=k. 


If we decompose H into canonical cycles we have 
either 


Gila Ee semeal( Ri Rosas Bp) OF 
(ii) H* = (Riko + + + Rm) (Rmythmy2 °° * Rmtm*) * 

with m <n 
where ki, ke,:-+*, Rk, is a rearrangement of 
(1, 2,---,). In case we have (ii), let K be the sub- 


set of H consisting of the pairs (hike), (Rek3),--- , 
(Rm—1Rm)(RmRi) and let L be the set of the remaining 
pairs of H. Then the pairs of L are of type (k,k;) with 
4>m and j>m. Consider the point 

ae pyaar aes eee 


Then P is in D(f). However, 
Pe ata) Mt gia AN Cua fa, Mt, 


[PE = Vi Vko * 


and 
PED ete sia > mand Sj. mM, 
and hence 
PE a SS, a) and PE€& p( oe vat . 
(uv) EK (uv)EL 


This is a contradiction since 


p( > vat uv ( Ds =) 
(uv)EK (uv)EL 
a p( » a) = DG): 
(uv) EH 


Therefore case (ii) cannot occur, and hence we must 
have case (i), 2.e., H* is an n-cycle. 

Conversely given an arbitrary n-cycle (qig2 + + + dn), in 
view of the complete symmetry in the subscripts 
(1,2,---,n), Lemma 7 tells us that 


, , 
Maqihas  Vasag to 8 gta, te XonHoy 


is in Zsp(f). 
If a and 6 are two distinct permutations of 1, 2, -- - 
n, then for some 12, a(z)#b(z) so that 


’ 


%1%a(1) + XoXa2)’ a5 an eae Ate AON 


A H1%p (1) $+ + + anteny’. 
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‘There are (4—1)!=3!=6 permutations of 2, 3, 4. They 
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Given an n-cycle H*=(hi, he,+ ++, Rn) we may 
cyclically permute the k’s so as to have ki=1, and with - 
this stipulation the n-cycle H* will uniquely determine 


ho, Rs, - + + , Rny and (ho, k3, - + - , Rn) is then a rearrange- 
ment of (2, 3,- +--+, 2). Conversely given a permutation 
€ Spain 1 ) 
bad Boles 
of (2, 3,-++, m) it uniquely determines an -cycle 
(1. ko ks -> +R). Thus, there are .@—1)! (distinct 
n-cycles. 


We summarize the results of this section in the follow- 
ing theorem. 

Theorem 4: Consider the Boolean expression 
Ff Bay 09 + 5p Bm) ia = en oar a= ees ee 
There are N=(n—1)! different permutations a, ds, 


-+, @y of 1, 2,: +--+, ” which are n-cycles. For the 
n-cycle a; let 
(AOE Gar es) Xn) = LF OP RGAY ela X2%a,;(2)" a eS PRS 
Then gi, ge, - ++, gy are all the members of Zsp(f) and 


LZsp(f) =2n—1. One way of enumerating the different 
n-cycles is to enumerate the VN=(n—1)! permutations 


w= ( ee) = 1 es 
Ro Rs-*+Rn 
of 2,3,---,m,and let a; be the n-cycle (1 ke k3 + + - Rn), 
1.€., 
DS ORs REC RATR ee n 
on (| elhsa ) 

so that 
gi(X1, shea Xn) 

= Likh,: te Cagkkg he eieteg se Eh ee 


Example: Take n=4 so that 


f(s, %2, %3, 44) = (~12%3%4 + 44’ X9' x3/x4')". 
are 
‘ @ 3 3 ; @ 3; -) ; 
Te ES ke SE NOs as ee 
2.3en 2, 3, 
ba = ) ) bs a ( 
Che ae 4, 2, 


The corresponding 4-cycles are 
Cn (Ge 2, 3, 4), ph = is 2, a; 3); 
a4 = Cb 3, a 2), CS as a; 2, 3), 


dg =a, 3,-2.4), 
de = "(ty Ay 32)3 
There are 6 minimal sp expressions gi, go, 3, 4, 5) £6 


equivalent to f, their length is 7, and we can write them 
down, using the 4-cycles a;, in the following two ways: 


/ 
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= 419! + xox’ + agay! + aay’ 


— X 1X9" Se Xox3" se X34" Ss ati. 


81(%1, 2, X3, 4) 
go(H1, Ho, X3, X4) = LyXo! + wong’ + aang’ + x3! 
= 41%0' + xox," Si x3ay! apehaes 
= 1X3" a tebe az xa" Peay 


= 4103 + Xora! + wgxe’ + agay’. 


3(X1, Xe, 3, Xa) 


ga(X1, Xe, Xs, %4) = ares’ + xgxa! + are’ + xox! 


= m1%3 + Xety’ + wees’ + vane! 
£5(%1, X2, X3, Xa) = Kita AG X49" ate XoX5" | 301" 
= XyK4' 1 XaNx3'  XgX1' + wane’. 
S6(X1, Xo, X3, Xs) = Viva aks A Sete. + Kary 


= eyXa" = Lox" ce 3X9" = 4%". 


Section 9 
_ Theorem 5: Consider the Boolean expression 
f(%1, eerie) Xn) 
= 4jte > Ly 0 ee * 9m Lt km ee * ~~ On 


with 1<mSn. Let ai, as, - - 


-, au be the M=(m-—1)! 
permutations of 1, 2,- +--+, m which are m-cycles. Let 


k;=a;,(2) and let 


| h;(x1, ty Xn) 
urs ay) (ta Xe.) = ~~ (Xm, 4 Hey) PmaXere «> + Xn. 
Then hy, to, - - - , Aw are all the members of Zps(f) and 
mLZAps(f)=n--m—1. 
Proof: Let 
H(i, Xo, °° * > Xm) = tik 2 + Am + 0y'%2' + + Bm. 
By Theorem 4, we know that for 
U(r), 85, 2 +» Xm) == (Hike >> hm Kal te! Kye’) 
the members of Zsp(F’) are 
3S ag B= 12> 3-, Me = mm — 1)! 


t=1 


Dualizing this we conclude that the members of Zps(F) 


ALC 


Now there are only two points Xx -- 
| X1!%—! °° 


™ 
jC ARE ea ie 
i=1 
ity, and 
Lin! Xm4iim42 °* * X, in the point-complex of 


f (x1, X2, + + * , %n) in the m-cube S,, and these two points 


“are contained in the m-cube S,, given by 


The point-complex of F(%1, %2, ° 


H (21, 2, °° * 5 Xn) = Xmpimp2 °° * Xn. 


- +, Xm) in this m-cube 
S,, is the same as the point-complex of f(*1, 2, + + + , Xn) 
in S,. Therefore, by Theorem 1, the members of 
Zps(f(x1, 2, ° + > » Xn)) are obtained by multiplying each 
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of the above M members of Zps(F(x1, x2, - - + , X%m)) by 
F(x, %2,+ ++, Xn); that the length of each of these 
products is »+m-—1 is clear; and this proves the 
theorem. 


Example: Take m=4 and n=6 so that 
f(a, 2, X3, X4, Xs, Xe) = XiXwX3XaN 5X6 oe 41! x9! 09! x4 L5G. 


Then using the previous example, we have the following 
(4—1)!=6 minimal ps expressions equivalent to f, and 
each is of length 6+4—-1=9. 


hy(x1, Xo, X3, X4, V5, x6) 


= (a1' + a%2)(aH2" + a3) (a3" + a4) (aa! + x1) x56. 


ho(x1, %2, %3, X4, Xs, Xe) 

= (a1' + %2)(a2" + 24) (x3’ + 21) (a4! + 23) U5. 
hs(x1, 2, %3, X4, Xs, Xe) 

= (a1! + 3) (a0! + 24) (a3" + x2) (aa! + x1) x5xX6. 


ha(a1, %2, 3, 4, Xs, Xe) 

= (a1 + 23) (a2’ + 2x1) (43" + 24) (a4’ + x2) 5X5. 
hs(x1, %2, 3, %4, Xs, X¢) 

= (x1' + x4) (a2! + 3) (a3" + 21) (axa! + 2x2) x56. 
he (x1, %2, X3, X4, Xs, X¢) 

= (ax1' + x4) ("2" + 41) (a3" + a2) (v4’ + 23) X56. 


Theorem 6: Let the notation be as in Theorem 5. Then 
LZsps(f) =LZps(f) =n+m—1 so that Zps(f) CZsps(f). 
If m<n then Zps(f) =Zsps(f) so that 1, kn, + - - , hy are 
all the members of Zsps(f), and if m=n then f, Mn, 
ho, + + + , hy are all the members of Zsps(f). 

Proof: Since D(f) consists of two points (say 
Py =%X1%_ + + > x, and P2=x1'%o! + - > Xm! Xm41 > * + Xn), the 
considerations of Section 5 tell us that as possible candi- 
dates for Zsps(f) we haee to consider the following two 
cases 

(i) h~f with h of type product-sum; 
Gi) h=q+@~f with gq: and g of type product-sum, 
D(q) =P; and D(q) =P. 

By Theorem 5, the minimal length for h of type (i) 

is ntm-—1. In case (ii) Theorem 3 tells us that 


L(h) = L(qi) + Lge) + 1 
=(m—-1)+(m—-1)+1=x+n-1. 


Therefore, 
LZsps(f) = min (n+ m—1,n+n— 1) 
=n+m—1 
= LZps(f). 


Hence Zps(f) <Zsps(f), 7.e., In, he, + ++, hw are all the 
distinct members of Zsps(f) which arise from case (i). 
Now case (ii) can arise if, and only if, n+m-—1 
=n+n—1, 1.e., =m, and in that case Theorem 3 tells 
US Gi =X1%_ + + + Xn and go=X1/Xq" « - os 
so that h=f. 


uz Mins Mudeien re a apc 
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Example: First take m=4 and n=6 so that 


if / / / 
f(%1, %2, %3, %4, Xs, Xe) = LiXoVgxaXex— + X1' xe! Hs' x4’ Coxe. 


Then the (4—1)!=6 expressions My, he, hz, ha, hs, he given 
in the previous example are exactly the minimal sum of 
product of sums expressions equivalent to f and their 
length, 2.e., LZsps(f) =g. 

Second, take m=n=4 so that 


f (41, %2, %3, %4) = H1%2%3%4 + 44/%9! x3/ x4". 


Then there are seven minimal sum of products of sums 
expressions equivalent to f, and the length of each of 
them is 4+4—1=7, namely: 


f(%1, %2, %3, %4) = %yXex3%4 + Hy He! x3' 44". 

hy(a1, %2, %3, %4) =(x1' + 2X2) (x2! + ws) (x3' + 24) (x4' + 21). 
hio(%1, %2, %3, %4) = (ay! + x2) (x0" + 34) (x3/ + 21) (aa! + 23). 
Mi3(%1, %2, %3, 4) = (41 + 43)(ae2" + 24) (a3! + xe) (a4! + 21) 
Iog(%1, %2, 3, 04) = (41' + x3) (a2! + 21) (a3’ + 44) (x4' + x2). 
his(%1, X2, %3, 4) = (41' + 24) (a2! + 23) (%3’ + 41) (x4’ + x2). 
he (x1, %2, %3, X4) = (ar! + x4) (a0’ + x1) (a3! + x2) (a4’ + 2x3). 


Section 10 


Let f(«1, - + + , Xn) be a Boolean expression such that 
D(f) consists of two points P; and P:. We want to find 
Zsps(f). If P: and P2 form a 1-cell then Zps(f) is given in 
Theorem 2, and as will be shown in our forthcoming 
paper,® Zsps(f) =Zps(f). Now assume P; and P2 do not 
lie on a 1-cell, z.e., that P,; and P2 are the only cells in 
C(f). Let S;, be the m-cell spanned by P; and P2; then 
Ge Webel, Ce hp) <5. Kg) = Sm) 5 banda, gaa) 
=I(Sn) and (ji, je, > * + , jt) =J(Sm) so that stt=n—m. 
Let 
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Va = Xkq for’ gly Zc, 00 
Vatq = hig LORG = 12a sane 
Yn+s+q = Lay for i Ly 2, Pale be 
Then 
F(x, © * ©, Bh) Rp ki, OF ee, ee 
7 ky Kha” eae Xm Viyig ots MisNjy Xin pa aa 
= yiya ss Int yy! + + Ym Imp rYmt2 * * + Ins 


and we are thus reduced to Theorem 6. 

Example: First, take n=3 and suppose that f(x, x2, x3) 
= 04X9'%3'+2x1X2’x3. Then both the “points” x x2’x3’ and 
XX2'x3 lie on the 1-cell xix2’, and hence f itself is the only 
minimal sps expression equivalent to itself. 

Second, take n=6 and suppose that 


f(%1, Xo, 3, La, Xp, Xe) = LoX—UsXexsK1 + Xo'xgxQ X4' X15 
Set 
Vi = X2, Vo = Xs, V3 = Xe, V4 = Xs, V5 = X38, Ve = X11. 
Then 
Hf (a1, %2, %3, Xa, Xs, Xe) = f* (1, Yo, V3, V4) Vo Vo) 
= yiyoysyaysye + yi'yo'y3' Va Yoyo. 


Hence, we can apply the first part of the previous exam- 
ple with y, in place of x;. In other words, we can apply 
that example with the following relabelling of the x,: _ 
X1—X6, X27 X1, X37 X5, Ns — Xs, X5—7X2, Xpe—X3. 
Third, take »=4 and suppose that 


f(%1, %2, %3, %4) = LeeyKaxe + x5! x4) x9! a1’. 
Then 
{G5 %2, X3, %4) = HyXoxgH4q + x1! H 9! x5 24" 


and this is just the second part of the previous example. 
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Summary—When a character printed in magnetic ink passes be- 
neath a magnetic reading head it generates a waveform. A method 
is described here to determine this waveform from the geometry of 
the printed character. The character shape is divided into elementary 
vertical units and the height of these units is then tabulated. A 
single set of experimental data is obtained in the laboratory by 
passing a magnetically printed bar beneath the same reading head 


| which will be used to read the characters. Formulas are derived in 
| this paper for combining the geometrical data obtained from the 


character with the experimental data obtained from the laboratory 
run with the head. A method is then described for programing a 
high-speed digital computer to derive the waveforms. The discussion 
explains why these waveforms are often superior to what might be 
obtained in an actual laboratory run of the printed characters. 


INTRODUCTION 


HE purpose of this article is to describe the theo- 
retical aspects and the computational problems 
involved in synthesizing waveforms generated 


from magnetic ink printing which passes beneath a 
_ magnetic reading head. 


This problem arises from investigations in character 


recognition. Printing characters in magnetic ink for 
automatic reading and recognition during document 
: processing is attractive from several points of view. 
Magnetic printing is highly resistant to most kinds of 
mutilation. The application of fountain pens, ball point 
pens, pencils, stamp pad, ink, dirt, or grease to the sur- 


face of magnetic ink printing has practically no effect 
upon the waveforms it generates. Also, folding, bending, 


_crumpling or otherwise abusing the document has little 


effect upon the waveform generated by the magnetically 
printed character. Furthermore, the methods for print- 
ing such documents can be those commonly employed in 
the printing field. 

Experimental waveforms derived from magnetic 
printing are not determined solely by the originally 
designed geometry of the characters. These waveforms 
are affected by the chemical composition of the ink, the 
squeeze-out around the impression, the fibre structure 
of the paper, the absorbency of the paper and other 
similar factors. These factors produce a three-dimen- 
sional configuration substantially similar to the original 
design but with contaminations which produce aberra- 
tions in the waveform. To facilitate design and to study 
individually the effect of the variables of printing, it is 
desirable to know the waveform generated by an “ideal” 


* Manuscript received by the PGEC, February 27, 1958; re- 
vised manuscript received, June 19, 1958. 
+ Remington Rand Univac, South Norwalk, Conn. 
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_ A Method for Synthesizing the Waveform Generated 
_ by a Character, Printed in Magnetic Ink, in Passing 
Beneath a Magnetic Reading Head’ 
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printed character free from the contaminations of print- 
ing, paper and inks. This article describes a procedure 
for deriving these waveforms by using digital computer 
techniques, and also presents the rationale for this 
procedure. 

Such a synthesizing procedure also enables the de- 
signer to alter the characters to qualitatively change the 
waveforms generated by them in a predictable manner. 


ANALYSIS OF THE PROBLEM 


It is required to find the voltage induced in the coil of 
a reading head while a character printed in magnetic ink 
and previously magnetized passes beneath the head. 
This voltage is directly proportional to the rate of 
change of the magnetic flux through the coil of the read- 
ing head. The total flux of the core material passes 
through the coil. The flux in the core is determined by a 
number of factors, the most important of which are the 
geometry of the core and the characteristics of the 
material from which it is made, especially its permeance. 
Thus, as a fixed magnetic pattern is placed beneath the 
core, the flux passing through the core will be deter- 
mined by the factors just mentioned—its geometry and 
physical characteristics. 

The pattern surrounding a magnetic geometric form 
is not easily represented in an analytic equation. In the 
very simple case of a magnetized rectangular area, for 
example, the flux at any point in the immediate vicinity 
surrounding this rectangle (even assuming point mag- 
netic sources) is represented in terms of a vector sum of 
two quadratic terms. 

To start with a given geometric form and to determine 
analytically the voltage induced in the coil appears to 
be a rather difficult problem. The first step, to find the 
flux concentration around the geometric form, is a task 
in itself. In addition, as mentioned previously, the flux 
picked up by the core depends on the geometry and 
physical characteristics of the core. 

The approach here is, therefore, to start with empiri- 
cal data. Such data would be in the form of waveshapes 
experimentally recorded in the laboratory when some 
basic unit of area passes beneath the reading head. It is 
next necessary to decide what basic unit will provide the 
best basis for a computational approach. 

The waveshape produced by a rectangle or bar pass- 
ing beneath a reading head is well known. It consists of 
a positive pulse generated by the leading edge of the bar 
and a negative pulse generated by the trailing edge of 
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BAR 


Fig. 1—(a) Printed bar. (b), (c) Two printed steps, AA(0) and 
—AA(t,), which together make up equivalent response as bar 
(a). 


the bar. (The direction is of course arbitrary.) Such a 
bar might comprise a fundamental unit. It was felt, 
however, that a step would provide a better approach. 
The bar can be simulated by the sum of a positive and 
a negative step, as illustrated in Fig. 1. 

It can be easily seen that as the step approaches the 
reading head, more and more flux lines will pass through 
the core. The rate of change of these flux lines will be- 
come maximum when the leading edge of the bar enters 
beneath the gap. As the full area of the bar passes under 
the gap, the flux concentration becomes maximum. The 
rate of change remains zero while the full area of the bar 
is under the gap. What is generated by a step is therefore 
a positive pulse which is centered at the leading edge of 
the step. The unit differential of the area will be the 
basic unit for finding a waveform of the character. The 
step and the voltage it generates are shown in Fig. 2. 
Such a step generates a unit pulse which is indicated by 


AA(0) — e(é). (1) 


A multiple of the unit step will generate a pulse which 
is a multiple of the pulse generated by a unit step. This 
might be indicated as 


kAA(0) — ke(t) (2) 


where AA(0) indicates the unit differential of area and 
(0) indicates that the leading edge of the step passes 
beneath the gap at time zero. The arrow is used here to 
indicate that the area on the left-hand side generates 
the voltage on the right-hand side. This is not an equa- 
tion, but rather a device used by mathematicians to in- 
dicate the mapping correspondence between two classes 
of functions. 

The time relationship, which states that a unit step 
passed beneath the head at a later time than that of 


IRE TRANSACTIONS ON ELECTRONIC COMPUTERS 


December 


ae...) 


(a) 


e@(t) 
(b) 


Fig. 2—(a) Unit step, AA (0). (b) Voltage, e(¢), generated 
by unit step. 


AA(0) will generate a pulse delayed by the same amount 
is indicated by 


AA(r) > e(t — 7). (3) 
Eqs. (2) and (3) are combined as 
kAA(r) — ke(t — 7). (4) 


The final necessary relationship states that two steps of 
different size and occurring at different times produce 
voltages which are linear and additive; it is expressed by 


k, AA (és) + Ro AA (te) z=? kyie(t == t1) -- Roe(t = to). (5) 


This last expression indicates that the three factors— 
proportionality, a linear time relationship and an addi- 
tive voltage relationship—are all necessary for our com- 
putational process. These factors can be easily derived 


from the principle of superposition and the two equa- 
tions 


¢= F/R (6) 
and 
do 
jy SES 7 
dh (7) 
PROCEDURE 


The character in which we are interested consists of 
one or more bounding curves and is a connected closed 
‘area. The important aspect of this area is its direction 
of travel with respect to the reading head. This will be 
called the ¢ axis. Such a form is shown in Fig. 3(a). Its 
area, integrated normal to the ¢ axis, is shown in Fig. 
3(b). The units to be used here must correspond to the 
definition of the unit step. This area is now divided 
along the ¢ axis into slices of length At. At will correspond 
to that used later in reference to the voltage curve e(t). 
The number of differentials of area in the ith slice A; is 
now found for integral multiples of Af. A plot of such 
values of the A,’s is shown in Fig. 3(b). 

The difference in area between one slice and the rest 
will generate a voltage pulse. The amplitude of this 
pulse is proportional to the difference in area between 
the two slices. The difference between the ith slice and 
the (¢—1)th slice may be called 


6; = Ai — Agi. (8) 
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(a) 


(b) 
Fig. 3—(a) Printed pattern. (b) Plot of area of (a). 


Because of the properties of linearity, proportionality 
and additivity, these relationships may be super- 
imposed. 

_ The data taken in the laboratory for the voltage de- 
rived from the printed step are now “calibrated,” using 
_the same value for At as above. The amplitude is plotted 
against time as in Fig. 2 and the maximum is reached at 
time ¢=0. This “calibration” yields a series of significant 
values of e(¢) with 


e = (Al) (9) 


where & can vary from —L to + M; e; is approximately 
zero for values of k less than —L and greater than + M. 
The first area slice or step area generates a voltage 
which may be called E,(¢) and which will be equal to 
5;e(t —At). It is illustrated in Fig. 4(c) for the area shown 
in Fig. 4(a) whose differentials are in Fig. 4(b). At the 
_time jAt the voltage £,(jAt) will be called £;;. Then 


Eur = 61° €0 ; (10) 
and 
Exo = 61°€1 Eyo = 61-1 (11) 
Or 
Hi 6-74 (12) 
where 
eee FM (13) 
and 
Ey, =0forj<—-L+1,j>M+1. (14) 


F(t) and E;(t) are similarly defined and illustrated in 
Fig. 4(a) and 4(e). The general term E;,(#) is defined in 
the same manner: 


E;,(t) = dxe(t — RAL). (15) 
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E33 
Fig. 4—(a) Area for waveform. (b) Differentials of area. (c) Re- 


sultant voltage for 6). (d) Resultant voltage for 62. (e) Resultant 
voltage for 63. 


Define 
Ei, (Ai) = Eyy (16) 
so that 
Exx = 64° €6 (17) 
and 
Expt = 04°61, Cg ees eS (18) 
or 
Exuj = Oxej—x (19) 
where 
stele ee (20) 
and 
Ca OO a ee ee (21) 


The output voltage Er from the coil, which is to be 
synthesized, is found by adding together all the voltages 
derived at one given time from all the steps which 
create an output at that time. The total voltage pro- 


duced at time 7 is given by 
peel att ont arta (22) 


so that 
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(a) 


(b) 


t=0 
Fig. 5—(a) Ei, E2, and E; of Fig. 4 superimposed. (b) Sum. Er. 


Erm = En t+ Ea +s: + + Engin (23) 


and 


N 
Er: = >) Ex where Ex1 = Ofork>L+1 (24) 


k=1 


and in general 


N 
Eq; = S. Exj-L4i<j<N+M 
k=1 


(25) 


where 


dij —= 0 for k<j—M and k>7 + L.- (26) 


The curves for Fi, E, and E; for the area of Fig. 4(a) 
are shown superimposed on Fig. 5(a), and the result of 
adding these curves together is shown in Fig. 5(b). 


COMPUTER CALCULATION 


The data to be supplied to the computer will be the 
G+M-+1 constants ep, e¢_1u,-- €0, @1,° °° 
€m—1, €m. Lhese are the constants obtained from the ex- 
perimental data. Also supplied to the computer as raw 
data are the 6,’s, which are obtained by the following 
method. An enlarged image of the character whose 
waveform is to be synthesized is drawn on quadrangle 
paper. The grid to be used would consist of lengths At 
in the horizontal direction and AA in the vertical direc- 
tion. From this is prepared a chart giving the boxes for 
each slice ending at kAt, which is found by simply 
counting the squares in the vertical column bounded by 
(k—1)At and RAt. 

The computer then makes these calculations: 


7 RANg 
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OE I A; = Ap-1 (27) 
Ex; — On€j—K (28) 
N 
Er; = >) Bij (29) 
c— 
where 
Ey; =O0fork>j+L and k<j—M (30) 
and 
-L+1<j<N+4+M. (31) 


The calculations required make use of subtraction for 
(27), multiplication for (28) and a series of additions for 
(29). These calculations can be simply programed on 
any modern computer. 

The result of these computations can then be graphed 
on the same time scale in units of At. Notice that this 
output voltage has values which range from —L-+1 to 
N-+M. This is because the voltage in the coil begins be- 
fore the printed form reaches the head gap and remains 
after the form has passed beneath the reading head gap. 


ILLUSTRATIVE EXAMPLE 


As an example of how this method is used, the wave- 
form for the printed character “0” has been synthesized. 
The character itself is shown in Fig. 6(a). In Fig. 6(b) 
the area is plotted along the ¢ axis; each of the vertical 
values is one of the A;’s. The increments of area desig- 
nated as 6, appear in Fig. 6(c). Although these 6,’s are 
plotted, they are not necessary for the computations 
and are not supplied in the form of raw data. The result 
of the maximum computations appears in Fig. 6(d). 
Notice that this output voltage consists of a positive and 
negative pulse followed a little later by another positive 
and negative pulse. These pulses correspond to the large 
increments in area. The first pulse is caused by the left- 
hand edge of the zero; the first negative pulse is caused 
by the right-hand trailing edge of the left hand bar of the 
zero; second positive pulse is caused by the left-hand 
edge of the right-hand bar of the zero; the second nega- 


tive pulse is caused by the right-hand edge of the char- 


acter. It is generally true that the waveshape generated 
by any character printed in magnetic ink will be deter- 
mined by the larger increments in an area. 

Although the two characters “0” and “8” appear to 
be quite different as shown in Figs. 7(a) and 7(b), the 
waveforms generated by these two characters are al- 
most identical. An automatic character-reading system 
could not use the characters whose generated waveforms 
correspond to such a degree. This difficulty could be 
overcome by using a different figure “8,” as shown in 
Fig. 7(c). This character would generate four sets of 
positive and negative pulses instead of the two sets of 
positive and negative pulses which would be generated 
by the character in Fig. 7(b). 

Through the study of printed characters and the 
waveforms which they generate, character designs can 
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(a) 


(b) 


(c) I "" i " 


(d) 


Fig. 6—(a) Sawaal “0.” (b) Vertical components, Ax. (c) 
k for “0.” (d) Er() for “0.” 


be made which are as easily distinguishable from the 
waveforms that they generate in magnetic reading as 
they are distinguishable in human reading. This is one 
of the needs that the synthesis approach fulfills. 


DISCUSSION OF RESULTS 


The results obtained from such an automatically pro- 
gramed machine computation agree very closely with 
the waveforms taken in practice in the laboratory from 
printed samples. The question arises, why not simply 
take the waveforms in the laboratory? The reason for 
not using the laboratory samples is that printing intro- 
duces certain distortions which cannot be closely con- 
trolled. To average the results taken from a given batch 
of printing would not be correct, since some of the dis- 
tortions are due to long term parameters. Such parame- 
ters would be present in a given batch and would change 
only from one batch to another. An example of such par- 
ameters would be the type of impression used by the 
printer or the amount of ink in the fountain of the print- 
ing press. 

The theoretical results use printed samples to derive 
only the function e(¢). A square of magnetized metal 


Flores and Ragonese: Waveform of a Magnetic Ink Character 


281 


(a) 


(b) 


Fig. 7—(a) Character “0.” (b) Character “8,” same waveform 
as (a). (c) Character “8,” waveform different from (a). 


may give a more exact measure of e(¢) . This function 
may be studied and also smoothed to reduce the eccen- 
tricities of the printing process, if printing is used to 
derive it. The waveforms predicted by using such a 
smoothed function, e(¢) , will be more suitable theoreti- 
cally than those derived by averaging the laboratory 
waveforms. 

Somewhat better results are obtained when the edge 
phenomenon is taken into account. Although the edge 
of the bar under consideration may be designed to print 
as a straight line, it looks more like a topographical map 
of the Rocky Mountains when viewed under a micro- 
scope. The ink is transferred to the paper from the 
printing roll. The way the ink prints on the paper de- 
pends on the roller pressure, the amount of ink on the 
roller, the viscosity of the ink, the absorbency of the 
paper, the fibre structure of the paper, etc. The ink dis- 
tribution at the edge which is designed to print as a 
vertical line may be approximated by a probability 
density distribution with very good results. The dis- 
tances being considered are in the neighborhood of frac- 
tions of a thousandth of an inch. This is well below the 
resolution of the head. The probability density distribu- 
tion which was arrived at from physical considerations 
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(b) 


(c) 


Fig. 8—(a) Ink density at edge. (b) Vertically integrated ink dens- 
ity. (c) Probability density approximation to vertically inte- 
grated ink density. 
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is very close to the actual flux distribution prevailing in 
the head. Fig. 8(a) illustrates the edge effect described, 
(b) the vertically integrated area at the edge, and (c) 
the probability density distribution. 


CONCLUSION 


A method has been discussed for deriving the wave- 
forms generated as a character printed in magnetic ink 
is passed underneath a magnetic reading head. The 
method employs the geometry of the character and the 
empirical properties of the magnetic reading head, the 
ink, the printing process, etc., involved. 
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On the Minimum Logical Complexity Required for a 
General Purpose Computer’ 
S., P. FRANKELT 


Summary—A definition is provided for the term ‘‘general purpose 
computer” (gpc) which is compatible with usage and is analogous 
to, but distinct from, Turing’s definition of a ‘‘universal computer.” 
A gpc is presented in functional and logical design which seems to 
approximate the minimum complexity consistent with this definition. 
No specific definition of complexity nor a bound on required com- 
plexity is presented. 


DEFINITIONS 


HE term general purpose computer (gpc) has been 
applied to a variety of digital computing instru- 
ments having the following characteristics: 1) The 
ability to perform various arithmetic operations on 
numbers in digital form. These operations usually in- 
clude addition, subtraction, multiplication, and divi- 
sion. Numbers are usually held in decimal or binary ex- 
pansion to a convenient precision (number of decimal or 


* Manuscript received by the PGEC, May 22, 1958. 
+ Consultant, Logical Design. 1764 Redondo Ave., Long Beach, 
Calif. 


binary digits) and are usually signed. That is, the arith- 
metic operations are performed on algebraic numbers 


- expressed to limited precision. 2) The instrument in- 


cludes a memory store, which holds both the numbers 
taking part in a calculation and the instructions, which 
control the execution of elementary operations. The se- 
quence in which instructions held in the memory take 
control of activity is dependent (in some situations, at 
least) on the numbers produced by arithmetic opera- 
tions. One familiar form of this branching is the selection 
of one or another subsequent chain of instructions de- 
pending on the sign of a number appearing in the calcu- 
lation. 3) Numbers produced in the course of the 
calculation may be stored in the memory and thus made 
available for use at later stages of the calculation. 

The appellation “general purpose” implies that by 
suitable choice of the numbers and instructions initially 
set in the memory, the computer may be enabled to 
carry out any one of a broad class of calculations. A 
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uantitative limitation on the breadth of this class is 
posed by the capacity of the memory, i.e., by the 
mount of instructional and numerical information it 

able to hold. To remove this uninteresting limitation 
e term general purpose will be used here in a strict 
nse requiring that the memory capacity be arbitrarily 
tensible without change in the logical design. By add- 
g the clearly reasonable restriction that the informa- 
on initially set in the memory (the program) be of 
inite extent, the present definition is made similar to 
the Turing definition of a computer. If, moreover, the set 
{ elementary operations is of ordinary flexibility, the 
eneral purpose computer has the essential property of 
Vuring’s universal computer,' the ability to perform any 
yomputation which could be done by any computer. It 
5 the object of the present paper to display a gpc logical 
lesign of nearly minimal complexity compatible with 
this practical approximation to the requirements for a 
iniversal computer. 

At each step of its operation the computer proper oc- 
tupies one of a finite number of active states. The term 
computer proper describes those parts of the computer 
tubject to immediate change; as opposed to the memory, 
most parts of which are quiescent in each step of the 
computation. One measure of complexity considered 
nere is the logarithm (to base two) of this number of 
active states, z.e., the number of binary digits required 
(0 specify the coneligion of the computer proper. 
| The computer described below is of the synchronous 
‘ype. It makes use of a clock which presents timing 
pulses by means of which all activities are synchronized. 
The time interval between successive clock pulses is 
-alled a digit period. In each digit period the computer 
oroper receives from the memory a body of information 
which is typically very small compared with the full 
memory capacity. This information is received as a dis- 
“rete set of binary digits (bits). Similarly, a few bits of 
nformation are deposited (recorded) in the memory in 
each digit period. These recorded bits are functions, as 
Jescribed in the logical design, of the current active 
state and of the bits received from memory. At the time 
sf each clock pulse the active state may change. The 
1ew active state is also a function, described in the logi- 
sal design, of the prior active state and of the bits re- 
seived from the memory during the past digit period. 

Various constituents of the complexity of a logical de- 
sign are: 1) the number of bits specifying the active state; 
2) the number (or mean number) of bits received from 
memory in each digit period; 3) the number of bits re- 
sorded in the memory in each digit period; 4) the com- 
plexity of the functions determining the bits recorded 
and the new state. These constituents may be combined 
(and the fourth given a quantitative measure) to pro- 
vide an over-all measure of complexity. Since no formal 
minimization of complexity is claimed for the design 


1A. M. Turing, “On computable numbers,” Proc. London Math. 


Soc., series 2, vol. 42 (1936-1937), pp. 230-265. 
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Computer of Minimum Complexity 


presented below, a choice of such an over-all measure 
need not be made here. 


FUNCTIONAL DESCRIPTION OF M’AC 


The computer described here is designated M’AC 
(from MicrocephalAC). Its memory organ is a magnetic 
tape which presents and receives information in several 
channels. Each of two channels is served by two heads 
spaced by integral multiples of the distance correspond- 
ing to a digit period. One of these channels acts as a cir- 
culating register: the first of its heads records a bit in 
each digit period, the second (in the direction of tape 
motion) presents to the computer proper the bit which 
was recorded in a correspondingly earlier digit period. 
The other of these channels holds the main memory. 
The distance separating its two heads is an integral 
multiple of the distance between the circulating register 
heads. In some digit periods both heads read from the 
main memory; in others one (the memory writing head) 
records while the other reads. To avoid the necessity of 
reversing the direction of tape motion, with consequent 
disturbance of the storage of information in the circulat- 
ing register, the tape is given the form of a loop. 

Since the tape loop is of finite size the memory capac- 
ity is finite; thus, M’AC is not a computer in the sense 
of Turing. However, the loop size need not be fixed by 
the logical design but may be treated as a parameter 
available for choice in the planning of each computation. 
(In the same way the memory capacity of a typical gpc 
might be left unspecified in its logical design, thus justi- 
fying the use of the term “general purpose” in the 
stricter sense described above.) 

In addition to the two channels described above, the 
tape may hold permanently recorded signals used to 
generate the clock pulse and two bits of timing infor- 
mation. Successive groups of m digit periods are called 
“word periods.” The last digit period of each word 
period is distinguished by a signal, denoted Z. The ab- 
sence of the signal L is denoted Z and marks bit periods 
other than the last of a word period. Another signal, NV, 
and its complement, WN, distinguish alternate word 
periods. Word signals in which JV is present are called 
number word periods, those in which N occurs are instruc- 
tion word periods. The mechanism by which the signals 
N and W are generated may consist of a flip flop oper- 
ated as a scale-of-two counter, changing its state at the 
end of each digit period in which Z occurs. For reasons 
described below the signal N, unlike L, is not presented 
directly by the tape. 

The distance through which the tape moves in a word 
period will be called a word length. The separation be- 
tween the two heads of the circulating register is two 
word lengths, so that if the bit put on the tape by the 
recording head (denoted R’’) were always that being 
presented by the reading head (denoted R), then an arbi- 
trary sequence of 2” bits would be recirculated indefi- 
nitely. In fact, the equality R’’=R holds often but not 
invariably, in keeping with customary usage of the term 
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circulating register. Similarly, the separation of the two 
main memory heads is an integral multiple of two word 
lengths so that a word recorded by the recording head 
in a number word period will be read by the reading 
head in another number word period, several word 
periods later. The number of bits in a word, 7, is not 
fixed by the logical design but, like the length of the tape 
loop, may be treated as a parameter available to the 
user’s choice. 

The 2n bits stored in the circulating register consti- 
tute two words: that presented (or recorded) during the 
N word period is a binary number resulting from previ- 
ous computation and held for use in the next arithmetic 
operation; the bits presented or recorded in the NV word 
period are related to the instruction awaiting execution. 
Thus, the circulating register plays the roles of both the 
accumulator and instruction circulating register of a 
simple binary digital computer. Similarly, information 
read from or recorded into the main memory is divided 
by the signals NV and W into two categories: bits read 
during an N word period will be regarded as forming an 
instruction; those read or recorded during an N word 
period (recording into the main memory occurs only 
then) will be called a number. This distinction would 
appear to separate the main memory into disjoint halves; 
an instruction memory which can only be read and a 
numerical memory which can be recorded into or read. 
Such a separation would uncomfortably restrict the 
programmer’s abilities. Accordingly, it is avoided by 
imposing the requirement that the length of the tape 
loop be an odd number of word lengths. Then each word 
of the main memory appears in successive tape revolu- 
tions in the guise of a number and an instruction 
alternately. 

Two phases of operation are distinguished by a flip 
flop, denoted £. Its on state marks the execute phase in 
which an instruction is executed and a new instruction 
then read. An execute phase typically consists of a 
number word period followed by an instruction word 
period. Following this the delay phase marked E may be 
entered to await the appropriate time for the execution 


of the newly read instruction. Since the memory of - 


M’AC is entirely serial, an address consists only of the 
specification of the time of execution. It is presented by 
the main memory as a delay number, an integer de- 
scribing the duration of the delay phase to follow in 
units of two word periods. (If the delay number is zero 
no delay phase is entered and the execute phase is con- 
tinued for another pair of word periods.) 

During the instruction word period of the execute 
phase, 7.¢., in the circumstance VE, a new instruction is 
set into the circulating register. It is derived from the 
bits presented by one of the two main memory heads, 
denoted M and M* respectively. In subsequent in- 
struction word periods of the delay phase the delay 
number is reduced by one with each circulation until it 
becomes negative, 7.e., until a one appears in its most 
significant bit position. At that time the execute phase 
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is entered. (If the delay number as read from main 
memory was negative this condition arises immediately, 
hence no delay phase intervenes.) 5 

The elementary operations of M’AC are intended to 
provide a minimum set of operations into which the 
activities of a computation can be broken. These clearly 
do not include the extraction of a square root, which in 
most computers must be programmed, 1.e., built up 
from more elementary operations under the control of a 
sequence of instructions. Some quite powerful comput- 
ing instruments have similarly dispensed with a wired- 
in division procedure. Carrying this spirit of parsimony 
a step further, the programmer may also be required to 
build up multiplication from the simpler arithmetic 
operations, which would thus seem to be only addition 
and subtraction. (Operations of a less sharply arith- 
metic character, like the familiar extract, can be ac- 


complished if means for doubling and branching on the — 


most significant bit are provided.) However, the list can 
be reduced still further by the omission of addition. By 
twice subtracting a number from itself, its negative 
may be formed. (In M’AC, as in many binary com- 
puters, the most significant bit of a word is regarded as 
having a negative position value.?) By use of this com- 
plemented number and the operation of subtraction, 


the effect of adding the original number may be ac- — 


complished. 

The basic set of operations required for a gpc thus 
appears to be 1) subtract, 2) record in memory, and 3) 
branch. It is not necessary, however, that these opera- 
tions have separate orders. In M’AC each instruction 
execution accomplishes all three of these operations— 
a number read from memory (by head M) is subtracted 
from that held in the circulating register, the result is 
simultaneously recorded in memory (by head M*), 
and the next instruction is then read either from M or 
from M*, depending on the outcome of the subtraction. 
(Specifically, the criterion chosen is the state in which 
the borrow flip flop, B, is left following the number 
word of the execute phase.) Typically, one of these 
three operations will be the desired activity; the other 


» 


two will be unneeded by-products. Thus, the simplifica- — 


tion of reducing the order list to one item is obtained at 
the cost of less efficient use of memory capacity and 
greater complexity of the programmer’s task. 

Since only one order is provided, no bits of the in- 
struction need be used for order specification. Accord- 
ingly, the entire instruction word is used for the delay 
number. The arithmetic unit, consisting of the borrow 
flip flop, B, and associated circuits, is used to effect the 
reduction of the delay number in instruction word 
periods of the delay phase, marked WE. This is accom- 
plished by presetting B to the on state just prior to 
these word periods, at the times NEL, and proceeding 
thereafter as if the number zero were being subtracted. 


2S. P. Frankel, “The logical design,” IRE TRANS. ON ELECTRONIC 
Computers, vol. EC-6, pp. 5-14; March, 1957. 
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Three of the constituents of the complexity of the 

‘AC design may now be summarized: except during 
he word periods marked NE, three bits of information, 

, M, and M*, are received from memory and one bit, 
R’’, returned. In word periods NE two bits, R and M, 
re received and two returned, R’’ and M’’. The active 
state is determined by the two bits, E and B, together 
vith the two bits of timing information, L and JN. 
| 


LoGIcAL DESIGN 


In every word period the arithmetic unit is engaged 
nm performing a subtraction, the result (remainder) of 
which is recorded as R’’ except in word periods WE. 
[he remainder is also recorded as M’’ in word periods 
VE. The minuend consists of the bits R at all times. 
The subtrahend bits are zero during the delay phase, 
M in word periods VE, and R in word periods VE (to 
snsure that no change in the state of B occurs then). 
Thus the subtrahend, denoted s, is expressed by 
#(MN+RN). The borrow flip flop, B, is ordinarily 
turned on by a one-bit in the subtrahend unmatched in 
the minuend, off by a one in the minuend and zero sub- 
crahend. Otherwise it remains unchanged, except that 
't is turned on at the end of a number delay word period 
(to provide for the diminution of the delay number) 
and off at the end of any instruction word period. Alto- 
gether, the condition for setting B on is sR+NEL, for 


285 


turning it off is 3R+NL. The remainder bits, denoted 
r, may be described as R@s@B where © represents the 
exclusive or function. 

The bits, R’’, recorded in the circulating register are 
r except in word periods marked NE, when they are M 
or M* as B is on or off. They are thus described by the 
expression, r(V+F)+NE(BM+ BM*). At the end of 
each instruction word period, hence on NL, the bit R’’ 
is set into flip flop EZ. The condition for recording in 
main memory is denoted w, and has the value VE. The 
bits, M’’, recorded then are the remainder bits, 7. 

This description of the logical design is summarized 
in Table I. The condition for turning on flip flop B is 
denoted B’, that for turning it off is denoted B’, etc 


TABLE I 
M’AC Logic 
ak Logical Equations Definitions 
N’=NL B’=sR+NEL s=E(MN-+RN) 
N’=NL B’=3R+NL r=ROsOB 
R"'=r(N+E)+NE(BM+BM*) 
E'=R"'NL 
E’=R"NL 
w=NE 
M"'=r 


Iterative Combinational Switching Networks— 
General Design Considerations” 


E. J: McCLUSKEY, Jr. 


Summary—An iterative network is a combinational switching 
sircuit which consists of a series of identical ‘‘cells” or sub-networks; 
‘or example, the stages of a parallel binary adder. A formal design 
nethod for iterative networks is presented. This is similar to the 
dow table technique for designing sequential circuits. 


1. INTRODUCTION 
| N iterative network is a combinational switching 
| 


circuit which consists of a series of identical 
“cells” or sub-networks. Some typical iterative 
4etworks are shown in Fig. 1 and Fig. 2, and the general 
terative network structure is shown in Fig. 3. There 
are two situations in which it is appropriate to design a 


* Manuscript received by the PGEC, June 10, 1958; revised 
manuscript received, July 30, 1958. 
+ Bell Telephone ‘Labs. , Whippany, N. J.; formerly with Princeton 
Jniversity, Princeton, N. aie 


circuit as an iterative network. The first is when the 
conditions for an output can be conveniently stated in 
terms of the number of inputs which are active,+~* for 
example, when there is to be an output only when 2 or 
5 of the inputs are active. The second situation arises 
when the conditions for an output can be stated in 
terms of the relative positions of the inputs which are 
active. Examples of the second situation are: 1) an out- 
put is desired if any two adjacent inputs are active and 


1 Jn this case the desired switching function is said to be totally 
symmetric. 

2 W. Keister, A. E. Ritchie and S. H. Washburn, “The Design of 
Switching Circuits,” D. Van Nostrand Co., New York, INE Yes Dp: 
55-64; 1951. 

3S. H. Caldwell, “Switching Circuits and Logical Design,” John 
Wiley and Sons, Inc., New WOrkaiNn Yor Chandi moo. So. 

fSy rl. Washburn, “Relay ‘ trees’ and symmetric circuits,” TRANS. 
AIEE, pt. I. vol. 68, pp. 582-586; 1949. 
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Ly = 2° Ayo t2! Xy +22 Lyo +23 Ly 
Xe = 29 Lap t2! Xa} + 2% Laat2? Lag 
Z=%+%2 


Fig. 1—An iterative network for adding two binary 
numbers, x; and xs. 


Fig. 2—A relay iterative network which is closed when an odd 
number of relays are operated. 


Ly Lay Cry Zin Zan rn 


ne ~ TRI 


Zi Za Zt Zin Zan Ztn 
Za tH 

aju—jT™ INpuT To LTH ceLL 

2ji-J™ output FROM LTH ceLL 

yji-J™ CARRY TO LTH CELL FROM (U-1)ST CELL 


Yju-J™ CARRY FROM LTH ceLL To (U+1)8T CELL 


YGU= Yj ttt 


Fig. 3—General structure of an iterative network. 


2) an output is desired if X1=1 and X2=0 for any cell 
1, and there has been no preceding cell 7 for which X:;=0 
and X2,=1. A basic characteristic of the second situa- 
tion is that some ordering among the inputs is assumed. 

An iterative network usually requires less equipment 
and is easier to wire-up than a standard combinational 
circuit having the same switching function. This econ- 
omy in equipment usually carries with it the penalty of 
a longer operating time.’ Because use is made of the 
symmetry of a switching function, it is easier to design 
an iterative network than to carry out a general com- 
binational circuit design, and it is possible to design 
functions of large numbers of variables. 


5 For electronic networks. 
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There are two general types of iterative networks: 1) 
Those which have outputs from each cell (Fig. 1), 
which will be called cell-output networks; and 2) those 
which have outputs only from the last cell, (Fig. 2) 
which will be called circuit-output networks. In general, 
the same design techniques apply to both types of net- 
work. Specific differences will be discussed elsewhere. 

The idea of an iterative contact network was first 
discussed by Keister, Ritchie, and Washburn® (in this 
book the term reiterative circuits is used). A suggestion 
for a formal design procedure using sequential circuit 
flow tables was published by Huffman,’ and a detailed 
discussion of both approaches is given by Caldwell.® 
These presentations are concerned almost entirely with 
contact networks and lack generality since they place 
restrictions on the nature of the allowed combinations 
of carry signals. In this paper a general method for de- 
signing both relay and electronic iterative networks is 
presented. This method makes use of an analogy be- 
tween iterative networks and clocked sequential switch- 
ing circuits. 


II. COMPARISON OF ITERATIVE AND NONITERATIVE 
COMBINATIONAL CIRCUITS 


Before starting the discussion of iterative network 
design, a comparison of iterative and noniterative net- 
works will be made. This will be done by considering 
the design of a parallel binary adder. Since such an 
adder is a combinational circuit, it can be specified by 
means of a table of combinations? such as Table I(a). 
For adding together two two-bit numbers such a table 
has sixteen rows; in general for adding two n-bit num- 
bers, 2?” rows are required. It is clear that it would be 
quite laborious just to write down the table of com- 
binations for a circuit to add together two numbers of 
useful length. The algebraic expressions corresponding 
to Table I(a) are listed in Table I(b). Fig. 4 shows the 
two-stage diode circuit which corresponds directly to 
these expressions. This circuit contains an “and” gate 
with four inputs. In general, a circuit to add together 
two n-bit numbers would require gates with 2n-inputs. 


- Gates with large numbers of inputs are difficult to 


realize physically, and may require the addition of more 
stages of logic. In an iterative network such as that 
shown in Fig. 1, the size of the gates is independent of 
the size of the numbers being added. There is an ad- 
vantage to the network of Fig. 4 in that the only asso- 
ciated delay is that of a signal passing through two 
gates (an “and” gate and an “or” gate) in series. In the 
network of Fig. 1 there can be a delay due to two gates 
in series per cell so that the high-order output bits may 


6 Ob. cit., pp. 55-64. 

7D. A. Huffman, “The Synthesis of Iterative Switching Circuits,” 
Quart. Prog. Rep., Res. Lab. of Elec., Mass. Inst. of Tech., Cam- 
bridge, pp. 63-67; January, 1955. 

SOD. ctt., chal 1, 

* W. Keister, A. E. Ritchie and S. H. Washburn, of. cit., p. 89. 
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TABLE I 


SPECIFICATIONS FOR A Circuit To App Two 2-Bir BINARY 
NuMBERS Xy =X 11X19 AND Xo= XoXo AND PRODUCE 
AN OUTPUT Z=Z2Z,Z 


(a) 
Xe Xy Le 
Xo Xa Xu = =X10 Z2 Zy Zo 
0 0 0 0 0 0 0 
0 0 0 1 0 0 1 
0 0 1 0 0 1 0 
0 0 i 1 0 1 1 
0 1 0 0 0 0 1 
ee 1 0 1 0 1 0 
0 1 1 0 0 1 1 
0 1 1 1 il 0 0 
1 0 0 0 0 1 0 
1 0 0 tf 0 1 1 
1 0 1 0 1 0 0 
1 0 1 1 1 0 1 
1 1 0 0 0 1 1 
it 1 0 1 1 0 0 
1 iL 1 0 1 0 1 
1 1 1 1 1 1 0 
(b) 


Zo=XX10' +X2'X10 

21 =XyX11'X19' + XuX20/ X11! + Xo1/X20/Xu +X! XyX10! 
+Xo1'XoXu' Xt XaXwXuX10 

Z2=XyXuitXurXwX19 +X 9XuX10 


Lai 

Loy 

Lao 

Lo 

x 

of pst tt 

Tor 

=a TT Cit 
eal i 


as Z10 Za 


Fig. 4—Circuit to produce the sum of two two-bit binary 
numbers. 


lag the inputs by an appreciable time. In summary, a 
two-stage electronic network is generally faster than an 
iterative circuit, but may require an excessive amount 
of equipment. The two-stage circuit is also quite labo- 
rious to design, although this design may be mechanized 
on a computer.!? 

By factoring the algebraic expressions for a two- 
stage logic network, it is possible to obtain other circuits 
which may require fewer diodes and which will have 
more of a delay through the circuit. One such factoriza- 
tion for the binary adder is shown in Fig. 5. In general, 
for any table of combinations there will be one factori- 


10, J. McCluskey, Jr., “Minimization of Boolean functions,” 
Bell ae ga J., vol. 35, pp. 1417-1444; November, 1956. 

H. Urbano and R. K. Mueller, aN topological method for the 
eich iston of the minimal forms of a Boolean function,” IRE 
TRANS. ON ELECTRONIC COMPUTERS, vol. 5, pp. 126-132; September, 
1956. 
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General Design Considerations 


Zo = Lao Lo (Lay + Ly) + Lay Lys 
Zy = (HS, + Ly) (Loi + Li) (Lo + Xho) + Lay Ly (Loy + Lfy) ('ay + X41) 
Zo =(X20t Xo) (X20 + Lio) 


Fig. 5—Circuit to produce the sum of two two-bit binary 
numbers. 


zation which will correspond exactly to the iterative 
circuit realization. The difficulties with designing a cir- 
cuit by factorization are that no straightforward tech- 
nique is known for determining just how to factor (this 
is still very much of an art, depending strongly on the 
skill of the designer)" and the expressions to be manip- 
ulated are extremely unwieldy. When applicable, the 
iterative design technique is useful for avoiding the 
difficulties mentioned in connection with the direct 
combinational circuit design method. 


Ill. SYNTHEsIS 


Since in an iterative network all of the cells are 
identical except for the end cells, it is only necessary to 
specify enough information for the design of a typical 
cell and for the modifications in the end cells. The typ- 
ical cells will be considered first. A (typical) cell must 
“determine” what signals to place on its output leads 
and on its carry leads to the next cell: The cell deter- 
mines this on the basis of the signals which it receives 
on its input leads and on its carry leads from the pre- 
ceding cell. The carry leads from the 7th cell to the 
t+1st cell will be called the 7 output-carries and the 
carry leads from the 7-1st cell to the zth cell will be 
called the 1 input-carries. A cell can be specified by a 
table which lists the signals on the output leads and 
the output-carry leads for each possible combination of 
signals on the input leads and the input-carry leads. 
Table II shows such a table for the binary adder of 
Fig. 1. The input lead combinations are listed across the 
top of the table and the input-carry lead combinations 
are listed along the side of the table. The corresponding 
output lead and output-carry lead conditions form the 
entries of the table. Such a table will be called an output 


uT, J. Beatson, “Minimization of components in electronic 
switching circuits,” Trans. AIEE, vol. 37, pp. 283-291; July, 1958. 
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TABLE II 
Output TABLE FOR A BINARY ADDER 


X1X2 
Ga Or On Ort, tide deo 
One $080-O9le 1000004 


1 OF 1 ASO tO 


Ve 
V =y(x2-+01) axe 
Z = x1'x2'y +1 xy! x xe"! -boxey 
= (01 +%2+y) (VY! +-a1x2y) 


table. The algebraic expressions for the outputs and 
output-carries can be determined directly from such a 
table by standard techniques.!? This table is identical 
with the output table for a clocked (synchronous) se- 
quential circuit in which the y-variables represent pres- 
ent-state variables and the Y-variables represent next- 
state variables. Thus an analogy exists between iter- 
ative networks and clocked sequential circuits, and 
many of the techniques of sequential circuit design will 
be shown to apply to iterative networks. 

An output table specifies the typical cell completely. 
For the first cell, an assumption must be made of what 
the input-carries should be. The proper assumption is 
clear from the requirements of the problem. In the case 
of the binary adder, the carry lead represents the arith- 
metic carry of ordinary addition: Clearly the carry to 
the first cell will be zero. This knowledge of what the 
carry must be can be used to simplify the design of the 
first cell. Similarly, in the last cell, the elements which 
generate the carry outputs can be eliminated since no 
carry outputs are necessary. 

It has been shown that an iterative network is de- 
termined completely by an output table for which the 
input carries to the first cell are specified. The problem 
of designing an iterative network therefore can be re- 
duced to that. of discovering an output table which 
corresponds to a network having the desired properties. 
The following material discusses the formulation of 
output tables. 

The first step in writing down an output table is to de- 
cide what “information” must be passed from cell to 
cell via the carry leads. In the binary adder the signal 
sent from the zth cell to the 7+1st cell specifies whether 
an arithmetic carry has been generated in the ith cell. 
No formal procedure can be given for determining what 
information the carries must represent. Actually, in 
some sense no formal procedure is possible; writing the 
carry-output table corresponds to specifying precisely 
what is the desired behavior of the circuit. Similarly 
no “formal technique” can be given for writing the 
table of combinations for an ordinary combinational 


2D. A. Huffman, “The synthesis of sequential switching cir- 
cuits,” J. Franklin Inst., vol. 257, pp. 161-303; March, 1954. 

18M. Phister, “Logical Design of Digital Computers,” John 
Wiley and Sons, Inc., New York, N. Y.; 1958. 
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circuit. Some consolation can be derived from the fact 
that if more carry “information” is assumed than is 
actually necessary, the superfluous “information” is 
easily eliminated by a formal procedure, Section IV. 
Once the information which the carry leads must convey 
has been decided upon, it is a straightforward process to 
fill in the table. For each combination of input and in- 
put-carry lead signals, a decision must be made as to 
the correct output and output-carry signals. The entire 
procedure of forming an output table is best illustrated 
by an example. A circuit is to be designed which has 
inputs Xn; Xn=1, « 7, Xv andsoutputs 2.) Zanes 
Z, such that Z;=1 if and only if” 


Se = atror 2 


j=1 


It must be possible to determine from the carry leads 
the following pieces of information: 


0) None of the inputs (X,;) to preceding cells are equal 
to one. 

1) One of the inputs (X,;) to a preceding cell is equal 
to one. 

2) Two of the inputs to preceding cells are equal to 
one. 

3) Three or more of the inputs to preceding cells are 
equal to one. 


By analogy with sequential circuits each such possi- 
bility will be called a carry state. Since the carry leads 
have binary signals, at least two leads are necessary to 
convey these four possibilities. At this stage of the de- 
sign process it is unwise to assume a specific number of 
carry leads or to decide which signals on the carry leads 
will represent the different carry states. Instead, each 
carry state will be assigned a decimal number. A table 
will be formed which lists the input combinations across 
the top and the input-carry states along the side. The 
output-carry states and circuit outputs will be the 
entries of the table. (See Table III.) This table is called 
a state table or flow table. If the input-carry state is 1 


. (one previous X;=1) and X;=0, the next carry state is 


again 1 and Z,=1 since there has been one X; equal to 
1. If the input-carry were 2 and X;=1, the next carry 
would be 3 and Z;=0 since then there would be three 
X;=1. The input-carry state will be represented by sy 
and the output-carry state will be represented by S;, where 
s, and S; represent decimal numbers. 

The output table is derived directly from the state 
table by replacing each decimal number by an arbitrary 
binary number.’ The binary numbers will correspond 


“The addition here is ordinary addition rather than Boolean 
addition. 

© There are many possible choices for the binary numbers. Dif- 
ferent choices will result in different networks. The problem of 
choosing the binary numbers with which to encode the states is 
discussed at the end of this section. 
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TABLE III 


TABLES FOR AN ITERATIVE CIRCUIT FOR WHICH Z;=1 IF, AND 
ONLY IF, >, §-,X;=1 or 2 


(a) State Table 


Sk 0 1 
0 preceding «=1 0 ONON 1d 
1 preceding x=1 1 Poe 
2 preceding x=1 2 Pell Weed) 
3 or more preceding «=1 3 SD 80) 
Sa 
(b) Output Table 
xX 
Viy2 
0 1 
01 Oil Mik eal Vy =x'y1 tayo 
11 dheailahe il Oka Vo=yo(x’ +41’) 
10 TOF 10) ONO) Lala 
00 00,0 00,0 


to the actual signals on the carry leads. Each occur- 
rence of the same decimal number must be replaced 
by the same binary number. The output table of Table 
III(b) results from the state table of Table III(a) when 
the replacement of decimal by binary numbers is made 
according to the following table: 


SE YiY.2 
Sk 1y2 
0 O1 
1 1B ai 
2 © 
3 0 0. 


The following algebraic expressions result from the 
output table of Table III(b): VYi=xi/yi tx, Yo=ye 
(x;’+y1'), Z;= V1. The corresponding circuit is shown in 
Fig. 6. The typical cell is designed, using standard de- 
sign techniques,'® as a combinational diode logic circuit 
having inputs x;, x;’, ¥1, ye and outputs VY, Ye, and Z,. 
It is usually possible to simplify some of the end cells 
(cells 1, 2,---and m, m—1,--- ) so that they require 
less equipment than typical cells. For the first cell, the 
input state is known. In the circuit whose design is be- 
ing described, this first input state is 0: None of the 
preceding cells have x;=1, since there are no preceding 
cells. State 0 was coded as y:=0, y2=1 so that these 
values can be substituted in place of y, and yz in the 
equations for Y%, Yo and 21: Yu=m'-0+x%,-1=%, 
Vo =1(%,'/+1)=1, Z:= Vu=x. These equations lead 
to cell 1 of Fig. 6. For cell 2, yi2= Yu, yo= Ya=1, so 


16 Caldwell, op. cit., ch. 9. 
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eS Sear Ss a ee eee, 


Zn 


Fig. 6—Network for the output table of Table III. 


that Vy = X0/ Vip toe: 1 = Vy +x, Yoo = 1(%2’ + Vu’) 
=%2'+ Vi’. Cell 3 will be a typical cell since yy3= Vi», 
yo3= Yo. For the final cell (7) it is only necessary to 
generate Z, since there is no next cell to which Y,, and 
Yon must be connected. However, Z,= Yin requires 
both y; and y. so that cell (7—1) must generate both 
Wann ennal J Sr. 

This completes the design of the network, and illus- 
trates the general design procedure. One problem con- 
nected with the design procedure which was not dis- 
cussed is the assignment problem—the assignment of 
binary numbers to the states. For gate networks, any 
choice of the binary numbers will lead to a realizable 
network. This is not generally true for relay networks. 
It is always possible to obtain a physically realizable 
relay network by assigning to the states only binary 
numbers which contain only one 1. For example, a cod- 
ing for the table of Table III(a) which would lead to a 
realizable relay network would be: 


Sx VaYoVesya 
Sk Jie yy 2 
0 1+/0,,05.0 
1 Opto 0 
2 Od0 E1250 
3 Or OE 


The resulting output table and carry functions are 
shown in Table IV and the network is shown in Fig. 7. 
For this table it turns out that 4 is not necessary since 
neither Z nor Y1, Y2 or Y3 depend on it. 

A more complete discussion of the assignment prob- 
lem will be presented elsewhere. 

Just as it is possible to design sequential circuits from 
a state diagram!” !8 rather than a state table, it is also 
possible to use a state diagram in carrying out an itera- 
tive network design. The state diagram corresponding 


17 G. H. Mealy, “A method for synthesizing sequential circuits,” 
Bell Sys. Tech. J., vol. 34, pp. 1045-1080; September, 1955. 

18 E. F. Moore, “Gedanken-Experiments on Sequential Ma- 
chines,” Automata Studies, Annals of Mathematics Studies, no. 34, 
Princeton University Press, Princeton, N. J.; pp. 129-153; 1956. 
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DERIVATION OF CARRY FUNCTIONS FOR RELAY REALIZATION 
OF STATE TABLE OF TABLE III (a) 


IRE TRANSACTIONS ON ELECTRONIC COMPUTERS 


(a) Output Table 


be 
yiyeyays 0 1 
1000 1000,0 0100, 1 
0100 0100,1 0010, 1 
0010 0010,1 0001,0 
0001 (0C01,0 0001, 0 


Vi V2V3 Vs, Z 


(b) Function 
Vi=x'y1, Yo=x'yotxyn 


/3=x'y3 tayo, Z= Vet V3 
Y2 See 
= x Zz, 
Fig. 7—A typical cell for the network for Table III. The end- 
cell connections are shown by dotted lines, 


ce) 1/1 


i 
Zi =11F AND ONLY IF oy xj =10R C4, 
yet 


Fig. 8—State diagram for an iterative circuit. 


to the state table of Table III is shown in Fig. 8. Only 
the state table approach will be discussed here. 


IV. REDUCTION OF THE NUMBER OF STATES 


It is possible in forming the state table to assume 
more carry-lead information than is actually required. 
This leads to superfluous states in the state table. The 
purpose of this section is to discuss formal procedures 
for eliminating such unnecessary states. 

One reason that it is desirable to reduce the number 
of states to a minimum is that the reduced state table 
is less cumbersome and easier to manipulate. An even 
more important reason for reduction is that the required 
number of carry leads is roughly proportional to the 
number of states. Number of carry leads >loge (number 
of states). In an electronic network it is usually neces- 
sary to provide an amplifier per cell for each carry lead. 
Reducing the number of carry leads thus reduces the 
number of amplifiers required. For both relay and elec- 
tronic networks a reduction in the number of carry leads 
reduces the number of outputs which must be generated 
by the combinational circuitry in each cell. This may 
lead to a reduction in the combinational circuitry re- 
quired; but such a reduction is not guaranteed since 


December 


reducing the number of carries may at the same time 
complicate the combinational circuitry required for 
each output carry. 

In summary, minimizing the number of states re- 
duces the number of carry lead amplifiers required in 
an electronic network and may simplify the combina- 
tional circuitry in electronic and relay networks. 

One type of superfluous state is a state which can 
never actually occur. Such a state is called an imnacces- 
sible state. It is possible to have, in a state table, an in- 
accessible state which occurs as an input-carry state, 
but which does not appear as an output-carry state and 
is not an initial state. Since the circuit does not start 
in the inaccessible state and does not go to the inacces- 
sible state from any other state, the inaccessible state 
never actually occurs and can be removed from the 
state table specification of the circuit. It is also possible 
to have an inaccessible state which can be reached from 
other states, but cannot be reached starting from the 
specified initial state. For example, in the following flow 
table, states 3 and 4 are inaccessible if the specified 
initial state is 1. 


Xx 


s. 


INITIAL STATE 


OS) hoe ey bY 


0 
1 
2 
Ss 
4 


OQ e RF Ne 


It is also possible to reduce the number of states by 
merging groups of equivalent states.!?:!7:18 Two states 
of a circuit A and B are said to be equivalent if it is 
possible to replace A by B, eliminating A, without 
changing the performance of the circuit. The following 
experimental procedure can be carried out to determine 
whether two states A and B are equivalent. Two copies 
of the circuit being considered are constructed. The 
input carry leads to a typical cell, 7, are disconnected 


from the preceding cell, and signals corresponding to 


state A are placed in the input carry leads to cell z of 
one circuit, and signals representing state B are placed 
on the corresponding leads of the other circuit. All pos- 
sible inputs are applied to the two circuits and the out- 
puts are compared. If the two circuits always have 
identical outputs when they have identical inputs, 
states A and B are equivalent; otherwise they are not. 
The outcome of such an experiment may depend on the 
number of cells following typical cell 7. For two states 
which are “unconditionally” equivalent, the outputs 
from the two circuits used in the experiment will be 
identical for any number of cells following cell 7. Two 
states which produce identical outputs for only a re- 
stricted number of cells following cell z are said to be 
equivalent for an experiment of length J, where / equals 
the number of cells having identical outputs. Two states 
which are identical for an experiment of length / can be 


1958 


_combined only for cells which are within / cells of the 
final cell of the circuit. Combining such cells is one 
_method of determining the end cell modifications for 
_ the final cells of the circuit. 

Of course, it is not necessary to actually carry out 
such experiments in order to determine which states are 
| equivalent. The outcome of such an experiment can be 
| determined by inspection of the state table. For ex- 
| ample, in the table below, states 1 and 2 are equivalent 
| since they have the same outputs and go to identical 
output states. States 3 and 4 are not equivalent since 
they have different outputs. The definition of equiva- 
lent states given here is identical to that used for se- 
quential circuits. Equivalent states will not be treated 
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further since detailed discussions have been given else- 
where -?:17—18, 


P;* 0 1 
iGreen 
De ho dene 
Sonat 4,0 
ee ee ae 


A general method for designing iterative networks 
has been presented. Subsequent papers will present 
techniques for obtaining economical networks using 
various criteria of economy. 


Some Properties of Boolean Equations* 
N. ROUCHE} 


Summary—Solubility conditions for a set of Boolean equations 
are established, first with respect to one variable, then with respect to 
all variables. By consideration of relations between minimal terms, 
a simple matrix form is deduced for Boolean equations. Using finite 
group theory and the properties of the characteristic equation of the 
matrix, a classification is introduced for Boolean mappings and their 
iterations, to which corresponds a classification of sequential ma- 
chines. 


INTRODUCTION 


HE problem of solving Boolean equations has 
eens been dealt with in the literature, but it 
does not seem that results of great generality have 
been obtained. See, for instance, a recent paper by 
Zemanek! which makes extensive use of truth-tables. 
It appears, however, that there exists a general criterion 
for solubility, expressible in terms of a determinant of 
ordinary algebra. Furthermore, there exists a matrix 
representation of Boolean equations which is free from 
difficulties previously encountered concerning multi- 
plication and iteration of matrices.” It will be shown how 
the theory developed here throws some light on the 
problem of the structure of sequential machines. 


SOLUBILITY OF A SYSTEM OF BOOLEAN EQUATIONS 
WITH RESPECT TO ONE VARIABLE 


Let 


Sine gi, EO eo En) oy) k) (1) 


> 


* Manuscript received by the PGEC, June 10, 1958. ; 

+ Lovanium University, Leopoldville, Belgian Congo, Africa. 

1H. Zemanek, “Die Lésung von Gleichungen in der Schaltalge- 
bra,” Arch. elect. Ubertragung, vol. 12, pp. 35-44; January, 1958. 

2 J. O. Campeau, “The synthesis and analysis of digitial systems 
by Boolean matrices,” IRE TRANs. ON ELECTRONIC COMPUTERS, 
vol. EC-6, pp. 231-241; December, 1957. 


be a system of k Boolean equations defining the k 
Boolean variables 4, yo,-+-, yz in terms of the n 
Boolean variables x1, x2, +++, x,. It is well known that 
(1) may be written as follows: 


Vi = Xibir + 21’ di2 (2) 


where + indicates union, the absence of symbol indi- 
cates intersection, and x,’ is x, complemented. Both 
ga and $i are Boolean functions of x2, x3, + + + Xn. 

We now consider under what condition it will be 
possible to express x; as a function of the y,. 


The expression in the right member, reduced to canon- 
ical form, yields 


ms >» mi(Y1, V2, * 
J 


V1 = f(y, V2) ee 


S 2 Yn), (4) 


the right member being the sum of a certain number of 


minimal terms of the set of variables 1, yo, +--+, Vp. 
But making use of (2), we have 
m(Y1, V2, °° * 5 Ye) 
= xymj(ou, a1, ++ + , ber) + X1'm;(b1, $22, >> + , dre). 
Substituting into (4), we obtain 
f= 41D) mi($ur, $21, +++ 5 bus) 
+ ay! >) mj(P12, $2, °° +, Pk2) 5 (S) 


>) being in (4) and in (5) a summation on minimal 
terms of the same ranks. Therefore, we may rewrite (5) 


f = if (b11, 621, * + * » ber) + X1'f(b12, 22, °° > , bee). (6) 


Considering (6), we deduce the following theorem: 
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Theorem 1: In order that x; be a given function f(y, 
Yo, ++ +, VR) it ts necessary and sufficient that 


f(b11, 621, °° * bei) = 1 f(biz, G22, - > + » bea) = O. (7) 


The system (1) may also be written in the following 
form: 


grt 
Si = am amy, VEE Gey Os bu) 
j=1 


Qt 
+xi( Dy miler, * 


Pas 1 


£6 *»)bu) (a aa Lf 2, e255 -,k) (8) 


where each coefficient 5;; is 0 or 1, and mj42”1(a2, + -- , 
Ky) =I; (Ho, * +, Xp) tor all7<2"4, The minimal term of 
rank J of the set (1, ye, --- , Yx) may be written, ac- 
cording to (8), 


milyiy V2, °° * 5 Ve) 
gr 

= al Dd) milxe, +++, en)mi(bij, ++ bn) | 
om 


+ wi] Ss M;(%X2, iy ae) Xn)mi(b13, eae par: in) |. (9) 


ona es 1 


Any function of 1, ye,--:, ¥x, being a sum of some 
minimal terms of the form (9), admits of the following 
development: 


f(y, Sys a sea Yn) 


gr-l 


= nf 5, mj( 2X2, eae) Hn) f (Big, mae aC bu) | 


Qn 
+p ay | MS; 
y= 


os) 


PACTS ers nrg caer in) | aah) 


The necessary and sufficient conditions for 


a1 = flys, Yo, * + + 5 Ye) 
are obviously 
on-1 
Ds m;j(x2, ete ah Oipa- ae tg Op, at 
sate (11) 
DEN PH ys, UR) 0 
jan" 4d 


and these two relations must be true, whatever the 
values given to %2, +: , X,. Therefore we may state the 
following theorem: 

Theorem 2: In order that x, be a given function f(y, 


Yo, + + +, Vx), 2b ts necessary and sufficient that 

fbig, = = Dee ek 1 oe 
jaueee bp re ae) 

J Oasys OY Eee 


This set of conditions, compared with those of 
Theorem 1, has the advantage that it depends only on 
the constant coefficients b;; and not on the quantities «;. 

It is noteworthy that in the present problem the 
minimal terms of the sets of variables 
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bij, ba;, cA US ars bn; 1 Sp < 2” 


play an important role. We shall see shortly that it is 
rewarding to dispose these minimal terms in an ordered 
rectangular array according to the following rules: ali 
the minimal terms of the set of subscript 7 are written 
under each other in a column. These columns are juxta- 
posed, from left to right, in order of increasing j. In 
each column the minimal terms are arranged down- 
wards in the order suggested by the increasing se- 
quence of the first 2" binary numbers. For instance, the 
minimal terms of three variables a, 6, c, would be writ- 
ten downwards in the following order: 


a'b’c’, a'b’c, a’bc’, a’bc, ab’c’, ab’c, abe’, abc. 
The rectangular array is characterized by the fact that 
it contains one 1 per column, all other terms being 0. 
We shall call it the matrix of the Boolean mapping (1). 


This term will be justified later. We represent the matrix 
by 


B = [61)]. 
In the case k=n=2, the equations of the mapping are 


yr = byyxy' x0! + Dion’ x2 + bizseixe’ + ee (13) 


— if if 
Nb pa boxy’ x2" + booxy' xe a 230149! + boats 


and 
bii/bo1’ — By0' Bon’ B3'b23’ by 4’ Boa" 
hens biv'be1 DB 12’bo2 OB 3'b23 Ba’ a4 (14) 
birder’ ~—2be0’ ~—Bb23 BD yadoa’ 
biib21 ~~ Byebe2 ~~ 3b23 — yado4 


THE EQUIVALENT MAPPING ON MINIMAL TERMS 


Momentarily, in order to prove the following theorem, 
we make use of a matrix similar to the A matrix intro- 


duced by J. O. Campeau? 

Pee As 25. 

A= fa,| 
ee a ee Oe 


which has k rows and 2* columns. The elements of 


column 7 are uniquely defined by 
4 = jhe a2" * + Goj2" +.-.-. + G52”. (15) 


In the case k=2, we have 


E Oatky 
A= if 
Oa tes ORaee 


Now, the 0;;’s of (8) are given by 
big = Do anby; 


U « 


and therefore 
Jo = > Se a1 15; 
ay 


if &; is the minimal term of rank 7 of the set (x, -- - 


| 
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x,). In the same way, we obtain 
Ve. = 3 Dai Buti. 
Gat 


A typical minimal term of the set (1, - +: , Vx) is 


Diya! Woot °° Me 


iS (3 Dy onda) ( Dye oa bi) its or 


Delt Ton Bb 


( ye De ond) ( > 3, oot top Bis) rAd 


i) lg j lg +1 


( Paap ola, 3) . 
I 


j 


This may also be written 


“yy! ut VaNgit 2 > Vb 


ms 2 


11,0’ 21, Ps Batylg Ling a) CP POS TOR GIS 
Uy ,12, +++, lg lg4i,++ ile 
( De 81,8) ( »y 8198) zie 
( X Pui) ( 2 Bist) or ( aL, anti) (16) 


Among all products of a, in this expression, the only 
one which is different from 0 (and equal to 1) corre- 


sponds to 
R=t= he], SH tH =] 


where / is the rank of the minimal term which is com- 
puted. Therefore (16) may be reduced to the very 
simple form 


= i 
1 x Bik; (17) 
Using matrix and vector notations, we write the formula 
CBE: (18) 


This relation shows why we have called B a matrix. 
Using the rules of Boolean algebra, we have shown that 
(18) can be deduced from (1) or (8). Conversely, it is a 
very simple matter to show that (1) or (8) can be de- 
duced from (18). We have 


Gagk = ye) AL 
or, using (18) we have ‘ 
Gi,¥2, 5 y= ABE. (19) 
It is readily seen that 
bir Die peal 
AB= (20) 


Dia Duos + eit 
and therefore (19) is equivalent to (1). 


Theorem 3: To each mapping of n Boolean variables 
X1,° °°, Xn, to k variables yi, ++, Ye corresponds an 
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equivalent mapping of the 2” minimal terms of the x;’s 
into the 2" minimal terms of the variables y;. This equiva- 
lent mapping may be represented in matrix form, the 
terms of the matrix being simple combinations of the 
coefficients in the equations of the original mapping. 
Any Boolean mapping may be studied on the equivalent 
mapping on minimal terms. It is noteworthy that (18) 
may be interpreted in ordinary algebra as well as in 
Boolean algebra. This is due to the fact that the vectors 
of minimal terms never have more than one component 
equal to 1, all others being 0. Therefore, the only typ- 
ical Boolean relation 1+1=1 never appears in (18). 
We shall make use of this fact later. 


Br-UNIFORM MAPPINGS OF ” VARIABLES 
INTO 2 VARIABLES 


In this case, the matrix B is a square matrix. For the 
mapping to be bi-uniform, it is necessary and sufficient 
that the correspondence between minimal terms of the 
two sets of variables be a one-one correspondence. We 
therefore obtain the following theorem. 

Theorem 4: A system of Boolean equations like (1) 
with k=n may be solved for x1, %2, + --* , Xn Uf and only if 
the corresponding matrix B is nonsingular in the sense of 
ordinary matrix theory. 

In this case, B contains one 1 per row and per column. 
It is a permutation matrix, and Theorem 5 results. 

Theorem 5: Any bi-uniform mapping of n variables 
X1, X2, °° * y Xn, into n variables y1, yo, + + + 5 Yn 1S equiva- 
lent to the equality of the vector of minimal terms in 1, Yo, 

- ++, y, and the vector obtained by a certain permuta- 
tion of the components of the vector of minimal terms in 
Ma. Ons 

As the vector of minimal terms has 2” components, 
we may also state the following: 

Theorem 6: There exist exactly (2")! bi-uniform map- 
pings of n variables into n variables. 

An interesting necessary condition of bi-uniformity is 
suggested by the form of B. For a system like (1), with 
k=n to be solvable with respect to x1, %2, °° * , %n, it is 
necessary that all the right members of the equations 
consist of 2”~! minimal terms. 

Obviously the set of all matrices corresponding to bi- 
uniform mappings of variables into 2 variables forms 
a group of order (2”)!, isomorphic with the symmetric 
group of degree 2”. We shall call this group G. 


SUBGROUPS OF THE GROUP G OF 
Bi-UNIFORM MAPPINGS 


At this point, it will be interesting to survey some of 
the subgroups of G to obtain a classification of all pos- 
sible bi-uniform mappings. 

Let us call centro-symmetric every permutation 
matrix which is symmetric with respect to its center. In 
order to deal more easily with centro-symmetric ma- 
trices, we momentarily introduce a new convention con- 
cerning subscripts of rows and columns for permutation 
matrices. If we are concerned with a matrix of even 
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order n?, the subscripts 7 for the rows will be all the odd 
multiples of 0.5 taken in increasing order and satisfying 


n—1 n—1 


(21) 


2 2 


A similar convention is adopted for the column sub- 
scripts. The following is an example of a centro-sym- 
metric permutation matrix of order 8? with the sub- 
scripts explicitly indicated. 


a5 02 — De Ok — 1 OS OnO TOO Iu. 0) home 
=.) 1 0 0 0 0 0 0 0 
—2.5 0 0 0 0 1 0 0 0 
— [9 0 1 0 0 0 0 0 0 
= 059 0 0 0 0 0 1 0 0 
0.5 0 0 1 0 0 0 0 0 
ES, 0 0 0 0 0 0 1 0 
Zio 0 0 0 1 0 0 0 0 
3.5 0 0 0 0 0 0 0 1 


In the case of a matrix of odd order n?, the subscripts 
will be all the integers, including zero, belonging to the 
closed interval determined by (21). 

The condition for a matrix [a;,] to be centro-sym- 
metric is then 


(22) 


Gipsy Fe 


The condition for symmetry (in the usual sense, 7.e., 
with respect to principal diagonal) is 


(23) 


ai; = Azz. 


The condition for secondary symmetry (i.e., symmetry 
with respect to the other diagonal) is 


(24) 


Cig = Oui 


Eqs. [(22)-(24)] show that any symmetric centro- 
symmetric matrix displays secondary symmetry. 
Theorem 7: The set of centro-symmetric permutation 


matrices of order n® is a subgroup of the group of permu- ~ 


tation matrices of the same order. 

The only axiom to be verified is the axiom of closure. 
Suppose that A and B are centro-symmetric. Let us 
prove that C= AB is centro-symmetric. We have 


(n—1) /2 
kt = AD; 
j=—(n—1)/2 


(25) 


and 


(n—1) {2 
Ck,-1 = Gp, 505,-1. 
j=—(n—1)/2 


(26) 


But since A and B are centro-symmetric, (25) yields 


(n—1) /2 


Ca = bye 


J=— 1) [2 


(27) 


0_1,—jb_j,-1. 
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The summation being commutative and the interval 
covered by the summation subscript being symmetric 
with respect to 0, we may change —j into j in (27) 
which gives 

(n—1) /2 


(28) 


Cn OPT 
j=—(n—1) /2 
Comparing (26) and (28), we get C_x, -1= Cir. We shall 
call this subgroup Gi. It is readily verified on an exam- 
ple that G; is in general nonabelian. 

Theorem 8: The order of G, is N=2"(n/2)! To con- 
struct all centro-symmetric matrices of order n’, we 
first put a 1 in the first row. There are 7 possibilities. 
At that moment, because of the condition of centro- 
symmetry, the last row is already occupied. There re- 
mains ~—2 possibilities for the second row, etc. Thus 


N = n(n — 2)(n — 4) +++ 4.2 
N -(2 Ne ) 
= — —— 1 —-—2 a 
Qnl2 DS D) D, 
») 
N = 2n/2{—) 1}, 
c 


In the study of Boolean mappings, we are interested in 
matrices of order 22”. For such matrices 


N = 27" (2"-»)!, 


Ck = 


then 


and 


(29) 


(30) 


Let us consider again the group G of all permutation 
matrices of order 27" corresponding to the permutations 
of a vector of minimal terms with 2” components. Re- 
versing the complementation of certain variables in the 
minimal terms is equivalent to permuting them in a cer- 
tain manner. Thus G contains a certain number of 
matrices whose application to a vector of minimal terms 
merely reverses the complementation of some variables. 
We shall call these matrices, matrices of complementa- 
tion reversals. 

Theorem 9: The matrices of complementation reversals 
form an abelian subgroup Gy of order 2” of G,. The order 
that we have adopted for minimal terms in a vector of 
minimal terms is such that the following property holds. 
Two minimal terms in a symmetric position with re- 
spect to the center of the vector may be obtained from 
each other by complementation of all variables. And, 
conversely, if we reverse the complementation of some 
variables in all the minimal terms, the same property 
obviously still holds. Thus, all matrices of complemen- 
tation reversals are centro-symmetric. 

The axiom of closure is satisfied because it is clear 
that the product of two complementation reversals is a 
complementation reversal. There are 2” distinct ways 
of reversing some complementations; thus the order of 
the subgroup is 2". Reversals of complementations are 
commutative, 7.e., the subgroup is abelian. The product 
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of two identical reversals of complementation is the 
identical operation. Thus, every matrix of Gy is also 
‘symmetric, and therefore secondary symmetric. This 
| also shows that Gy is generated by 7 subgroups of order 
i. 

One of the characteristics of the matrices of Gy is 
that, except for the identity matrix, the elements at 
both ends of the principal diagonal are always 0, the 
| first and the last minimal terms of a vector changing in 
every complementation reversal. 

Let us now consider another type of permutation 
/ matrix. A certain number of the matrices of G corre- 
spond to a permutation of the original Boolean varia- 
bles. There are m! such permutations, which are equiv- 
alent to changing the names of the variables, and do 
not alter the structural character of the mapping. The 
matrices of this type will be called variables permutation 
matrices. 

Theorem 10: The variables permutation matrices con- 
tained in G form a nonabelian subgroup Gy of order nl of 
 G,. First, it is obvious that these matrices form a sub- 

group of G because the product of two permutations of 
the variables is a permutation of the variables. Further- 
more, this subgroup is isomorphic with the symmetric 
group of degree », and therefore is nonabelian. 

Gy is a subgroup of G,;. All variables permutation 
matrices are centro-symmetric because two minimal 
terms in a symmetric situation with respect to the center 
of the vector are transformed by a variables permuta- 
tion into two minimal terms having the same property. 

One of the characteristics of the matrices of Gy: is that 
the elements at both ends of the principal diagonal are 
always 1, since the first and the last minimal terms of a 
vector never change in a variables permutation. This 
property and the fact that Gy; has the inverse property 
prove that Gy and Gi have but a single element in 
common. 

Theorem 11: The product Gz=GuGt of the complemen- 
tation reversals group and the variables permutation group 

is a subgroup of order 2"n! of the centro-symmetric group. 
For all the subgroups so far identified in G, we compute 
Table I. This shows that for »=1 and m=2, G, is an im- 
proper subgroup of G,, and for n>2, G2 is a proper sub- 


group. 


TABLE I 
Number of os 1 2 3 4 
variables 
Order of G Oo” )! 2 | 24 | 40, 320 | 20, 924, 789, 888, 000 
Order of G; | 27° 74(27-1)/ | 2 8 384 10, 320, 920 
Order of Gz | 2°n!/ 2 8 48 384 
Order of Gu | 2” 2 4 8 16 
Order of Giz | 1! 1 2 6 24 


It is possible to develop still further the theory of 
permutation matrices and of singular Boolean mappings 
square and rectangular matrices from an abstract point 
of view, but it will be more convenient for our purpose 
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to go in on this direction after having introduced the 
notion of sequential machine which will constitute a 
suggestive example of application of Boolean mappings 
to a concrete case. 


AUTONOMOUS SEQUENTIAL MACHINES 
AND FLOW-GRAPHS 


We consider a machine with m terminals, each of them 
being capable of two states (for instance low and high 
potentials). The machine receives a discrete sequence 
of advancing pulses, and it may change the states of its 
terminals at each advancing pulse. If none of the 
terminals is used for arbitrary input signals introduced 
before the advancing pulses, we call such a machine an 
autonomous sequential machine. As has already been 
shown,? the y,;’s of (1) with k= may be interpreted as 
giving the state of the machine after advancing pulse f, 
if the x,’s give the state of the machine before the same 
pulse. All the properties of the machine are then con- 
tained in the matrix B of the mapping. The properties 
of the matrix may be visualized in a flow-graph, each 
state of the machine being represented by a small circle 
containing the number of the state, that is, the rank of 
the corresponding minimal term. If, before an advancing 
pulse, the machine is in state V, its state after the ad- 
vancing pulse will correspond to the row subscript of 
the single 1 of column V in the matrix.® 

For instance, suppose a machine is defined by ihe 
following: 


/ 5 Pas EO 2 A 
Yi = X102%3 + Xox3 + Xa x3 


Yo = y'X9'H3' + X1x3' + x1'Hst 5 (31) 
ys = (x1'x2' + 41%2) x3 
its matrix 1s 
0 
B 
| (32) 


OS Re OS oo Oo: © 

SS a SSS ale = CORSO 
Se Si Oy a) Or OR oe 
Gy en SO Se OO a * 
Gy me. oO Se Orr 
oS). © eo? OF G'S 
(oly Me Eth =? ila tegen Hs) 
So, OS eS: 1 o"OR Or iO 


and its flow-graph is shown in Fig. 1. 


REVERSIBLE AUTONOMOUS SEQUENTIAL MACHINE 


We call reversible every machine whose matrix B is 
nonsingular. To every 1 situated on the principal di- 
agonal, there corresponds a distinct part of the flow- 


It will be recognized here that the matrix B is similar to the 
transition matrix defined in the following reference. 1 E. Hohn, 
S. Seshu, and D. D. Aufenkamp, “The theory of nets,” IRE TRANS. 
on ELecTRoNic Computers, vol. EC-6 »pp. 154-161; September, 
£9545 
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Fig. 1—Flow-graph of a machine. 


Cm) 


Fig. 2—One—state circle. 


graph consisting of a single state circle, as shown in 
Fig. 2. It is possible to find a permutation matrix P 
which will bring up this 1 into leading position on the 
principal diagonal if B is transformed by P. The matrix 
B is then transformed to 


B, = PBP- (33) 


If there is another 1 on the principal diagonal, a trans- 
formation of B, by a permutation matrix will bring this 
1 into second position on the principal diagonal, and so 
on. When all the 1’s of the principal diagonal are ex- 
hausted, suppose there is in the matrix a nonsingular 
principal minor of order 2. Transformation by an ap- 
propriate permutation matrix will bring this minor next 
to the line of 1’s already treated. This minor corresponds 
to a double state circle in the flow-graph, as shown in 
Pig..3; 

More generally, we rearrange the matrix as a succes- 
sion of nonsingular principal minors, as shown in Fig. 4, 
each principal minor containing no nonsingular minors. 
Each separated principal minor of order & in this repre- 
sentation has a corresponding k-state circle in the flow- 
graph. 

To every such principal minor are associated k roots 
of the characteristic equation of the matrix, and these 
roots are given by 


(—A)* -—-1=0 
if k is even, or by 
(Nhe l = 0 
if k is odd; that is, in both cases, the k roots are 


Pri 1k) 


Nee deri ae 


From this we deduce Theorem 12. 

Theorem 12: The multiplicity of the root P™1'* (where 
k is an integer) of the characteristic equation of B is equal 
to the number of circles of the flow graph whose degrees 
admit k as a divisor. 

In particular, the multiplicity of the root 1 is equal to 
the number of distinct circles of the flow graph. 

We shall call simply-connected any machine for which 
the multiplicity of the root 1 is equal to 1. 
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Fig. 3—Two—states circle. 


Fig. 4—Ordered permutation matrix. 


IRREVERSIBLE AUTONOMOUS SEQUENTIAL MACHINES 


Theorem 13: Every matrix B has at least one principal 
minor different from 0. 

The matrix B is characterized by the fact that it has 
one and only one 1 per column. Suppose the determi- 
nant of B is 0. Then there is at least one row of zeros. 
Let us remove from the matrix this row and the column 
of same rank. We are left with a principal minor. If this 
minor is 0, it contains at least one row of zeros. We re- 
move this row and the corresponding column, etc. If all 
principal minors found by this procedure are 0, we are 
at last left with a principal minor containing only one 
element, and this element is 1. 

Let us now look for an ordered form for B in case it is 
singular. We choose one of the nonzero principal minors 
of least order, say Ao, and transforming B by a permuta- 
tion matrix, bring this minor into leading position on 
the principal diagonal (Fig. 5). If there is no 1 in the 
rectangle A; we bring the next nonzero principal minor 
into second position on principal diagonal. If A; contains 


- a certain number of 1’s we bring them close to Ao in Az 


using a permutation matrix (Fig. 6). If there is no 1 in 
A;3, we bring the next nonzero principal minor with its 
upper left corner at point x. Otherwise, using a permu- 
tation matrix, we bring all the 1’s of A; in the left part 
of A; that we call A,, etc. 

At last, we are left with a matrix B, ordered as shown 
in Fig. 7, where each square on the principal diagonal 
is as shown in Fig. 8. 

It is clear that ordering the matrix B has not altered 
its flow graph. Theorem 14, which follows, remains true 
for irreversible machines. The number of rows of zeros 
in B is equal to the number of transients that the 
machine may exhibit, z.e., to the number of states that 
may be adopted as initial states in the operation of the 
machine, but will never be attained in course of opera- 


Fig. 6 


tion. It is evident from Fig. 8 that the number of rows 
of zeros in B is equal to the multiplicity of the root 0 
of the characteristic equation of B, and Theorem 14 
results. 

Theorem 14: The multiplicity of the root 0 in the char- 
acteristic equations of B is equal to the number of possible 
distinct transients the machine may exhibit. 


ISOMORPHISM OF MACHINES 


As has been shown, pre or post-multiplication of B 
by a matrix of the Group G: only changes the names or 
complementation of the variables but does not alter the 
structural character of the Boolean mapping. Thus, ex- 

cept for some possible inverters at inputs and outputs, 

the circuit may be the same for B and for B multiplied 
by a matrix of G2. We shall call isomorphic in a restricted 
sense two machines whose matrices B; and B, may be 
obtained from each other by pre or post-multiplication 
by a matrix of G). Two such machines may be con- 
structed with the same arrangement of elementary 
logical circuits. 

On the other hand, transformation of B by a permu- 
tation matrix does not alter the character of the cor- 
responding flow graph. We shall call zsomorphic in a 
broad sense two machines whose matrices B, and B,; may 
be obtained from each other through transformation by 
a permutation matrix. Two such machines will have the 
same flow graph, except for the names of the states. On 
the other hand, they may correspond to two altogether 
different Boolean mappings. 

For instance, the ordered form of the matrix (32) is 
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Fig. 8—Typical principal minor of ordered singular matrix. 


(34) 


ay ka oa ee (Se = 
S OSS S&S] iS SS) 
SS Qos eS SiS 
Sh SS] ey Sa es) SS 
Ses) ha eS) SS) 
Se eta S&S] 2 > 
SS OS 1S SSS 
ay er Sey =) eer eS) 


The corresponding equations are 


V1 = X1%3 


(35) 


Se 1 , 
Yo = Hy Xo'x3 + 2X yXox3! 
a Vhs 
Y3 = XyX2%3' + Xoxg + x0! Hs" | 


and the flow graph is shown in Fig. 9. 
Thus we may say that, to some extent,* the problem 
of minimizing a sequential machine whose flow graph 


* The minimization of a sequential machine may be considered 
as consisting of three steps. 

1) Minimization of the flow graph by reduction of equivalent 
states. (See D. D. Aufenkamp’s and F. E. Hohn’s, “Analysis of 
sequential machines,” IRE Trans. ON ELECTRONIC COMPUTERS, 
vol. EC-6, pp. 276-285; December, 1957.) 

2) Determination among the possible sets of Boolean equations 
of the set which lends itself to the best minimization. 

3) Actual minimization of the set of Boolean equations. 
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Fig. 9—Flow-graph of machine defined by matrix (34). 


is given, is reduced to that of finding the permutation 
matrix P which, transforming B, will yield another 
matrix B corresponding to the simplest possible set of 
Boolean equations. We shall make no attempt here to 
solve this minimization problem. 


SEQUENTIAL MACHINES WITH EXTERNAL INPUT 


We take as representing the equation of such a 
machine the following set of Boolean equations: 


pA eS files, Sty May, May t1, °° * he) 


(36) 
Vin Fy (1, a 


7) Unyy Vny41, * * * 9 Ung 


where x; (4=1,---, m) are the states of the output 
terminals before advancing pulse number m, and 
x; (¢=m+1,---, m) are the states of the input 
terminals, established before advancing pulse number 
n. The y,’s are the state of the output terminals after 
advancing pulse number zn. 

Eq. (36) may also be written 


271 
it DD mila, worst peg) es en pots ears > Xng) | 
—s 
(37) 
PASE 
By 7 es mim, sis tee) Wa) Lins Saret PS oo) Xing) 
7 


With (37) the machine appears similar to an autono- 
mous sequential machine with the important difference 
that the flow-graph may be changed by exterior signals 
before each advancing pulse. 

Eq. (36) may also be written 


2" 
Sine DS MU Lees 3 


t=1 


One Xn.) M(x, cae) ny) 
(38) 


"e 
Mi Do mil nH, sis 


i=1 


: Kay aC, 6K 


a) ni) 
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With (38), the machine appears to give a definite re- 
sponse to a given excitation except that it is capable of 
changing its response by taking the previous response 
into account. Both types of machines are defined by (36) 
and correspond to the same machine as seen from two 
different viewpoints. 

Considering (37), we may define a matrix B as for the 
case of an autonomous machine. But now, B will oc- 
casionally change before each advancing pulse. The 
operation of an autonomous machine was given by the 
iteration of the matrix B. More precisely we had 


E(m + k) = BYé(n) 


where é is the vector of minimal terms of the terminals, 
and the argument of £ is the rank of the advancing pulse 
(z.e., the time). In the present case, we have 


En + k) = Ban +k —1)--- Bin)é(n) 


and B has to be given as function of time. 
Suppose this function is periodic of period k. Let 


By= BQH A) BG (39) 


Then 
E(n+ 7k) = Bi(n) (7 =1,2,---). 


Theorem 15: The operation of a sequential machine 
with a periodic input of period k, when it is observed at 
every kth advancing pulse, 1s the same as that of an auton- 
omous machine with matrix B=B,. 

The complete operation of the machine is given by 
the matrices B,, Bz,---, By, corresponding to the k 
circular permutations of the terms in the second member 
of (39). 


CONCLUSION 


It has been shown how every set of Boolean equations 
corresponds to another set of equations in which the 
variables are the minimal terms of the variables in the 
initial system. This equivalent system of equations 
lends itself very easily to a matrix representation, the 
characteristic equation of the matrix giving important 
informations concerning the structure of the Boolean 
mapping considered and its iteration. 

It seems likely that the methods developed in the 
present paper will lead to more results concerning classi- 
fication, analysis, and synthesis of sequential machines. 


1958 
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Analysis of Sequential Machines II* 


D. D. AUFENKAMP} 


Summary—Mealy’s model of a sequential machine is assumed 
and a relation of ‘‘compatibility” of states is introduced to further 
the analysis of such machines. In the event that input restrictions 
exist it is often possible to effect combinations of states under this 
relation in addition to those permitted under equivalence of states, 
a relation previously studied. Compatibility of states is analyzed by 
an iterative technique, rigorously established, which makes it possi- 
ble to determine readily connection matrices of simpler ‘‘compati- 
ble”? machines. 


I. INTRODUCTION 


alence” of states in Mealy’s model of a sequential 
machine. That analysis yielded an iterative tech- 
nique for finding all sets of equivalent states. 

For machines in which there are input restrictions it 
is often possible to effect combinations of states in addi- 
tion to those possible under equivalence without alter- 
ing the relationships between inputs and outputs, as 
long as only allowable input sequences are presented to 
the machine. In this paper we investigate the relation 
of “compatibility” of states which permits this addi- 
tional simplification of state diagrams. The notation and 
conventions established in the earlier paper! will also be 
observed in the following discussion. The problem dis- 
cussed here has also been studied by Huffman.? 

We recall the definition of Mealy’s model of a sequen- 
tial machine and certain other definitions which play an 
important role in the subsequent development. 

Mealy’s model of a sequential machine* assumes that 
the machine may be described in terms of: 

1) a finite number of states, 51, 52,°-° +, Sn; 

2) a finite number of distinct inputs symbolized by 


iz a previous paper,! we investigated the “equiv- 


M1, X2, Sa ’ Xmy 
3) a finite number of distinct outputs symbolized by 
EVID 2 3 ki eas 


It is assumed that the present output and the next 
state are uniquely determined by the present state and 
the present input. Each such machine M may be repre- 
sented by a net*® called the state diagram of M in the 


* Manuscript received by the PGEC, June 20, 1958. This work 
has been carried out under the Lockheed General Res. Program. 

+ Missile Systems Div., Lockheed Aircraft Corp., Palo Alto, 
Calif. P 

1D. D. Aufenkamp and F. E. Hohn, “Analysis of sequential 
machines,” IRE TRANS. ON ELECTRONIC CoMPUTERS, vol. EC-6, pp. 
276-285; December, 1957. nhs ee. 

2D. A. Huffman, “The synthesis of sequential switching circuits, 
J. Franklin Inst., vol. 257, pp. 161-190, March, 1954; pp. 275-303, 
April, 1954. ‘ eee 

3G. H. Mealy, “A method for synthesizing sequential circuits, 
Bell Sys. Tech. J., vol. 34, pp. 1045-1079; September, 1955. 

4F.E. Hohn, S. Seshu, and D. D. Aufenkamp, “The theory of 
nets,” IRE Trans. ON ELEctRonic Computers, vol. EC-6, pp. 
154-161; September, 1957. 

5 Mealy, op. cit., p. 1052. 


following way: the set of vertices vj, v2, +--+, U, of the 
net corresponds in 1-1 fashion to the set of states 51, 
So, + * +, Sp Of M. For each input which effects a transi- 
tion of the machine from state s; to state s;, there is a 
direct branch from the vertex v; to the vertex v;. The 
ordered pair consisting of the input symbol concerned 
and the symbol for the corresponding output is assigned 
to this branch as its weight. When there is more than 
one such branch, say k branches from a vertex v; to a 
vertex v;, we simplify matters by employing only one 
branch and assign it a weight (x,,. Yr) V (Xre, rg) Vt 
V/ (Xr. rz) (where \/ means “or”) which defines the 
union of the weights of all the original branches from 
U; tO U;. 

The connection matrix of the machine is defined as 
follows :°:7 


C = [cis] 


where c;; is the above defined union of the branch 
weights on all branches from vertex 7 to vertex j. If no 
branches connect vertex 7 to vertex j, then c,;;=0. 

By an input sequence we mean a sequence of input 
symbols, for example, x:, x,,---, x; We say that a 
machine has input restrictions® if its output and next 
state are not defined for every input. Any input sequence 
which, if applied to the machine when it is in state 1, 
violates no input restriction of state 7 or of any subse- 
quent state, is called an allowable input sequence for 
state 1.° 


II. FUNDAMENTAL DEFINITIONS 
Definition 1 


Two states, s; of a machine M and t; of a machine JN, 
are called compatible states if and only if, starting with 
M in S;and N in t;, the two machines will yield identical 
output sequences on being presented with any given in- 
put sequence from the intersection of the set of allowa- 
ble input sequences associated with state s; and the set 
associated with state ¢;. The output sequences associ- 
ated with distinct input sequences need not, of course, 
be the same. An important special case occurs when the 
intersection is the null set, 7z.e., when the two states 
have no common allowable input sequences. At the 
other extreme, if the two sets of allowable input se- 
quences are identical, then compatibility of states be- 
comes equivalence of states.® 


6 Hohn, Seshu, and Aufenkamp, op. cit., pp. 156, 158. 
7 Aufenkamp and Hohn, op. cit., p. 276. 

8 Tbid., p. 277. 

9 Ibid., p. 278. 
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Definition 2 

Two machines M and WN are called compatible ma- 
chines if and only if for each state s; of M there exists 
at least one compatible state ¢; of N, and for each state 
t; of N there exists at least one compatible state s; of M. 

Except in special cases, compatibility of states and of 
machines are not equivalence relations. The relations 
are reflexive and symmetric but they are not, in general, 
transitive. 


Definition 3 

A set of states, S= {s1, Gannon EA se} of a machine 
is called a set of pseudoequivalent states if and only if 
states s; and s; of S are compatible for all 7 and J, 2.e., 
if and only if the relation of compatibility is an equiv- 
alence relation for the set S. 


Ill. EXAMPLES 


As an example, consider the machine M, having three 
inputs and three outputs, which is defined by the fol- 
lowing connection matrix: 


alt S| beech 
(0, v) | 0 (c, B) 
(Gay. G0 (c, B) 


The state diagram is shown in Fig. 1. The absence of 
certain inputs in the rows of the connection matrix is 
assumed to result from restrictions on the inputs. 


(a,¢) 


Fig. 1—State diagram of a machine in which there exist input 
restrictions. 


For two states, 7 and 7, to be compatible it is necessary 
that identical outputs be associated with each input 
common to rows 7 and j of the connection matrix. 
Therefore we see that states 1 and 3 cannot be compati- 
ble since identical outputs are not associated with the 
common input symbol a. States 1 and 2 as well as states 
2 and 3 may be compatible. Examining states 2 and 3 
we observe that all input sequences in the intersection of 
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the two sets of associated allowable input sequences 
must start with the input symbol c. However, input ¢ 
takes both states 2 and 3 into state 3 and gives the out- > 
put 6. Consequently, states 2 and 3 are compatible. 
States 1 and 2, on the other hand, are compatible if and 
only if states 1 and 3 are compatible, since input 0 takes 
the machine from state 1 to state 3 and from state 2 to 
state 1. Since states 1 and 3 are not compatible, neither 
are states 1 and 2. If the states of the machine are 
partitioned into the two sets, S:= {1} and S,= {2, 3} 
as indicated by the dashed lines in the connection ma- 
trixof M, then the states in each set are pseudoequivalent 
after definition 3. Letting the set Me 3} be represented 
by a single state 2’ and the set {1} by 1’, we obtain the 
following connection matrix of a two-state machine M’ 
which is compatible to the original machine under the 
correspondence 1<>1’; 2, 3<02’: 


0 (a, a) V i 
(6, vy) V (a, ) (c, B) 


The corresponding state diagram is given in Fig. 2. 


(a,e)v(b,7) 


(a,7) v (b, 7) (c,6) 


Fig. 2—State diagram of a machine compatible to the one 
defined in Fig. 1. 


That the two machines are compatible can be ascer- 
tained in this example by comparing the input-output 
behavior of the two machines (a more formal treatment 
appears later). In other words, as long as only the 
originally allowable input sequences are presented to 
the two-state machine given in Fig. 2 it can be repre- 
sented as the three-state machine described in Fig. 1. 
The connection matrix of this latter machine is readily 
obtained by replacing each submatrix of the partitioned 


‘connection matrix of machine M by the union of all its 


entries, a rule which will be formalized in Section VI. 

It is important to note that the connection matrix for 
M’ no longer contains all the information concerning 
the input restrictions of M. This information is given 
only in the original state diagram and the associated 
connection matrix. That the machine described in Fig. 1 
could be represented by the compatible machine in Fig. 
2 depended upon the existence of a certain partitioning 
of the connection matrix as illustrated in the connection 
matrix M. A detailed development of this observation 
is given below. 

As a second example, consider the machine N, having 
4 inputs and 4 outputs, defined by the following connec- 
tion matrix, 
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One (Che IOI Siete (d;cy) 
eee Mey PricgaeneyG 
Cee ob 
| 
eameen | Goteiinn fetid: comenen) egueues  ecieceer commen bonds famnial | (3 Fe SF ES 
(Choo? a Osatiled aac a VS 
ee aa OO 


The associated state diagram appears in Fig. 3. From 
an inspection of the connection matrix we see that 
states 1, 2, and 3 form a set of pseudoequivalent states 
as does the set consisting of states 4 and 5. Indeed, the 
only common allowable sequence in the first case is 
ccc +++. In the second case there is none. Therefore, 
under the correspondence 1’<1, 2, 3; 2’<04, 5, a com- 
patible machine N’ which has only two states is ob- 
tained. The connection matrix of N’, which is given 
below, can be obtained as in the previous example. 


Leavian 0 


The state diagram is in Fig. 4. 


Fig. 3—Example of a machine in which there is more than one 
way to combine states under compatibility of states. 


(d,y) 


(c,o) 


(c,4)v(d, 7) 


Fig. 4—A two-state machine compatible to the one defined in 
Fig. 3. 


The above partitioning of states into pseudoequiva- 
lent sets is not the only one possible in this case. States 
2, 3, and 4 also form a set of pseudoequivalent states 
since there is no common allowable input sequence. 
States 1 and 5 are not compatible since the input c is 
associated with a different output in the two cases. We 
have then as another partitioning the sets S;= {1}, 
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S»= {2, 3,4}, Ss= {5}. (S, and S, cannot be combined 
since state 1 is not compatible to state 4.) Partitioning 
the connection matrix with respect to the sets Sj, S», 
and .S3 gives 


0 | (c, a) 0 0 | (d, y) 
Fh. Osea eee 
0 : (rayon 0 0 7 0 
(d, y) | 0 | 0 
pct wim opt la pe pra ee paises 
Om Wie OF Bic Meno anne 


The connection matrix of a compatible 3 state machine, 
N’’, under the correspondence 1/1; 2'<2, 3, 4; 3/065 
is given below while the state diagram appears in Fig. 5. 


0 (,«@) (, ) 
(d,y) (ca) 0 
OR CS) eee 


Fig. 5—A three-state machine compatible to the one defined in 
Fig. 3, but not compatible to the machine described in Fig. 4. 


A check of the connection matrix of N’’ reveals that 
no further nontrivial partitioning under pseudoequiv- 
alence exists. Since machines NV’ and N”’ do not have 
the same number of states we see that different simpli- 
fications of a given state diagram with respect to 
pseudoequivalence do not necessarily lead to represen- 
tations involving the same number of states. 

We also note that although machines N’ and W are 
compatible, as are machines N and N’’, machine N’ is 
not compatible to machine N’’. This illustrates the as- 
sertion that compatibility of machines is not necessarily 
transitive. 

To verify the noncompatibility of machines NV’ and 
N"’, observe that the input sequence ¢, d, ¢ presented to 
machine N’ when it is in state 1’ produces the output 
sequence a, y, B. If the same input sequence is applied 
to machine N’’ the above sequence of outputs is not 
obtained under any starting condition. 

The allowable input sequences of a machine MM are, 
as has been noted,® obtainable from the powers of the 
connection matrix C. The same holds true for the con- 
nection matrix C’ of an arbitrary reduced (with respect 
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to compatibility of states) form M’ of M. If a state s; 
of M maps into a state t; of M’, then the allowable input 
sequences of length 7 associated with s; all appear in 
row j of C’(r), along with any additional input sequences 
of length r which may be allowable for state t; of M’ by 
virtue of the combination of states of M. It is significant 
that allowable sequences of length r may well appear in 
association with states of M’ whereas these sequences 
are allowable for no state of the original machine. This 
re-emphasizes the fact that a reduced machine M’ is not 
ordinarily equivalent to the original machine M, but 
rather performs the same as M only for sequences allow- 
able for M. 


IV. FUNDAMENTAL THEOREMS 
Definition 
A matrix A whose elements are input-output poly- 


nomials® is called a generalized r matrix if it has the fol- 
lowing properties: 


1) All nonzero entries of A are homogeneous and of 
a common degree 7, 

2) All nonzero terms in each row of A have distinct 
input sequences, and, 

3) If an input sequence appears in more than one 
row of the matrix all associated output sequences 
are identical. (If all input sequences which appear 
in any given row also appear in every other row, 
then a generalized r matrix becomes an 7 matrix.) !° 


Lemma 1: The union® of two generalized r matrices 
A and B of the same order is again a generalized 7 matrix 
whenever the terms of corresponding rows of A and B 
contain no identical input sequences. 


Lemma 2: The product® AB of a generalized 7 matrix 
A and a generalized s matrix B is a generalized (r+s) 
matrix. 


Theorem 1: Let the submatrices of a symmetric par- 
titioning of a given connection matrix C be denoted by 
Ci, 4,7 =1, 2, --*,q. Let the submatrices of an identi- 
cal partitioning of C’, where 7 is any positive integer, be 
denoted by Ci;, 2, 7=1, 2,---+,9¢. If now each Cy isa 
generalized 1 matrix and if in each row of submatrices 
Ci +++ Cig no two distinct submatrices have any com- 
mon identical inputs, then each C;; is a generalized 
ry matrix and in each row of submatrices Cy + + + Cig 
no two distinct submatrices have any input sequences in 
common. 


Theorem 2: If a connection matrix admits a sym- 
metrical partitioning in which each submatrix is a gen- 
eralized 1 matrix and in which no two distinct sub- 


10 Thid., p. 280. 
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matrices in any given row of submatrices have an input 
symbol in common, then the states contained in each 
set of the corresponding partitioning of states are- 
pseudoequivalent. 

The proofs are similar to the proofs of the correspond- 
ing lemmas and theorems in the paper on equivalence 
of states. 


V. ALGORITHM FOR DETERMINING PARTITIONINGS 
INTO SETS OF PSEUDOEQUIVALENT STATES 


Consider the connection matrix C of any machine M. 
The following procedure will determine partitionings 
into sets of pseudoequivalent states. As we have already 
noted in Section III, the partitioning of the states of a 
machine into sets of pseudoequivalent states is not 
necessarily unique; therefore the number of sets in the 
resulting partitioning will, ordinarily, depend upon how 
the partitioning was obtained. 


1) Partition the states of M into disjunct sets Si, 
So, + + +, Sq such that the rows of C associated with the 
states in each S; form a generalized 1 matrix and such 
that this condition is not preserved if any two sets S;, 
S;, are combined. Note that this partitioning is not 
uniquely defined. Ordinarily, there will be more than 
one such partitioning. 

If the only partitioning is the trivial partitioning (into 
sets of order 1), then there are no compatible, and hence 
no nontrivial, sets of pseudoequivalent states in VM. 

2) If there is a nontrivial partitioning S,, S:, - +--+, S¢, 
reorder the states according to this partitioning, if neces- 
sary, so that the rows corresponding to states in the 
same S; are adjacent in the matrix, and hence it is possi- 
ble to partition the corresponding connection matrix C 
symmetrically with respect to the S; into generalized 
1 matrices. Indicate the partitioning by C= [C,;]. 

If the Ci; are all generalized 1 matrices so that for all 
jand k Ci; and Cx have no nonzero common inputs 
where C;; and Cy are any two submatrices in the ith 
row of submatrices, then the process terminates at this 
point because the states in each S; are pseudoequivalent 


_by Theorem 2. 


3) If the Ci; are not all generalized 1 matrices so that 
Ci; and Cx have no nonzero inputs in common, then 
partition the states associated with each row of sub- 
matrices.Cy, Ca, + 2 Cighe =o ie et ng intaneies 
junct sets Si, Se, +--+, Sa,;so that for a given S;; each 
corresponding set of rows in Cy, Cx, +++ , Cig is a gen- 
eralized 1 matrix; so that no two of these associated 
submatrices corresponding to a given S;; have any non- 
zero inputs in common; and finally so that these con- 
ditions are not preserved if any two sets Si;, Siz are com- 
bined. Again, there will be, in many cases, more than 
one such partitioning. If the only partitioning is the 
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trivial one for all rows of submatrices Ca, Cx, + + + , Cia, 
-4=1, 2,:--, g, and for all partitionings indicated in 
step 1, then there are no compatible states. 

_ 4) If there is a nontrivial partitioning, Su, Se, ++ -, 
Sinz, resulting from step 3, then repartition the connec- 
tion matrix with respect to the S;;, again reordering the 
states, if necessary, so as to have the rows of states in 
the same set S;; adjacent to each other in the matrix. 
Denote the new partitioning by 


Ce [C3] = [ D.;**]; Ci = [ Dij**] cj fixed) 


where the indexes k and ¢ are used to denote the sub- 
matrices in the partitioning of each submatrix C;;. 


If now the D,,*' are all generalized 1 matrices so that 
for all t and s D;;** and D;;** have no nonzero inputs in 
common, the process terminates since the states in each 
—S;;form a set of pseudoequivalent states after Theorem 
2. Otherwise, steps 3 and 4 are repeated. 

We have thus outlined a two-step process: first, par- 
tition the sets of states into disjunct sets so that the 
submatrices of the previous partitioning are generalized 
1 matrices with no two of the newly formed submatrices 
in the same row of submatrices having a nonzero input 
in common, and so that these conditions are not pre- 
served if any two of these newly formed sets of states 
are combined; second, repartition the connection ma- 
trix with respect to this new partitioning of states. These 
two steps are repeated until a partitioning of the con- 
nection matrix is obtained in which all submatrices are 
generalized 1 matrices with no two submatrices in the 
same row having an input in common. This is a finite 
process since the number of states is finite. When such 
a partitioning is obtained, then by Theorem 2 we see 
that all states in each subset of states in the final par- 
titioning form a set of pseudoequivalent states. If the 
subsets of states ultimately obtained are all of order 1 
for all possible partitionings indicated in the second step, 
then there are no compatible states, and thus no non- 
_ trivial sets of pseudoequivalent states. 

We have already noted that the number of sets of 
pseudoequivalent states ultimately obtained depends 
upon how the partitioning was effected at each stage. 
An unsolved problem of considerable interest is to de- 
termine the minimum number of sets of pseudoequiv- 
alent states without repeating the above process for each 
possibility and comparing results. 


VI. THe REDUCED MACHINE UNDER 
PSEUDOEQUIVALENCE 


After completing the process outlined in Section V, 
we define a new machine by replacing, in effect, the 
states of each pseudoequivalent set by a single state 
representing that set. To account for the transitions of 
all the states of one set to states of other sets under the 
various inputs to the machine M, we replace each sub- 
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matric of the final partitioned form C by the union of 

all its entries. For simplicity of notation denote the final 

partitioned form of C by 
C= [Ci], i,j = 1) 2; Ma Ma A 

and the corresponding sets of pseudoequivalent states 

by Si, S:,+-+, Sy. Then if dj; is the union of all the 

entries of Ci;, the gXq matrix 


D = [d3] 


is the connection matrix of a machine NV with gq states, 
which is compatible to M. That D is indeed a connec- 
tion matrix is seen by noting that because C is a con- 
nection matrix, the C;;’s are generalized 1 matrices, and 
no two matrices C;; and Cy in the same row of sub- 
matrices contain any common nonzero input symbols. 
Therefore, d;; and dx will contain no common nonzero 
input symbol, which is the only condition to be 
satisfied. 

The machine NV defined by D is called a reduced form 
of M under pseudoequivalence. (It may very well be 
the case that there exist compatible states in V and even 
sets of pseudoequivalent states.) It remains to show 
that M and WN are compatible machines. To prove this, 
for each positive integer 7 partition C" according to the 
indicated sets of pseudoequivalent states of M, and let 
C;; denote any of the resulting submatrices. The en- 
tries of C;; define the allowable input sequences of 
length r which take the states of set S; into those of set 
S; and also give the associated output sequences. 

Now the entry d;; of D is the union of the entries C;; 
and hence defines the set of all inputs which take some 
member of S; into some member of S;. On the other 
hand, D;;, the 7j entry of D’ , is given by 


Tr 


dj = We 


ky,kg,+++,ky—y=1 


Cit Cirke * 


oe deed 


Since each state of M appears in some S;, every term 
appearing in C,; has its factors appearing in some 
product of the form dix,disng * * + Cx,_,3, and hence is a 
term of d,;. 

From these observations we see that all the allowable 
input sequences associated with each of the states of 
M in S; are also allowable input sequences for state 2 of 
N (the converse is not necessarily the case), and the 
output sequences associated with a common input se- 
quence are identical. Hence, by Definition 2 of Section 
II, the machines M and WN are compatible. We then 
have the following theorem. 


Theorem 3: Given a machine M, there exists a corre- 
sponding compatible machine N (not necessarily 
unique) with minimum numbers of states. 

At this point, the minimum number of states can be 
determined only by examining all possibilities. 
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VII. AN EXAMPLE 


In order to illustrate the algorithm in Section IV we consider a machine having 6 inputs—a, 8, c, d, e, f; 4 outputs 
—a, B, y, 5; and 13 states, which is defined by the following connection matrix: 


1 2 3 4 5 6 Y) 8 9 10 11 iD 13 
1 0 (e, 8) 0 0 0 0 0 0): =(ata) Vile, py 2E0 0 0 0 
2 0 (a, B) (e,8) 0 0 0 0 0 (c, B) (d,a) (by) 0 0 
3 (e, 5) 0 (b,7) 0 0 0 0 0 (a, «) Oavir(ereyrn 20 0 
4 (a, a) 0 0 (7 0 0 0 0 (e, 8) 07" = (e620 0 
5 0 0 (do) 8) “Oy (0° VG By 0 0 0 0 0-9, eSB) 
6 (e, a) 0 Gay 70 O@ 2G; 6) nO 0 0 (6, 8) 0 (d,a) 0 
71 (c,B) V (e, 8) (b, y) 0 0 0 0 0 0 0 0 0 0 
8 0 0 0 0 0 0.0 (Eire) a0 0 04 sb 0 0 
9 0 0 (b,7) 0 0 0 0 (e, 8) 0 0 0 0 0 
10 0 0 0 0 0 0 0 (0,7) (c, B) (a,a) 0 0 0 
11 (c, B) (b,y) V (4,8) 0 0 0 0 2¢@5e)) 4.0 0 0 O.- —Wechye ane 
12 0 (e, 8) 0 “Gee 20 0 0 0 0 0 0 (ar a0 
13 0 0 Gy) 0 0 O= 16,.6).2 20 (d, «) 0 0 (ea) (a,8) 


There are two distinct, disjunct partitionings so that the corresponding sets of rows form generalized 1 matrices 
and so that this condition is not preserved if any two sets of states are combined: P,: {1, 3; 4 LOnceo. sie gt, 
{2, 5, 11, 8}, {6, 13}; Po: {1, 3, 4, 7, 9, 12, 10, 8}, {2, 5, 11}, {6, 13}. If the above connection matrix is parti- 
tioned with respect to P;, we have 


1 BWihaial a7 9 10) Sf) 8 2 soa kt eens 
1 0 Oe Ole Oa (eh ee O, fo ak0 al peo (e, 8) rey Oma geo 

3 (e; 8) (Coane (ie onhe (4, a) Oude Wie 0 OL n(ep ea One 

4 (4, a) ONG) 12.0 (e, 8) 9: MOR 0 Gee Oe 

Tle CeBy see 0s. On BO 0 OO a (, ») ree 6. oe 

9 0 Gao ane 0 OL atom Wiel 0 Guditard Ogun 
10 0 hee es (Gia)s 6 elaca) Moles) 0 fea ge eg ieas 
12 0 Cees ye n0 0 He et (e, 3) 080 Ooi an 

8 Oem eta - Un aaa) 0 on ee 0 Ov) (Bi) wie OS a One 
2 0 (e, B) 0 0 (c, B) (d,a) 0 : 0 (a, 8) 0 , 7) 0 0 

5 0 (d,a) (eB) (a, B) 0 Oy 540 : 0 0 (b, y) Ores a reO (c, B) 
11 (c, B) 0 0 (,a) 0 0 (e, B) : 0 (7) V (4,8) 0 0 | 0 0 

6 (Bei ACRE ar OE 0 6) eo QT AG aoe GT Oa ats nl 
13 0 (CA) 0 (0, 6) (d, a) Oe CAR | 0 0 0 0 0 (@,8) 


We note that the rows of submatrices corresponding to the two sets, ae ee ew me ee 12} and {8, Dies 11} are 
not all generalized 1 matrices so that any two submatrices in the same row have no common nonzero input sym- 
bols. Therefore, the states in these two sets must be partitioned in accordance with the algorithm. The only two 
distinct, nonsuperfluous disjunct partitionings are: 


Put (1, 9, 12}, 13, 4}, 17, 10}, {2, 12, 8, (Si (6, 13) Pee | t, 9, 12) 4347, on aloe qo 
{6, 13}. The partitioning P»; leads to the following connection matrix: 
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1 9 12 3 4 7 10 2 8 11 5 6 136 = 
1 0 CEVA Cea 0 0 0 (e, 8) 0 0 0-7 =0 0 
9 0 0 ele (b,7) 0 0 0 | 0 (3)a 26 a | 0 0 
12 0 0 (ala) 46 (b, 7) a 0 (e, 8) 0 pine 10 | 0 0 
gl oe OO cin no Siete att ear Gariee eee 
4 (ane) (e, 8) eats (b, y) esa 0 0 0° ra) = <0 0 0 
GAVE Den ETE 
10 0 (c, 8) pie eer 0 | 0 (a, a) 0 (,7) 0 6 0 0 
2/9 ae iol salt plain ae Cai ae lar i 3 a 6 EWE May aoa 
8 0 0 0 0 0 (d, a) 0 0 0 (6, y) 0 | 0 0 
11] &,8) 0 RoWes 0 Nasr TOTEM) Ulhok WB As OT eqokai Osea wCbebeD 
2 a aa ee cae marek eae eee 
MM ee LO Ga) Ga) Own 0 (G8 \ 610 Gol us OROn 
13 0 (d, a) (e,a) * Uy) —-” (; 8). 20 : 0 0 0 | 0 | 0 (4,8) 


The row of submatrices corresponding to set of states {2, 8, 11} are not generalized 1 matrices so that no two sub- 
matrices have a common nonzero input symbol. There are again two possible partitionings: Po: { RES 17 } : {3, 4}, 
fee 10) 2,8), 111}, 16n13}5 Poe: (1) 9,12}, (3, 4}, {7,10}, {2} {8,41}, {5}, 46, 13} The partitioning Pa, 


results in the following connection matrix: 


1 9 12 3 4 - 1p 2 8 ch es eg eae 
{or 0 exile 8) yOu 0 0 0 Oink fj ween) 0 0 Ory PRO 0 
0 0 ” piso rath oe 4) 2-H Gin. to PGI. 
12 0 0 en eas pe ia visae0 «0. Kigges) Oo SAO, EO Te SORE 
2 hcp Pao 5 eel aes ‘ws CREY Vee Olea Pe AMET ORI, ees 0 Siete 0 oe ene 
4 (a, «) (e, 8) REE hs aio 0 WEL Rear nme 
oo ; — . ea eee ee : ae 
10 0 (c, B) Stilo MlOpsi=a0: "Gra ee eo (Sy Beare Oey Re aay 
an PAE ACS reall tele oe 0 ee eo 3 0 Gap Xap). — 0 6,0 0 | 10 att “0 a 
8 0 0 0 Oia datk Oniocufdsc) Remo 0 O sextby\ 3s 0, Caen 
ee ee i Gi EG es 
Bk a a Oe) OP ED soitinenbert on Or iomi anobna erp Dita aeglle 
6 (e, «) 0 CO ant! IOS Dea ean ble | 
13 0 (d, a) (ey CPG ey) 1) OSH =(05. 8) LOH yy 0 0 0 Ou eve0 erpase) 


We observe that all submatrices are now generalized 1 matrices with no two submatrices in the same row having 
nonzero input symbols in common. Thus, each set of states in the corresponding partitioning of the states forms a 
set of pseudoequivalent states. The reduced form of the machine relative to this particular partitioning into pseudo- 


equivalent states is given by the following connection matrix: 
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1 3 
1 F(a, a) V (¢, B) (, ) 0 
3 | (a, a) V (e, 8) (, ) 
dal (err B) N/E n.0) 0 (a, a) 
2 (c, B) (e, B) (d, a) 
11 | (c, 8) V (e, 8) 0 (d, a) 
5) 0 (d,a) V (e,8) (a, B) 
6 L(d, a) V (e, a) (fv) (, 9) 


There are no compatible states in the above connection 
matrix. 

We do not include a corresponding development for 
the other possible partitionings. P22, for example, leads 
to the trivial partitioning. P, gives rise to a number of 
different decompositions into sets of pseudoequivalent 
states, none, leading however, to a compatible machine 
with fewer states than the machine resulting from Pou. 


VIII. ConcLusION 


The relation of compatibility of states and of ma- 


Theoretical Consideration 
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2 11 5 6 

(e, 6) 0 0 0 

0 (cB) 0 0 

(6, Y) 0 0) 0 

(a, 6) Chey de 18 0 

(6,v) V (2,8) 0 0 0 
0 O On¥y CG a8) 
0 0 0 (a, B) 


chines permits a simplification of state diagrams in addi- 
tion to that provided by the relation of equivalence of 
states and of machines. A method has been developed 
for reducing a given state diagram to its “simplest” 
compatible form. 
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of Computing Errors of a 


Slow Type Electronic Analog Computer* 


T. MIURA anp M. NAGATA 


Summary—When solving differential equations using analog 
computers, there are two causes of errors: one results from the 
integrators and the other results from the coefficient setting ele- 
ments. The former is independent of the order of the differential 
equation if the integrating time constants are all equal, while the 
latter is dependent on the setup procedure. The amount of the error 
from the coefficient setting elements is comparatively small in the 
operating frequency range from w<rad/sec. Therefore, theoretical 
analysis considering errors due to the integrators only is quite suffi- 
cient. However, higher accuracy is required in w= 10 to 100 rad/sec 
(for example in flight simulators) and consequently the error caused 
by the coefficient setting elements has become impossible to neglect. 

The authors have succeeded in formulating a generalized and 
practical equation developed from conventional formulas. 


INTRODUCTION 


HERE have been several reports, including 
Macnee’s, of analyses of errors caused by solving 
linear differential equations by means of elec- 
tronic analog computers. However, from the standpoint 


* Manuscript received by the PGEC, June 25, 1958. 
{ Hitachi Central Res. Lab., Tokyo, Japan. 


of practical usage, the analyzed computer setups in all 
these reports are not general; for instance, Macnee! 
treated the case in which all the constant coefficient 
multipliers are used to set coefficients without consider- 


ing the sign of the operational amplifiers and in which 


the constant coefficient multipliers are fed back through 
the adder to the integrator corresponding to the highest 
order of the differential equation. Marsocci? treated the 
case where the potentiometers and the sign changers are 
used as coefficient setting elements and the sign changers 
are used in order to align the sign. A further condition 
was that every coefficient k» was in the range 0<ky <1. 
Kurokawa?’ treated the computer setup excluding the 
constant coefficient multiplier used in Macnee’s case, 


1A. B. Macnee, “Some limitation on the accuracy of electronic 
differential analyzers,” Proc. IRE, vol. 40, pp. 303-308; March, 1952. 

?'V. A. Marsocci, “An error analysis of electronic analog com- 
puters,” IRE Trans. oN ELECTRONIC CompuTERs, vol. 5, pp. 207- 
212; April, 1956. 

*K. Kurokawa, “An Error Evaluation of Electronic Analogue 


Computers,” preparatory essay of the seminar of Japan Sci., Assn.; 
1957. (In Japanese.) 
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while Nomura‘ pointed out that a part of Macnee’s 
analysis was incorrect and showed the required cor- 
rection. 

In all the above cases except Nomura’s, the frequency 
characteristics are given in the form of the operational 
amplifier (closed-loop characteristics). The authors®-7 
already have introduced a practical and general equa- 

_tion of the frequency characteristics in the form of the 
| dc amplifier (open-loop frequency characteristics) on the 
assumption that errors caused by the constant coef- 
ficient multipliers and the adders may be assumed to be 
: negligibly small compared with the errors of the inte- 
| grators.® 

_ However, recently higher accuracy in the operating 
frequency range of w>10 rad/sec (such as in flight 
simulators) has been required, and it is very often neces- 
sary to connect constant coefficient multipliers in tan- 
dem. Consequently it is impossible to neglect the errors 
due to the constant coefficient multipliers and the errors 
due to the associated stray capacities. Therefore it is 
essential to develop general formulas considering the 
errors of the coefficient setting elements. 

This paper discusses theoretically the errors due to 
computer setups, including constant coefficient multi- 
pliers and potentiometers, efc. In the following discus- 
sions, computing impedances are assumed to have proper 
values. 


ANALYZED COMPUTER SETUP 


The circuit shown in Fig. 1 is a computer setup for 
solving a linear differential equation with constant co- 
efficients, and it is used in this analysis. 

The newly-developed formula has the following fea- 
tures: the general expression of the transfer function 
Km(p) for all coefficient setting elements is applied for 
all the cases of potentiometer or constant coefficient 
multiplier settings and the sign of the integrator is in- 
cluded. This also is true for all the other computing ele- 

-ments, and errors due to the adders are included in the 
sign changer as a first order approximation. The dif- 
ferential equation which is solved by the computer 
shown in Fig. 1 is 


pap T”p"kmy = 0 where k, = 1. (1) 


m=0 


Now, considering the transfer function I(p) of the 


4T. Nomura, “An Error of Electronic Analog Computer,” In- 
dustrial Lab., Tokyo University, Tokyo, Japan, Rep. No. 414; 
1953. (In Japanese.) : ; 

5 JT. Miura, “The electronic analog computer,” Hitachi Rev., 
(special issue for electronic tube and application of electronic tube) 
special issue No. 3, pp. 139-145; 1953. (In Japanese.) 

6 T. Miura, M. Aoki, M. Hibi, and T. Numakura, “On the pre- 
cision of a slow type analog computer,” J. IEE Japan, vol. 76, pp. 
815, 896-905; August, 1956. (In Japanese.) 

7T. Mirua, Z. Abe and M. Nagata, “On Computing Errors of a 
Slow Type Analog Computer,” essays presented at IEE Conference, 
Tokyo, Japan; November, 1957. (In Japanese.) ; 

8 In practice it is possible to neglect the errors due to the inte- 
grators when operating in the frequency range of w<10 rad/sec. 
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Fig. 1—Computer setup for solving equation (1). 


Integrating capacitor 


Integrating 
resistor 


: direct coupled amplifier 


chopper omplifier 


Fig. 2—Equivalent circuit of integrator. 


integrator and the transfer function Km(p) of the co- 
efficient setting element in Fig. 1, we have 


Ip)? = {1 + (1 + wd Z.(p) V(p)} Jus, (2) 


where Z;(p), Y(p) are the equivalent input computing 
impedance and the feedback computing admittance re- 
spectively. (In the following considerations we neglect 
the effect of the grid current and the drift voltage of 
the integrator, which will be treated in Appendix I.) 
The following relationships are obtained if we take into 
consideration the inductance and the stray capacitance 
across the integrating resistor, and the leakage resist- 
ance and the dielectric absorption across the inte- 
grating capacitor in Fig. 2. (For the dielectric absorption 
phenomenon, Cole-Cole’s experimental formula® which 
gives a good approximation for high molecule materials 
was applied.) 


Zip) ~ R{1 — C.Rp + (L/R)p} 


V(p) c( ran leaas ~) (3) 
aE. Clarita eGR s 


From (2) and (3) 


1p = T |— 24 ! Hee ah | 


pil Mi 
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where 
AG OR ip 2 
A ae aed Ala (e) Bere Wat 


The general formula for Km(f) is 
Kn(p) = bm{ 1 — (Kmn/tx) + ptm}. (5) 


The second and third terms of the right-hand side of 
(5) show the error terms caused by amplifier gain uw, and 
the stray capacitance across the computing resistors re- 
spectively. The determination of Km, 7m, can be seen by 
referring to Table I. 

If we used the transfer function for each of the com- 
puting elements, the characteristic equation for the case 
of computer setup of Fig. 1 is 


Do 1(p)-"Km(p) = F(p) = 0. (6) 
m=0 
Assuming that the characteristics of the operational 
amplifier are ideal, then from (4) and (5) we have 


I(p)*=T-p, Km(p) = Rm. (7) 


Substituting (7) into (6), we find (6) is identical to (1). 
We will consider the significance of (6) in the following 
section. 


INTRODUCTION OF THE GENERAL FORMULA FOR 
OVER-ALL COMPUTING ERRORS 


It may be adequate to think of the quantity ¢ for ex- 
amining the error which is given by the difference be- 
tween one root p,' obtained from the characteristic (6) 
and the characteristic root pa from (1). By applying 
Taylor’s expansion theorem, (6) becomes 


dF(p))  & P @F(p) 
IMA) = 1G € a 
(p.!) (oo + | estes | + 


S <|s FP al fe 
miL dp” In, 

In general d”F(p) /dp” has a finite value and the above 
series converges for any value of e. Furthermore, com- 
paring the second and the third terms, obviously the 
latter is smaller than the former at the computing fre- 
quencies in which case we may assume 1>>|e/p| to be 


correct. Therefore we can neglect the third and the fol- 
lowing higher order terms. 


e= — F(p.)/|dF(p)/dp|p,°. (8) 
Calculating dF(p)/dp from (6) 
dF (p) =e 
SOF ie DL mrp) al) | 
Sap fs : 
+p oe 5 I(p)- (9) 


The term —I(p)-?-dI(p)/dp, dKm(p)/dp within (9) can 
be reduced to 


9 We write pa=p hereafter unless otherwise stated. 
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TABLE I 
Tue Metuop oF DETERMINATION OF Km(p) AND Tn 


(a) Determination of Km(p) 


1) Coefficient setting po- | Km=0 
tentiometer. 

2) The tandem connection 2 
of N constant coefficient = x (1+Cm)* 


multipliers. 
where Cm is the multiplication of each 
constant coefficient multiplier; 7.e., 
for sign changer, Km=1+1=2. For 
the tandem connection of constant co- 
efficient multipliers with multiplica- 
tion 5 and 6, Km=6+7=13 


3) (Coefficient setting po- | The same as 2). 
tentiometer) + (constant 


coefficient multiplier). 


The determination of Km is possible 
from 2) on the assumption that pass- 
ing through one sign changer is equal 
to passing through one adder. 


4) The case through the 
adder. 


(b) Determination of tT» 


1) Constant coefficient 
multiplier. 


Tm = Tim — Thm 
siti =Tim 

CapRy =Thm 

where 

R;: input computing resistor 

R;: feedback computing resistor 

C.i: stray capacitor across the input 
computing resistor 

C.z: stray capacitor across the feed- 
back computing resistor. 


SS Tm 


2) Tandem connection of 
constant coefficient mul- 
tipliers. 


3 ™ =11—T» 
Ti=(1—km) RC) 
=km(1—km)R(Ci + C2) 

Ci, O~100~300PF 

where 

R: resistance between A and B ter- 
minals 

Ci: stray capacitor between A and C 
terminals 

C2: stray capacitor between C and B 
terminals 

km: dividing ratio 

A and B are each terminals of the co- 

efficient setting potentiometer and C 

is the slider. 


ey, 


Coefficient setting po- 
tentiometer. 


* See Appendix IT. 
, al (p) 
root AG) Ra ae 
dK» () 
dp 
Substituting (10) which has been obtained by the first 
approximation of (9) on 


dF ( re z 


~ T{1 + (1/us) + (dA(p)/dp)} 


Li Tm’ Rae 


(10) 


Tn 1 + — + —+ ay 
dp 


Mi 
A 1 

+m o(2 + Py — frat pm 
Mi p MEK 


al n 
+ — >) m(m — 1) Tp" 2m 
Ki m=1 


n=l 
+> Trakigli po. 


m=0 


(11) 
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In case the first approximation is possible for e, dF(p)/dp 
can be approximated by the following: 


AON ae a 
i ~ > mT"Rnp™. 


m=1 
Next putting (4) and (5) into (6), then 


(12) 


n 


A Ky, 
Pe) = refi em Ph (1 A8 + pra) hs 
m=0 Mi p UK 
at bea a m= 1 
bee mInps + 
Kil m=1 M1 
A Ron 
+ (m — 1) ak an (1 eet pra) 
p MK 
a A Ky 
~ LI pin Su sea 5 bral 
m=0 Mi p MK 
1 n 
+ — )imT"p"™ Rm. (13) 
Mil m=1 
Substituting (1) and (12) into (13) 
% m NA) an Som 
F(p) = 2 Tyre | ee + Prat 
m=0 Mi p UK 
len ae 
(6) ie 
Mil dp 


From (8), (13) and (14) 


>> T"p"kimKm 


1 al tall cones 
Saal) aval 
Bs ee me SS mT™p" Rm 
m=1 


n 


ye Tp" mT m 


= na? = 
>> mT" ph" Rin 
m1 
Rwhere 
ap) = — CRe + C/e +S 2 
pases SE COM Cee CR, 


Fae Om K=O 115) 


Eq. (15) is the generalized equation of errors. The first 
and third terms of the right-hand side of (15) are the 
error terms caused by the dc amplifiers of the integrator, 
while the second term shows the errors caused by the 
computing impedance of the integrator, and the fourth 
and fifth terms show the errors caused by the de ampli- 
fiers of the constant coefficient multiplier and the com- 
puting resistor respectively. 

From (15) it is obvious that the error of the integrator 
is quite independent of the order of the differential 
equation in the case where the integrating time con- 
stants are all equal. Therefore, the estimation of the er- 
ror of the integrator by the circle test (analysis of the 
second-order differential equation with no damping) is 
possible and, what is more, it may become the necessary 
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condition of the determination of the over-all accuracy 
when using this test. 

The approximating equation for ¢, taking the effect 
of the frequency characteristics of the amplifiers into 
consideration, can be easily obtained from (15) by treat- 
ing w; and pw, as functions of p and considering the effects 
of their frequency characteristics. A combination type 
dc amplifier which has an automatic balancing circuit 
is used. Assuming that time constants, which provide 
the gain and the frequency characteristics of the chopper 
part, and the direct-coupled amplifier part are as shown 
in Fig. 3, then 


sy Ho(1 + pro) 
(Lat pray (lh 4 pra) 


where po =p1' Ms, To=Ti/pa, Mi> 1. 
In (16) assuming p=,’ and substituting this equa- 
tion into (15), we can obtain by the method of succes- 


(16) 


be 


GAIN (db) 


Ve, loga@ 


y To 


| 
% 
Fig. 3—Frequency characteristic of the operational amplifier. 


sive approximations the following equation as a first 
approximation. 


~—_ Otro ( 3) 4 
ESS Gu Gyo 


Tt pek hp 
(1 + prei)(1 + prez) x P 


box(1 + prxo) f 


Dy mT™p” Rm 
m=1 
SS. De em 
m=0 
Ss (17) 
S Mp Rin 
m=1 


Eq. (17) is the approximating equation for ¢ taking into 
consideration the effect of the frequency characteristics 
of the amplifiers. 


EXAMPLES OF APPLICATIONS OF THE GENERAL 
FORMULA FOR € 
Circle test: Example 1 


The computer setup in this case is shown in Fig. 4. 
The equation to be solved by it is 


310 


{ 
K Ty 


Fig. 4—Computer setup to solve equation (18). 
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Fig. 5—Computer setup to solve the equation 


M 
prTNy + TJ Kny = 0. 


TABLE II m=1 
JA 
Constants OF Km(p) FOR USE WITH FIG TABLE III 
m km Km ConsTANTs OF Km(p) FOR USE wiTH Fic. 5 
2 1 0 K 
A 0 0 m km m 
0 K 1+K N 1 0 
N-1 0 0 
N-2 0 0 
; -- 0 0 
ay Vat pic 0. (18) 1 0 0 
M M 
The result of finding Km by means of Table I is shown in 0 I Km » (1+Cm) 


Table II. Putting the result from Table II into (15), the 
error C(p) of the coefficient setting elements is 


G = = 
(?) 2T2p = 2T? 
Koi Keep. OK ee 
Se en) ped. eT? 


Substituting the value p=j./K/T obtained from (18) 
into the above equation, we have 


1+K 1 
C(p) = — Tm 
(p) te uae 
p + K 1 
Se Ap) — — ptm. (19 
€ Apo (p) ii ee eee (19) 


The result of (19) is quite the same result that was 
obtained already.’ 


Tandem connection of N-integrators and M-constant co- 
efficient multipliers: Example 2 (Fig. 5) 


In this case the values of Km(p) are given by Table 


III. Therefore, from (15) 


M M M M 
De(t+ Cm) TL bm = Le tm: LD bn p 
m=1 


1 m==1 m=1 m= 
C(p) = SSS ——. (20) 
LK LVI PY 1 NTN pN 1 
In considering the relation 
M 
pale so ae ul Rn 
m=1 
1 p M M 
N HK m=1 m=1 


We can easily understand from (21) that the errors of 
the coefficient setting elements in this case are propor- 
tional to the reciprocal of the number of integrators and 
increase as the algebraic sum of errors of each of the 


constant coefficient multipliers. This result is very useful 
in estimating the computing errors of each loop in case a 
complicated circuit setup is being treated.1° 


Computing errors in the case of 0<Rm<1: Example 3 
In this case we can use only sign changers and co- 
efficient setting potentiometers as coefficient setting 
elements. This circuit has been analyzed by Marsocci 
(Fig. 6). In this case the values of Km(p) are shown in 
Table IV; therefore in (15) we put tTm0 and 
n—1 
2 >> bnkmpb™ 


dm—0 1f mis odd. 


C(p) = an , where (22) 
a D> mT™p™ bm Sm—1 if m is even. 
m=1 
Example 4 


Let us consider the computer setup shown in Fig. 7. 
This circuit was built up so as to obtain the same result 
as in the case of the circle test without any kind of minor 
loop by giving positive damping by loop K’ and nega- 
tive damping by loop K’’. If the characteristics of the 
sign changer are ideal, and the input characteristics of 
the two summing integrators are the same, we can 
examine the magnitude of the errors of the minor loop 
by comparing the setup having loops K’, K’’, with the 
one having no minor loop. We then calculate the incre- 
ment of the theoretical value of the errors due to the 
minor loop. In this case the constants of Km(p) are 
shown in Table V by virtue of Fig. 7. Substituting these 


values into (15) (in this case we may put we Tm&0), 
obtain 


_ 7° We can separate all kinds of circuits into several separate 
integrators and coefficient setting amplifiers, no matter how compli- 
cated the original circuits may be. 


Fig. 6—Computer setup to solve the equation 


> Tp kny = 0 


m=0 


where 0<km <1. 


Fig. 7—Computer setup to solve the equation 
pT’y + Ky = 0. 


TABLE IV 
CONSTANTS OF Km(p) FOR USE witH Fic. 6 


nt km Km 
nN 1 0 
n—1 Res 0 
i= 2 R23 2 
a— 3 Rees 0 
n—A4 Res 2 
2 ke 2 
1 ke 2 
1 ky 0 
0 ko p 


—2Tp+ K+ kK) 1 
2T*p MK 
1 Kii+ K) 
7% p. 
2uK 


C(P) 


(23) 


uk 
The first term on the right-hand side of (23) is due 


to the increment of the error from the minor loop. 


Example 5 


Next we consider a linear differential equation having 
a continuous oscillatory solution of order three. 


p Ty + 3p°T’y + 3pTy + 9y=0 p= +jV3/T. (24) 


We may think of Fig. 8(a) as a computer setup for this 
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Dip te 
V 


(a) 


pai Bane, 
FEY 
M7 


Fig. 8—(a) Computer setup to solve the equation 
pT*y + 3p?T2y + 3pTy + 9y = 0. 
(b) Computer setup to solve the equation 
pT*y + 3p?T*y + 3pTy + 9y = 0. 


TABLE V 
ConsTANTS OF Km(p) FoR USE witTH Fic. 7. 


m km Km 
2 1 0 
pfx’ path 1 0 
\K"" path —4 & 
0 ; JK 14K 


equation. The table for Km(p) is shown in Table VI re- 
ferring to this figure. 
Therefore from (15) and Table VI, we have 


1 18)°7* = 12p7. + 108 


CP) = 
uw 3p?T? + 6pT? + 3T 
Putting p? = —3/T? into above equation, we obtain 
Cen eaee, Ganai (25) 
Ao 4 uT 


Fig. 8(b) may be considered as a setup circuit for (24) 
and we can replace two sign changers a and 6 with a 
single adder. The constant, Km(p), in this case may 
also be obtained from Table V referring to the standard 
of Table I. 


CONCLUSIONS 


It is clear from the above examples that the new 
general formula is closely connected with practical com- 
puter setups and that the errors are easily calculated 
when the circuit is given. Moreover, the errors which 
are calculated by the conventional formula are included 
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TABLE VI 
CONSTANTS OF Km(p) FOR USE wiTH Fic. 8 


km Km 


CORN WwW 3 
OW We 
NLAO 


as a special case of the general formula which has been 
developed here. Several practical examples are shown 
as examples of the application of this formula. The 
authors believe that the general formula for the errors 
from the tandem connection of N integrators provides 
some estimate of the magnitude of the errors which arise 
in the case of complicated operational circuits. 


APPENDIX | 


EFFECTS OF THE DC AMPLIFIER GRID CURRENT 
AND DRIFT VOLTAGE ON € 


In the previous analysis of €, we assumed that the 
grid current 2, and the drift voltage eg of the dc amplifier 
may be neglected. Now if we consider the effect of 7, and 
ea, the following relation can be obtained between ym 
and Y,,-1 of the output of the mth and the (m—1)th 
integrator from the right-hand side in Fig. 1. 


Ym—1 = I(p)(y + ea — tgR)n — Cams (26) 


In (26), the subscript m of the bracket shows that the 
subscripts of the quantities in the bracket are all m. 
The computer setup shown in Fig. 1 is for the solution 
of the equation. 


» RmYm = 0. 


m=0 


(27) 


From (26) and (27) the following equation can be de- 
rived by putting Y,=y 


ym = Up)-my +  1(p)™(ee — igR)s 


i=m+1 


aay Da Lp)? COD ea: 


i=m+1 


(28) 


By substituting (28) into (27), one obtains 


n 


Do Rm = Ip)” D2 Rmyl(p)—™ 
m=0 0 


m= 


+> DY ballQ) (en — i,R)s 


m=0 i=m-+1 


= Dy oy tabi) Pers = 0, 


m=0 i=m+1 
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By dividing both sides of the above equations by I(p)*, 
one obtains 


n 


hndinl (6) + r(p)-"4 Se eee a nl (py 


m= 7=0 


=D (ex) D bar(gymt = 0. 2) 


m=0 


Since I(p)— is a differential operator, then by assuming 
that eg and 7, are constant and independent of time, the 
second term has only (eg—7,R), left and the third term 
vanishes. Eq. (29) becomes 


Sy RmI(p)—"y = — (€a — tgR)n = const. (30) 


m=—0 


For the ideal case I(p)-!=p 


> Rmp™y = const. 
m=0 
s 


Obviously the result of (30) shows that the grid current 
and the drift voltage have effects only on the initial 
condition voltage but not one. 


APPENDIX II 
CALCULATION OF Km(p) 


When the tandem connection of two constant co- 
efficient multipliers with multiplication constants C, 
and C; is used, the transfer function Km(p) is 


a a 
US eae ee RS er EN 


1+ C 1+C 
~cc,(1- == ‘\(a- te ‘ 
u bu 


i 


Kn(p) = “C1iC2 


Kee : 
= kn (1 _ ), Kn = >. Ge Ga) 
MK m=1 
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BIDEC—A Binary-to-Decimal or Decimal-to- 


Binary Converter* 
JOHN F. COULEURt 


Summary—Simple, high-speed devices to convert binary, binary 
coded octal, or Gray code numbers to binary coded decimal numbers 
or vice versa is described. Circuitry required is four shift register 
Stages per decimal digit plus one 30-diode network per decimal digit. 
In simple form the conversion requires two operations per binary 
bit but is theoretically capable of working at one operation per bit. 


INTRODUCTION 
Bes following article is a description of devices to 


convert between the binary and decimal form of 
numbers. The method of conversion imposes no 
limitation on the number of digits and, in proposed ad- 
vanced form, requires a time equal to the clock period 
times the number of digits to perform the conversion. 
BIDEC is the result of a need for a fast; simple, 
binary-to-decimal converter to be used as the output 
buffer of a scale factor computer to reduce missile flight 
data. A study of conventional techniques indicated that 
time available for the conversion ruled out the use of 
counters. Complexity, size, and cost made the conven- 
tional division by 10 undesirable, and word length 
eliminated matrix conversion. Several other special 
methods of converting binary to decimal and decimal to 
binary were found,! but they either had inherent speed 
limitations or became increasingly complex with the 
number of digits to be converted. 
The resulting search for an entirely new approach to 
the conversion of both binary to decimal and decimal 
to binary led to the formulation of the principles of 
operation of BIDEC, the original embodiments of which 
will be treated in this article. In the simplest form, 
either a binary or reflected binary number of any length 
-in serial or parallel form may be converted to a binary 
coded decimal number, in a length of time equal to 2X 
number of binary digits X shifting time. Circuit re- 
quirements consist of four shift register stages per 
decimal digit, one, 30-diode logic per decimal digit, and 
a clock to count shift pulses. A change in wiring enables 
the conversion of decimal to binary, and an increase in 
complexity permits either conversion in the same device. 


THEORY 


Making use of the fact that the decimal value of any 
digit in a binary number is 2"—!, where J is location of 
the digit in the number, the decimal value of a binary 
number may be determined by adding appropriate 
powers of 2 as indicated by the presence of 1’s in the 
binary number. The powers of 2 may be readily calcu- 

* Manuscript received by the PGEC, July 22, 1958. Revised 
manuscript received, August 7, 1958. 


+ General Electric Co., Syracuse, N. Y. 
1 Richards, Arithmetic Operations in Digital Computers, et al. 


lated by doubling one N—1 times or by multiplying the 
1 in the binary number by two N—1 times. 

A shift register is a convenient method of multiplica- 
tion of a binary number by 2 since shifting the entire 
register once towards the most significant end doubles 
the value of all binary digits in the register if the stages 
are assigned arbitrary values of multiples of 2. With as- 
signment of values to the stages of the shift register 
being somewhat arbitrary, there is no reason why the 
shift register may not be broken up into groups of four 
stages to be called decades. The four stages within the 
decades are assigned values of 1, 2, 4, and 8 ascending 
in the direction of shift, but the decades are assigned 
decimal values of 1, 10, 100, etc., ascending in the same 
direction. If a binary number of JN digits is shifted into 
the register (most significant digit first), the most sig- 
nificant digit will be shifted and supposedly doubled 
N—1 times, and the following digits will also be shifted 
and doubled an equally appropriate number of times. 

Obviously, the number in the register is not the deci- 
mal equivalent of the binary number and a quick exami- 
nation will show why. The appropriate power of 2 was 
to be computed by doubling at each shift, but when a 1 
was shifted from the “8” of any decade to the “1” of the 
next higher decade it only increased in value from 8 to 
10. (Remember, the assignment of values was arbitrary.) 
Accordingly, a value of 6 was lost somewhere, and if the 
computation is to hold true, the 6 must be replaced. 
Therefore, 6 is added to the decade from which the one 
was shifted. 

After the shifting is completed, only the necessity to 
clear out numbers over 10 in the decades by subtracting 
10 and adding 1 to the next higher decade remains, and 
the conversion from binary to decimal has been ac- 
complished. 

While the above process performs the desired conver- 
sion, it is quite cumbersome to instrument and several 
simplifications may be made. Remembering that shift- 
ing to the left doubles the number in the register, 3 may 
be added to a decade when a “1” is sensed in the “8” 
stage since the subsequent shift will transfer the “1” to 
the next higher decade and double the 3 to the neces- 
sary 6. 

To eliminate the possibility of a carry from the addi- 
tion propagating into the next higher decade with the 
resulting necessity for adding 6 again, and to eliminate 
the clearing operation at the end of the conversion, it 
is desirable to prevent a decade from ever containing a 
number greater than 10. If the number in the register is 
5, 6, or 7 prior to the shift, it will become a number 
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greater than 10 after the shift, and 5 should be sub- 
tracted from the decade and a 1 shifted into the next 
higher decade. In practice, this could consist of setting 
the “8” stage to “1” and subtracting 5, prior to the 
shift. This is the equivalent of adding 3, however, since 
+8—5 is equal to +3. A rule then develops: if the 
number in any decade is 8 or greater, add 3 before shift- 
ing, and if the number is 5, 6, or 7, add 3 prior to shift- 
ing. Or, in simpler terms, if the number in a decade is 5 
or greater, add 3 prior to shifting. Since the number 
after shifting can never be greater than 9, only five 
possible cases exist to which 3 must be added: 5, 6, 7, 8, 
or 9. There is no carry to a higher decade, and each 
decade can stand as a separate entity. 

It should be noted that there is no limit to the number 
of decades which may be cascaded. One decade is re- 
quired for each decimal digit in the final answer. 

The conversion from decimal to binary is based on 
repeated division of the decimal number by 2 with the 
remainders forming the binary number. With the deci- 
mal number entered in a shift register, divided in 
decades similar to that used for the binary to decimal 
converter, the division by 2 is performed by shifting 
toward the least significant end of the register. 

When a 1 is shifted from the “1” stage of decade to 
the “8” stage of the next lower decade, however, a di- 
vision by two is not performed since the value of the 1 
changes from 10 to 8 rather than 10 to 5. Obviously, a 
value of 3 has been gained, which must be subtracted. 
If, after a shift, a decade has a 1 in the “8” stage, it can 
only have come from the next higher decade and, ac- 
cordingly, the rule may be formulated: if, after the shift 
pulse, the number in a decade is 8 or greater, 3 must be 
subtracted. For obvious reasons there are again five 
possible numbers from which 3 must be subtracted: 8, 
9, 10, 11, and 12. There is no borrow to the next decade, 
and again each decade can stand as a separate entity. 
The number of shift pulses to be applied is equal to the 
number of binary digits in the largest number which can 
be contained in the register. 


BINARY-TO-DECIMAL CONVERTER 


The process of converting a number in binary form 
to its decimal equivalent consists of shifting the number 
into a shift register having 4N stages grouped to 
form WN decades, the contents of each decade forming a 
decimal digit in binary coded form, and the number of 
N decades required being determined by the decimal 
equivalent of the binary number. The conversion proc- 
ess consists of shifting the binary number into the 
register one bit at a time, most significant digit first, 
testing the content of each decade prior to each shift 
and adding binary 3 to any decade which contains 
binary 5 or greater. 

Fig. 1 is a block diagram of the converter, Fig. 2 is 
a block diagram of the logic circuitry associated with 
each of the decades, and Fig. 3 is an illustration of the 
operation of the converter. 
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100'S DECADE 10's DECADE UNITS DECADE 


8 4 2 \ 8 4 oe ' 8 4 2: ' 


SHIFT PULSE BUS 


BINARY FROM NH STAGE é 


a) 
he (2) SERIAL BINARY 


7 (3) 


REFL 
GREY TO BINARY 
BINARY 


DIODE LOGIC NETWORK DIODE LOGIC NETWORK 


M2 
:) TEST PULSE BUS 


DISABLE BUS 


Fig. 1—Block diagram, binary-to-decimal converter. 


aE. | 
| pewleada 
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Fig. 2—Binary-to-decimal logic network. 


BINARY CODED DECIMAL BINARY 


HUNDRED TEN 
Peviee eT 
ee a a 10. Us Leonor Ve Pol 
HARE eens oat of0% SHIFT 
mk a | (0 a | 1 ijo ove SHIFT 
BEERS Saee 10 0°1 SHIFT 
Ea SnSSne 1oot ADD 3 TO UNIT DECADE 
a i a a SWNT) SHIFT 
Eee L [fol oO! SHIFT 
EERE EES ol ADD 3 TO UNIT DECADE 
aes | [tfolt] ' SHIFT ~ 
Se anos oro ' ADD 3 TC TEN DECADE 
| | | [+] Jofofolo| SHIFT 
LL] [1] [olofofo] ADD 3 TO UNIT DECADE 
fmol onomro mn ayice 
\ 


8 SHIFTS COMPLETE THE 
CONVERSION 


Fig. 3—Example binary to decimal. Note: prior to each shift, 
each decade is tested for 5 or greater, and, if so, 3 is added. 


Referring to Fig. 1, there is shown a shift register 
which, for example, consists of 12 stages, S1 through 
S-12. The four stages S1-S4 plus network M1 form the 
units decade of the shift register, stages S5—-S8 plus 
network M2 form the 10’s decade, and stages S9-S12 
plus network M3 form the 100’s decade. Each stage of a 
decade, of course, will contain a binary “1” or a “QO” 
and the stages within the decades have binary weights 
of 8, 4, 2, and 1 increasing in significance in the same 
direction as the decimal weights. For convenience, the 
stages within a decade will be referred to by their binary 
weight; 7.e., the “4” stage would refer to stages S3, S7, 
or S11. The shift register is connected to a clock, or 
source of shift pulses, in such a manner that each shift 
pulse shifts the number in the register one stage to the 
left, or toward a more significant value. 
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_ A second output of the clock preceding each shift 
pulse, is the test pulse applied to the diode networks 
own in Fig. 2. In essence, the networks determine if 
ithe number in a decade is 5 or greater and, if so, add 3 
by complementing the appropriate stages of the decade. 
|For example: 
0101 (5) 
0110 (6) 
0111 (7) 
1000 (8) 
1001 (9) 


should become 
should become 
should become 
should become 
should become 


1000 (8), 

1001 (9), 

1010 (10), 
1011 (11), and 
1100 (12). 


(Examination will show: 

fie wand: 1? are 15,8,” “"4.” and: “1” stages are 
complemented. 

4 ands 2 dre 1 sand “1” is 0; “8,” “4,".“2 ” and 
“1” are complemented. 

Meo 1S loan “1” is-0, “2” and “1” are’ comple- 
mented. 

If “8” and “1” are 1, “4” and “1” are complemented. 


| 


_ The network contains the necessary logic to generate 
the proper pulses to be applied to the complementing 
inputs of the shift register stages. Many other networks, 
or adders, are available where the shift register stage 
oes not lend itself to being complemented. 

The clock contains a counter, and after a number of 
est and shift pulses equal to the number of binary digits 
in the word has been applied to the converter, the 
rocess is complete and the decimal number is in the 
converter in binary coded form. 

The description of the conversion has been for a bi- 
uary number in serial form. Switch 1 permits the con- 
wersion of a binary number in parallel form or a reflected 
binary number in parallel or serial form. 

_ For the conversion of a binary number in parallel 
form, the binary number is read in parallel into the 
shift register with the most significant digit in the most 
significant stage (S12 in the example) and SW1 is put 
in position 1 connecting stage S1 to S12. Application of 
shift pulses causes the conversion to be accomplished 
s previously described. To prevent digits of the binary 
number from being converted prematurely, the net- 
works are prevented from generating complementing 
ulses until the parallel-entered binary number has been 
shifted out of the decade. This is accomplished by 
eans of the “disable” bus shown on the network in 
ee 2. The clock generates disable gates during the 
counting period in such a manner that the units decade 
metwork is “enabled” after three shift pulses, the 10’s 
decade after six shift pulses, the 100’s decade after nine 
shift pulses, etc., the number three deriving from the 
fact that three shifts are necessary to enter a number 
greater than 5 into a decade. 
| To convert either serial or parallel Gray code to 
decimal, the same converter is used with the addition of 
a serial Gray-to-binary converter at the input to stage 
S1. The Gray code converter consists of an “and” of the 
input stage and the test pulse. If the input stage is a 
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or Decimal-to-Binary Converter 


HUNDREDS DECADE TENS. DECADE UNITS DECADE 
"—onoetq"n0#"———~ 7 i 
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SHIFT PULSE BUS 
OlODE ee NETWORK DIODE LOGIC NETWORK 
M2 
TEST PULSE BUS 


== DISABLE ous (iF REQUIRED! 


CONNECT TO 8y 
TO STORE 
BINARY NUMBER 
BINARY OUT 


DIODE Logic NETWORK 
MI 


Fig. 4—Block diagram, decimal-to-binary converter. Note: 
shift register shifts to right. 
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Fig. 5—Decimal-to-binary logic network. 
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Fig. 6—Example decimal to binary. Note: prior to each shift, 
each decade is tested for 8 or greater, and, if so, 3 is subtracted. 


“1,” a pulse is generated which complements either a 
pre-input flip-flop, if the input is serial, or the last stage, 
if the parallel storage method is used. 


DECIMAL-TO-BINARY CONVERTER 


The process of converting a number in decimal (bi- 
nary coded decimal) form to its binary equivalent con- 
sists of shifting the number out of a shift register having 
4N stages grouped to form N decades, the contents of 
each decade forming a decimal digit. The conversion 
process consists of shifting the binary number, one bit 
at a time, least significant digit first, out of the register, 
testing the content of each decade after the shift and 
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subtracting three from any decade containing eight or 
greater. 

Fig. 4 is a block diagram of the converter, Fig. 5 is a 
block diagram of logic circuitry associated with each 
of the decades, and Fig. 6 is an illustration of the opera- 
tion of the converter. 

Referring to Fig. 4, there is a shift register and clock 
arranged in a manner similar to that of Fig. 1 except 
that the shift is toward the least significant digit, and 
the decimal number in binary coded form has been en- 
tered in the register in parallel form. A second output of 
the clock, following the shift pulse, is the test pulse, 
which is applied to the diode networks shown in Fig. 5. 

The networks determine if the number in the decade 
is 8 or greater and, if so, subtract 3 by complementing 
the appropriate stages of the decade. For example: 


1000 (8) should become 0101 (5). 
1001 (9) should become 0110 (6). 
1010 (10) should become 0111 (7). 
1011 (11) should become 1000 (8). 
1100 (12) should become 1001 (9). 


Examination will show that: 


Lio uisdiand.°47.and “1”. are 0;%8,"."4."and'“1 are 
complemented. 

Mis giswlood Mente and 24” sand<2”-are0, 48)" s°4,” 
“2,” and “1” are complemented. 

If “8” is 1. and “1” is 1, “2” and “1” are complemented. 

If “8” is 1 and “1” is 0, “4” and “1” are complemented. 


The network contains the necessary logic to generate 
the proper complementing pulses. 
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After a number of shift and test pulses equal to the 
number of bits of the largest binary number which could 
be represented by a decimal number of WV digits the 
conversion is complete, and the decimal number has 
been shifted out of the converter in serial binary form, 
least significant digit first. 

If the decimal number is in serial form, it may be 
shifted into the register prior to the conversion. If the 
binary number is desired in parallel form, switch SW-1 
is thrown to position 2 and the binary number circu- 
lated back into the shift register. However, disable 
gates must be generated to prevent complementing the 
binary number. 


CONCLUSION 


The circuit design of either of the converters is 
straight forward and there are many techniques avail- 
able to provide the shifting and testing functions. The 
model built used flip-flop shift register stages and dc 
logic. As the stages complemented in the testing opera- 
tion affect the logic, it is necessary to delay the stage 
inputs to the logic by approximately the width of the 
test pulse. 

Further study indicated that theoretically, the speed 
of conversion could be doubled at the expense of some 
complexity by eliminating the “test” operation. Essen- 
tially, the presence of a 5 or greater in a decade would 
block the shift pulse to the “8,” “4,” and “2” stages of 
that decade, and these stages would be set to the num- 
ber plus 3 as it would have been shifted. Using 1-mc 
stages, a 20-digit binary number or 6-digit decimal num- 
ber could thus be converted in 20 usec. 


Diodeless Magnetic Shift Registers Utilizing 
Transfluxors* 
NOAH 8S. PRYWES} 


Summary—Shift registers using magnetic cores have been con- 
ventionally designed with diodes to provide isolation between 
stages. It is possible to use transfluxors in a two-core-per-bit 
fashion so as to provide isolation between stages by magnetic means 
rather than with diodes. A design procedure for doing so is discussed. 


I. INTRODUCTION 


AGNETIC materials with a characteristic as 
Vi shown in Fig. 1 can be used for information 
storage. Their main useful features are two 
opposite remanent magnetizations at saturation and a 
minimal coercivity H, required for a domain change in 


a process of magnetization reversal. 
To explain further the magnetic characteristic, con- 


* Manuscript received by the PGEC, September 3, 1958. 
_ t Moore School of Elec. Eng., Univ. of Pennsylvania. Formerly 
with Remington-Rand UNIVAC, Philadelphia, Pa. 


sider, for instance, an initial situation of +B, and H=0 
in Fig. 1. For a positive or a negative field intensity 
smaller than H,, the nearly horizontal line at the top of 
the characteristic in Fig. 1 will represent the corre- 
sponding changes. At the removal of such a field, re- 
turn to the magnetization of +B, and H=0 will occur. 
Thus any change of flux at the application of such a field 
was reversible. However, if a negative field in excess of 
H, is applied, an irreversible flux change occurs along 
the nearly vertical part of the left of the major hysteresis 
loop shown in Fig. 1. The average rate of change of 
such flux along this line is approximately proportional to 
the amount by which the applied field exceeds the 
coercive field H,. If the applied field is removed before 
the complete reversing of magnetization along the major 
hysteresis loop, the return to H=0 will occur along a 


| 
| 
B 
| 


Te 


Fig. 1—Magnetic characteristic. B = flux density, H=field in- 
tensity, H,=coercive field. 

minor hysteresis loop corresponding to the average 
remanent flux density at the instant of removal of the 
applied field. This description of the magnetic char- 
acteristic should be sufficient for the purpose of this 
Hiscussion. Other secondary phenomena will be as- 
sumed to have a negligible effect. 

As bistable devices, magnetic cores may be assembled 
to constitute a shift register. It has been found that 
certain isolation must be provided between stages; this 
so far has been available through the use of diodes. [1] 
(The invention of the transfluxor [2],[3] reveals a 
method by which the desired isolation can be achieved 
without diodes. This paper will discuss some registers 
in which such advantage is taken. Fundamentally they 
are of the “two-core-per-bit” kind, [1],[4] with the 
isolation between stages provided by magnetic means 
rather than diodes. A design procedure for a shift reg- 
ister utilizing cores and windings only will be discussed. 
Production costs for such a register would be low with 
printed-circuit techniques. 


II. A THREE-APERTURE TRANSFLUXOR 


A short description of the operation of the transfluxor 
will be given here, but a complete explanation and de- 
scription is given by Rajchman [2], [3]. 

The transfluxor shown in Fig. 2 can be used for a 
stage of the proposed shift register. It is a core which 
consists of four legs with a small aperture between legs 
1 and 2, a small aperture between legs 3 and 4, and a 
large aperture between legs 2 and 3. It is made of ferrite 
material which has the loop characteristic shown in Fig. 
1. In addition to its characteristic described in the In- 
troduction, the requirement of continuity of flux lines 
(div B=0) will determine any flux changes. In order to 
change the flux around the big center hole, the ampere 
turns of the winding must produce a magnetic field ex- 
ceeding H, along that long path. Fewer ampere turns 
are needed to change the flux along the much shorter 
paths around the small holes. Thus it is possible to 
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Fig. 2—Three-aperture transfluxor. 


change the flux around one small hole without affecting 
the windings through the other small hole, provided 
enough ampere turns are applied to exceed H, in the 
short path around the small hole, but fewer than are 
required to reach H, along the longer path around the 
large hole. The required isolation is achieved, in that 
information can be read out by flux change around the 
right-hand small hole without affecting the input wind- 
ing through the left-hand small hole. When a shift of 
flux is required, a larger pulse applied to the input wind- 
ing through the left-hand small hole reads information 
into the area around the right-hand small hole. 

The transfer of a “one” and of a “zero” is shown at the 
bottom of Fig. 2. The arrow below each of the legs of the 
transfluxor indicates the direction of remanent flux at 
saturation in that leg. The arrows show the direction of 
the remanent fluxes after application of the pulse 
named. ; 


Transfer of a One 


The first line shows the remanent fluxes after a block 
pulse. This situation also corresponds to a zero state. 
The flux polarities are downward in legs 1 and 2 and 
upward in legs 3 and 4. The condition of continuity of 
flux in the ferrite is thus satisfied. 

A set pulse, representing a one, causes a clockwise 
field around the left-hand small hole. This field is of 
sufficient amplitude to exceed H, in the path round legs 
1 and 3. The flux in leg 2 can change only slightly since 
it is already in saturation in the downward direction. 
The change in flux in leg 1 must then be compensated 
by an equal and opposite flux change in legs 3 and 4. 
The path of flux change through leg 3 is shorter than 
that through leg 4, and thus will be preferred. Also, leg 
4 can be biased to prevent a change of flux in it. A flux 
change in leg 3 then necessarily compensates the change 
of flux in leg 1 to satisfy the continuity of flux condition. 

A drive pulse, applied after the set pulse, changes the 
flux around the small hole on the right, producing an 
output pulse. The input winding, however, is not af- 
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fected because of the limitations on the amplitude of the 
drive ampere turns. 


Transfer of a Zero 


The absence of a set pulse represents a zero. In the 
absence of a set pulse, the core remains in the block 
state. When a drive pulse is applied, the field caused can- 
not affect the flux in legs 1 and 2 because it does not 
exceed the value of H, along the path around the large 
center hole. It cannot change the flux in leg 3 as this is 
already at saturation in the direction of the applied field. 
No significant flux change in legs 1, 2, and 3 is possible. 
Hence, no flux change can occur in leg 4 because any 
change would have to be compensated by equal and oppo- 
site changes in other legs to satisfy the condition of flux 
continuity. Thus there is no significant change in flux 
linking the output winding, resulting in lack of input to 
the following core, or a zero output. The process of 
transfer and discrimination between a zero and a one 
has thus been demonstrated. 


Ill. A DiopELESs Two-CorE-PER-BIT SHIFT REGISTER 


Before analyzing the operation of a transfluxor stage, 
it would be well to present a general description of how 
these cores may be connected to form a shift register. 
Figs. 3 and 4 illustrate the operation of a two-core-per- 
bit shift register. The register shown at the top of Fig. 3 
is capable of holding four bits. The arrows show the 
pattern of remanent fluxes in all stages for the process 
of shifting “1100” into the register. Fig. 4 diagrams the 
timing of supply pulses and some input pulses. 

The cores in Fig. 3 are shown with rectangular aper- 
tures only for convenience of presentation. Legs 1, 2, 3, 
and 4, input, output, drive, cancellation, and block 
windings are labeled on the first core of the left of Fig. 
3. The windings on the other cores are similar. The out- 
put winding of each core is connected to the input 
winding of the successive core. The winding sense of all 
windings and the series grouping of windings, so that 
they may be fed by a single current pulse, are also 
shown. 

The windings differ somewhat from those in Fig. 2. 
The input goes to leg 2 instead of leg 1 in order to 
shorten the path of flux changes during a one input. 
When the input winding is on leg 2, an input pulse 
reverses the flux in a path surrounding only the large 
center hole. On the other hand, when the input winding 
is on leg 1, as in Fig. 2, the input pulse changes the flux 
around the longer path which surrounds both the large 
center hole and the small hole at the left. Placing the 
input winding on leg 2, therefore, requires fewer ampere 
turns in the input. The drive winding is wound on leg 3 
rather than leg 4 in order to lessen air flux linkages be- 
tween drive and output windings. A cancellation wind- 
ing, wound in opposition to the block winding, has been 
added. This cancellation winding raises the upper limit 
on the drive-pulse ampere turns. 

The lower part of Fig. 3 diagrams the process of 
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shifting a pattern of information. Again the arrows un- 
der each leg represent the direction of magnetization 
at saturation. With the timing shown in Fig. 4 a source 
of four current pulses is necessary. These pulses are 
Drive I, Block I, Drive II, and Block II. Since this is 
a two-core-per-bit register, Drive I and Drive II pulses 
are applied to alternate cores. The block pulses follow 
the drive pulses and reset the cores to the zero state. 
Block pulses are required because the transfluxor, by 
contrast with the single-aperture core, has nondestruc- 
tive readout. This seeming disadvantage can be turned 
to advantage if an additional small hole similar to the 
other two small holes is provided, as shown in Fig. 5. 
Then reading out without shifting the information is 
possible. 

A pattern of 1100 is shifted into the register to illus- 
trate its operation as shown in the lower part of Fig. 3. 
The one input into the first core is in the form of a posi- 
tive pulse at the time of Drive II, followed by a negative | 
pulse at the time of Block II. (See Fig. 4.) The output 
from a stage into the successive stage follows the same 
pattern. Thus, if the output of the eighth core is con- 
nected to the input of the first core, the information is 
propagated around a ring. The flux patterns of cores in 
a one stage are framed to assist in examination of the 
shifting flux pattern which represents the information 
introduced. 

The shift register may be modified to accept an input 
of all positive pulses, as shown by the alternative input 
wound on leg 1 in Fig. 3. 

At Drive II the one enters the first core and flux 
changes in legs 2 and 3. The Block II pulse, which fol- 
lows, reverses the direction of flux in legs 1 and 2. At 
Drive I the one enters the second core where flux 
changes in legs 2 and 3, while flux changes on legs 3 and 
4 of the first core. Block I reverses the flux in legs 1 and 
4 of the first core and legs 1 and 2 of the second core. 
On the next Drive II pulse the second one in the 1100 
pattern enters the first core of the register, and the first 
one in the 1100 pattern is shifted to the third core. 

With this concept of how information is shifted, a 
closer look can be taken at a single step in this process. 


IV. ANALYSIS OF INFORMATION TRANSFER 
Analysis of Flux Switching during Drive Pulse 


In the three cores shown in Fig. 6, information is 
transferred from left to right. This analysis deals with 
a cycle in which information is transferred from the 
middle core. Supply windings are shown on the middle 
core only. The block winding with Vz turns provides a 
magnetomotive force (mmf) of B in the clockwise di- 
rection around the large aperture. The drive winding 
on leg 4, with Np turns, provides an mmf of A in the 
clockwise direction around the small aperture at the 
right of the transfluxor core. An additional cancellation 
winding with Ng turns provides an mmf K in the 
counterclockwise direction. The K mmf is provided to 
cancel the effect of the A mmf around the large aper- 
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Fig. 4—Supply and input pulse pattern for the register of 
Fig. 3. 
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Fig. 5—Four-aperture transfluxor with nondestructive 
interrogation. 
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Fig. 6—Mmf’s and their paths in information transfer. 


ture; that is, to confine the effect of A to legs 3 and 4 
only. Hence K must be coincident with A and is sup- 
plied by arranging the cancellation winding in series 
with the drive winding. The mmf’s Q and J are pro- 
vided by the output winding on leg 4 with NV, turns and 
the input winding on leg 2 of the input core at right with 
N; turns, respectively. The box connected in series with 
the output winding and the input winding of the core 
at left represents the resistance R and the self-induct- 
ance L; of the input and output windings. 
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The process of reading a one from the middle core t 
the core at the right is illustrated in Fig. 6. The arrow 
in the diagram again represent direction of magnetiza 
tion at saturation. The line marked “one in” shows the 
state of the core prior to the drive pulse. The line marked 
“one out” shows the flux pattern of the core after the 
drive pulse. The line marked “block” shows the flux 
configuration after application of the block pulse. I 
the case of a zero transfer both cores will remain in the 
block state. 

To simplify the analysis it is assumed that the wall 
thickness is small enough to justify considering a uni- 
form magnetic field in the core walls. Thus the mag- 
netic paths and their lengths are depicted by lines Ly, 
L, L3, L4, and LZ on the left-hand core in Fig. 6. The sub 
scripts refer to the respective legs. 

When the drive pulse is applied, it is desirable that 
the flux be switched in legs 3 and 4 of the transfluxor 
core in the case of a readout of a one and that no ap- 
preciable flux change occur in the case of a zero. These 
requirements can be set into formulas with the aid of 
Fig. 6. x 

As seen in the flux diagram of Fig. 6, the following 
condition is imposed on operation during the drive 
pulse. In reading a one out of the center core, the path 
enclosing the large aperture and the small left aperture 
should not be affected. Hence the fields in this path 
should be smaller than H,. That is, 


ie Hdl < H,(2L + Ly). (1) 

Ls 

When its direction is upward, H is taken as positive. 
Equal magnetic fields are required for equal rate of 
change in induction at any point in the core. Since the 
flux is being switched at an equal rate along the path 
L; and Lu, the following can be stated: 


Hdl = H;3L3; 
Lz 
Js lGHL, = Hills, & H; = Hi> OZ, (2) 
Ing % 
and from (1): 
K < H.(2L + Ly) + H3Ls. (3). 


Also, the mmf’s in the path in which there is a flu 
change (legs 3 and 4) follow the relationship: 


Als + Hal, = Hz(L3 + Ls) = A — Qyne. (4) 


Let Ma,3 be the mmf existing in the path of legs 1 and 3 
when the field in this path is H,: 


Mes = H(Li + L; + 2LZ). 
Eqs. (3) and (4) can be written in the form: 


Konan = Me1,3 a L;3(H3 aay HT) (5) 


Aviax ar Qone = H3(L3 =| La), (6) 


/ 


| 

Both Qone and H; depend on load conditions and the 
ted of operation. Consider the case of the output con- 
tted in shunt to m input cores as ‘in Fig. 7(a). The 
trent in input and output windings is the same. The 
.e of change of flux and the fields in legs 2 and 3 of the 
ut cores are equal. Since the cross section of mag- 
He material in the sections denoted by lengths Z can 
‘made much greater than that in legs 2 and 3, the 
lied field can be essentially H,. Thus, we can write 
(expression FOU ana: 


AEA eee peat, A eG, 
N, iets Vi 2( 2+ 3) + a] (7) 


here the value of Hy, relates to leg 2 of the input core. 
Meo3=H(L2t+L3;+2L). Then (7) can be written 
the form: 


No 
Qone =n WN: [M.2,3 == (H. = H.)(L2 + L;)]. (8) 
In the case of a zero transfer, the drive pulse does not 
ct the cores. They remain in the block state. Again, 
ks requires that the field in the path around the large 
rture be less than H,. That is, 


tae Grero her K < A Ls + ZL =i L2) 


AS yer — K > 0 and Qero < one: (9) 


LA for Mosz=H-(Lo+L,+2L) and using the maxi- 
im values of A, K, and Qrero, 


ee a Tax Pag ees = M 2,4. (10) 


In the case of a zero transfer, the real requirement is 
» absence of flux change in the input core at the right 
Fig. 6. Even though the flux in leg 4 of the middle 
re changes slightly, producing a small pulse in the 
‘tput winding, this does not amount to an information 
insfer unless the flux also changes in leg 2 of the input 
re. Any small flux clrange in leg 4 of the output core 
erased by the following block pulse without affecting 
2 operation. To cause a flux change in the input core 
2 mmf Qzero Must provide a field in excess of H, along 
2 path consisting of Z2, L, L3, and L. The significance 
‘this observation is that in the case of a zero the flux 
ay change in leg 4 of the output core. Such a flux 
ange would be compensated by the voltage drop in 
2 resistance R and self-inductance L;, without causing 
‘irreversible flux change in the input core. The maxi- 
im mmf Qyero allowed in the case of a transfer of a zero 
puld then be: 


Ne (a1) 


Jecrotuax) eer, M 2,3 


i 
1ere n is the number of inputs of cores connected in 
unt to the output of the driven core. 

The rate of flux change in leg 3 of the output core is 
oportional to the difference between H; and H,; that 
leg 2 of the input core is proportional to the difference 
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between Hz and H,. Hence, the relationship between 
H; and His 


N (3 — H.) ='V.(A2 — Z.). (12) 


It is desired to find the upper bound on A. An exam- 
ination of the conditions leading to (11) shows that A is 
maximum when K also is maximum. Hence, in design 
the number of turns on the cancellation winding should 
always be such as to provide Kmax. There are five equa- 
tions [(5), (6), (8), (10), and (11)] with the five un- 
Enews aA ax, | Dex, (lis), regs and + Wins.  L eyceans ve 
solved for Amax and Kmax which are required as design 
information. Since Anax and Kmax vary somewhat for 
differently shaped cores, it is of little interest to write 
these expressions. However, certian approximations 
can be instructive. When legs 1, 2, 3, and 4 are ap- 
proximately equal, and where Z is large and n=1, 
Amax Can be approximated to be 


Amax = 3M ¢2,3. (13) 


The upper limit of drive mmf (A) imposed by (13) 
limits both the number of transfluxors (7) that a single 
transfluxor can drive and the speed of switching. When 
L, Ly, Ls, L3, and L4 are small and switching is extremely 
slow, » can be calculated. A minimum of one transfluxor 
and a maximum of ” would operate satisfactorily. When 
switching is very slow and n,/N;~1 then, 


JBla S ily Seb 
Qsero ee M.2,3 
Dire = nM -2,3. 


This shows the solution for mmax from (5) through (12), 
Nmox = 3. (14) 


To determine the switching time of a transfluxor driv- 
ing only one other transfluxor, as in a shift register, as- 
sume 


DIN=iL1 =[,= 13 = Ly 


then, 
Hp ad Az. 


Introducing these assumptions into (5) through (12): 
H; =) He = vale. (5) 


In an RCA transfluxor, a switching field of three 
times the coercive field provides 0.3-usec switching 
time. The time can be improved by making Z much 
greater than Li, Lx, Ls, and Ls. A shift register using 
RCA transfluxors has been operated above 500 kc with 
a switching time of 0.8 usec. 


Turns Ratio of Output to Input Windings 


To determine the ratio of input to output winding 
turns, the process of reading out a one will be analyzed. 
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Fig. 7(a) shows the transfluxor parts involved. Only 
legs 3 and 4 of the output core and legs 2 and 3 of the 
input core are shown since there is no flux change in the 
other legs. This simple situation, the transfer of infor- 
mation in a single-aperture core, has been treated pre- 
viously by A. Wang [1]. Fig. 7(a) shows the flux con- 
figurations before and after a drive pulse. 

In Fig. 7(c) the magnetic characteristic is repre- 
sented by an equivalent circuit. The requirement that 
field H, must be exceeded before flux change occurs can 
be represented as a current drain of M,/N,. The amount 
by which the applied field exceeds this value is used to 
accelerate the flux change. This excess can be repre- 
sented by an equivalent resistance which is propor- 
tional to the square of the number of turns [5]. In the 
figure, the current drain of (M,3,4)/N. represents the 
coercive current in the output core, and 2(M.2,3)/N; 
represents the coercive current in ” input cores. Also, 
the applied drive field reverses the flux in input and 
output cores at the rate determined by the equivalent 
resistance. This resistance is Ro times the square of the 
number of turns. 

The ratio, NV,/N;, should be greater than one in order 
to compensate for flux lost in the transfer. Also, V./N; 
has an upper limit, a result of the requirement that 
drive field must exceed coercive drain. 


A > M34 M 2,3 ; (16) 
No Nis N; 
This, together with the lower limit, forms: 
No A aa M 3,4 
dx ——————_ ; (17) 
N; nM -2,3 


The rate of flux reversal in the input cores depends on 
energy supplied to the resistance (V;?R,)/n. See Fig. 
7(c). This energy is maximized if N,/N; is near enough 
to unity to compensate for flux loss. In practice, such 
refinement may require so many turns that it rules out 
printed circuit techniques. In the shift register that was 
operated, the number of turns used was Ny=3 and 


N;=2. 


The Blocking of the Buildup of Zeros 


Fig. 7(a) illustrates the reversible and irreversible 
flux changes. The drive winding is magnetically linked 
to the output winding through the magnetic core, which 
is assumed to have only irreversible flux changes. 
Separate linear linkage m accounts for reversible flux 
changes. During the transfer of a zero, voltage V, is 
induced through the linear linkage m. The amplitude of 
voltage V, is 

dla 


Vo=m — 


dt te) 


where Jq is the drive current. If voltage V, causes an 
mmf in the input winding in excess of M...3, a small 
change of flux may occur in the path of legs 2 and 3 of 


IRE TRANSACTIONS ON ELECTRONIC COMPUTERS 


December 


INPUT CORES 
Lo L3 


OUTPUT CORE 
Ls La 


(a) aii 3 
4 : we ls = 

} A & } {J} } 

DRIVE ; 


2 


‘« 
DIRECTION OF FLUX 
: 
L 


eae | a 


ove tS 1g 


R Ly 
4, AM 100% 
Mosq Me23 N,2R 


; : : A : 

(c) SiNGae i Ro aa Fare n 

Fig. 7—Transfer of a one. : 

the input core. Successive drive pulses amplifying this. 

small change may eventually introduce a one. There are 

several methods for avoiding this undesirable process. | 

Resistance in the output circuit: A simple way of 
blocking zeros is to introduce a resistance, R, in the 

loops of the input and output windings such that 


ig 
Mee, > — N;. 19 
2,3 R ( ) 


The disadvantages of this method are 1) it compli- 
cates production, thus raising costs and 2) the resistors 
dissipate energy thus reducing speed. 

Use of a buck core: The use of a buck core is known as 
well [4],[6]. A corresponding circuit configuration is 
shown in Fig. 8. The buck core ™m this circuit is driven. 
at the application of each drive pulse inducing in each 
output a voltage in opposition to that induced by the 
transfer of a zero. During transfer of a one the blocking” 
pulse is small and of short duration so as to subtract. 
little flux from the output. Only the output circuits are. 
shown in Fig. 8, with the rest of the circuit remaining 
the same as in Fig. 6, | 

Cancelling zero output during block pulses: Flux 
change due to zero output can be cancelled during the 
block pulse. This is done by reversing the flux in the 
path of leg 2 and 3 of the input core by an amount equal 
and opposite to the change caused by a zero output. 
The block pulse must then induce the corresponding 
voltage in the output circuit. One means is to provide 
sufficient linkage between the block pulse circuit and 
the output circuit and to control the rise times of the 
drive and block pulses. This approach has been found 


very simple and satisfactory in operating the shift 
register. 


ORIVE 2 


_ Fig. 8—Suppressing zero output by the use of a buck core. 


BLOCK I 


Fig. 9—Another arrangement to cancel zero buildup. 


| An arrangement where the linkage between the out- 
ut circuit and the block pulse circuit is provided in the 
put core is shown in Fig. 9. 

| Reducing linear flux linkage: By winding the drive 
td output windings on legs 3 and 4, respectively, in- 
ead of both on leg 4, the linear flux linkage coefficient 
| is decreased greatly. The two ways of winding with 
‘e correct number of turns are shown in Fig. 10. The 
me number of ampere turns will link each of the 


»ertures in the two ways of winding if the following 
jationships are observed: 


| A=A'’; K’=A—K; Q=%. (20) 


Thus the previous analysis is valid for both ways of 
inding. In Fig. 2 the drive and output windings are on 
fferent legs. 

The validity of the analytical considerations in Sec- 
on IV has been verified experimentally. 


V. CIRCUIT VARIATION 


An increase in gain may be achieved by using a wind- 
g through the large aperture of each core that is 
orted during the drive pulse applied to the core but 
yen during the drive pulse applied to the previous 
age, causing information input into the core (Fig. 11). 
he upper limit on the drive pulse then is removed. 
he shorted windings are arranged in series alternately 
1 odd and even stages and normally short all of the 
res in the forward direction of the diodes. However, 
iring the drive pulse which reads information into the 
re, the short is removed by the pulse created by the 
yrresponding drive pulse. The process requires four 
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ings both on Leg 4, (b) Drive winding on leg 3, output winding 
on leg 4,£,A=A’, K’=A—K, 2=0’. 
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additional diodes per shift register, which in themselves 
do not add much to the component count. 


VI. SHirt REGISTERS AND DELAY ELEMENTS 
MADE OF A STRIP OR A RING OF 
MAGNETIC MATERIAL 


Production may be simpler and cheaper if the cir- 
cuitry is composed of a strip or a ring of magnetic ma- 
terial rather than of individual cores. This is particu- 
larly so when utilizing printed circuit techniques. Some 
circuits of this type are discussed here. 

The individual cores of the shift register of Fig. 2 can 
be laid out as a strip with spacer apertures providing 
separation between the cores. Such a system is shown 
in Fig. 12. At the bottom of this figure the flux con- 
figuration for transfer of a one is shown. In the system 
of Fig. 12 the block state in alternative stages is repre- 
sented by saturation in opposite directions. Comparison 
with Fig. 2 shows the significance of the difference of 
the flux pattern. 

This type of shift register can be connected in a ring 
or made out of a ring, connecting the output back to the 
input. 

Fig. 13 shows a scheme which eliminates the necessity 
for input and output windings between stages. Also, 
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Fig. 13—A strip shift register without input and output windings. 


less magnetic material is involved per stage. Even 
though the losses may be lessened in this way, it is 
likely that some flux will be lost in the process of shift- 
ing information. Since no flux gain can be expected 
from the turns ratio between output and input wind- 
ings, the operation may not be stable with information 
rotated in a ring. However, information may be stepped 
up through many stages before satisfactory discrimina- 
tion between ones and zeros is lost. 

This type of shift register would be particularly sat- 
isfactory then for use as a delay line. Furthermore, for 
this use only two drive pulses are required. A ring that 
may serve for this purpose is shown in Fig. 14. A further 
advantage of such systems is that a bit of information 
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transferred corrects the nonregularity in flux pattern 


remaining in the input stage after turning once around 
the ring. In doing that, the ring is returned to the 
blocked stage. This is illustrated at the right of Fig. 14 
where the corresponding flux pattern is drawn under 
the spread ring. 
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CORRECTION 


M. W. Marcovitz and E. Seif, authors of “Analytical 
Design of Resistor-Coupled Transistor Logical Cir- 
cuits,” which appeared on pages 109-119 of the June, 
1958 issue of these TRANSACTIONS, have requested that 
the following corrections be made to their paper. 

Eq. (17) on page 113 should be 


To1 


cul . 
b1 


V owt = 


On the same page, (21) should be preceded by the 
phrase “From Fig. 5.” 
The last term of (35) on page 119 should be 


+ i , 
The squared term in the numerator of the term on 
the left side of (36) should be J, and not 1. 


nV ce {= 
BCLs 


j 
q 


orrespondence 


Iding of Symmetric Functions* 
INTRODUCTION 


In this note, T(X, - + + , X,) isa Boolean 
ction of m binary variables X), +--+, Xn, 
resenting the transmission between two 
minals of a contact switching network. 
h a function T is said to be symmetric on 
|and X; provided 


‘ OMECES 9 CP, CD. Gr yen 


Xj, Xj, XG; hes Xn) 
TX, ae Xi, Xj, Xin, + ee 
| Xj, Xi, Xjuty re re Xe 


| T is said to be symmetric provided T is 
metric on all pairs of variables X; and 
in Xy,-++, X,. Following the conven- 
on of Shannon,! a symmetric function T 
i be denoted by Sp,,95,° * * p,,(X1, °* * » Xn) 
'T=1 only when any #; variables are 1 
hd #— >; variables are 0, fori=1,---,7. 
It is well known that a standard contact 
ttwork exists to produce any symmetric 
iction.!2 For example, in three variables 
fe circuit is as shown in Fig. 1. 
The terminals marked 0, 1, 2, and 3 pro- 
(ce transmissions from terminals a of So, 
, S2, and S;, respectively. We will say that 
- row of relay contacts with output termi- 
(1 z will be called the ith row, and that the 
Jumn whose output terminal is 7 will be 
ed the ith column. Shannon! has used 
se concept of folding for the purpose of 
mplifying symmetric functions Sp,,p5,° * * .p, 
fi - , Xx). 
Suppose this function can have its cir- 
jit replaced by one in which all contacts 
wove the $.41 row and on that row can be 
sted, with the remaining contacts 
bove the ~, row being connected so that 
fe output of X; is connected to the junction 
_X;’, Xiu, Xiu’, Xj on some row below 
‘e ps row, where j=ps+1,---, n. 
If this new circuit has for output 
17° *,p,(X1, +++, Xn), the original circuit 
‘i be said to fold down to the psth row. 
~ Furthermore, a circuit for Sp,,n9,° * *.7, 
Ki, ° + * ,Xm) is said to fold completely pro- 
‘ded it will fold down to the f: row; here 
‘is supposed that pi<po< +++ <$y. 
Shannon? and Caldwell! both remark that 
folding process will apply to any symmetric 
‘nction in which the subscripts are in arith- 
‘etic progression. This note gives a neces- 
ry condition for folding which includes 
larger class of subscripts. 


| THEORY 
heorem 1 


If the symmetric function Spas 
X,°-+, Xs) folds down to the pth row 
»><q<s), then there is an integer m such 
nat 


* Received by the PGEC, November 14, 1958. 

1C, E. Shannon, “A symbolic analysis of relay 
1d switching circuits,” Trans, AIEE, vol. 57, pp. 
3-723; 1938. 

2S. H. Caldwell, “Switching Circuits and Logical 
esign,” John Wiley and Sons, Inc., New York, N. Y.; 
58. 


3 Op. cit., p. 721 
4 Op. cit., p. 250. 
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Fig. 1. 


Big. 2. 

Someee. 
=m. 

Gia? 


Proof: Consider Fig. 2. It is known that 
none of the relays above the gth row and to 
the right of the sth column are needed, nor 
above the pth row and to the right of the 
gth column. 

The first fold will be accomplished by 
connecting node 1 to node 2. It can be ob- 
served that S, could originally be obtained 
only by going up the sth column. After fold- 
ing, g+1 open relay contacts are encountered 
up to the g+1 row. From node 2 up to the 
qth row are s— (q+1) more open contacts, so 
that S, is still possible, but now at the gth 
row terminal, together with Sz. 

The second fold will require the output 
of Xp41 on the sth column to be connected 
to some 7 column. It will connect just below 
Xp42. Since there are g—p contacts on the 
sth column above X41, after the first fold, 
the 7 column must end at Xp424[¢-@1)) 
=X. But then either g+1=s (in which 
case no folding is possible) or else in the fold- 
ing process, the 7; column must have its top 
contact connected to the 72 column. In this 
fashion, since this is a finite system, we can 
see that the sth column is connected to 
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the 7; column, that to the rz column, ++:, 
the 7m_2 to the 7m_1, the 7m_1 to the pth. The 
number of open contacts on each of the 7; 
columns (¢=1, +--+, m-—1) and on the pth 
column above the connection from the 7:1 
row and below the pth row is the same and 
is equal to g—(p+1). Then for Ss to be ob- 
tained at the pth row terminal, there must 
be a path starting at a on which there are 
exactly s open contacts. On the sth column 
there are p below the pth row. Above the 
pth row on all columns is a total of m open 
contacts. Hence, 


s—p=mq—(+1))+m 


or 
See? 
Cie. 
It can be observed that when m=2 the sub- 
scripts are in arithmetic progression. 
Corollary lie Tis Sota Caen 
folds to the pth row (where p<q<s<n) 
then there exists an integer m such that 


Sect 
Gia? 
. Proof: S; has as one minterm the quan- 
tity 
X1X2, “ectar XX oy1' Xey2', bab FS 9 DOME 


This minterm must be realized by the con- 
tacts on the mth column. After folding to the 
pth row, the first +1 open contacts are 
found on the nth column. Furthermore, that 
minterm can be realized only by a path 
up the mth column to the pth row. Then as in 
Theorem 1, the output of Xp, must be 
connected to the 7 row, that to the rz row, 
etc. Again, the part of the minterm 
Xp42, °* * »Xqu1, can be realized only by a 
path up the ™ column (since a move to the 
right would introduce X,’ for b<s); and 
similarly for the 72, 73,°+*, %m columns. 
The top contact on the 7 column (below the 
pth row) will be Xs. If it were not, the 
path X,--+>, Xsy,-+-,Xn’ would be 
open contrary to hypothesis. Hence, as in 
Theorem 1, s—p=m(q—(p+1))-+m and 


Sree 

Cie 
Coreary2: li Ss ao, Ane aaa) 
where pi<po<:+:<p,<s folds completely, 


then there exist integers m,7=0,1,2,---, 
r—3 so that 


=m. 


= mM. 


=. 


Pri an pri-2 A 
—_-—— = mi. 
Pri — pr—i-2 
Proof: It is known that 
Spine 92st maAG, sg pps) 
= Spy .p9.03(X1, cae Xs) +Spopyp,(X1, ere ser X;) 
Siete ’ +Sp,_op,_ 17,15 se Ta | Xs) 


If the function <S5).r...<9;X ites, oe) 
folds, then each of the functions in the sum 
on the right-hand side of the equation must 
fold. Applying Theorem 1, the result follows. 
G. P. WEEG 
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f the Institute for Computer Research. 
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| Takeo Miura was born in Kyoto, Japan, 
p 1925. He graduated from Kyoto Uni- 
‘ersity with the degree of Bachelor of Engi- 
fering (Kogakushi). He joined Hitachi, 
itd., Tokyo, in 1949, and was engaged in 
ine design of induction motors of Kameido 
Vorks. In 1950, he was transferred to 
{itachi Central Research Laboratory, where 
22 has since been working on the develop- 
went of electronic analog computers, regu- 
ttors, and servomechanisms. 

| Mr. Miura is a member of the Institute 
f Electrical Engineers of Japan. 


9GEC News 


CALL FOR PAPERS 


The 1959 Solid-State Circuits Conference, sponsored 
xy the IRE Professional Group on Circuit Theory, 
sIEE Committee on Electronics, IRE and AIEE 
>hiladelphia Sections, and the University of Pennsyl- 
ie, will be held on February 12 and 13, 1959, in 
hiladelphia, Pa. Plans for the 1959 Conference call for 
he presentation of papers dealing with circuit proper- 
ies, circuit philosophy, and techniques related to solid- 
tate devices in the following general areas: 


/ 


_ Novel semiconductor devices and conventional semi- 
conductor devices in novel modes of operation. 


PGEC News 


Minoru Nagata was born in Tokyo, 
Japan, in 1933. He graduated from the 
University of Tokyo with the degree of 
Bachelor of Engineering (Kogakushi) in 
1956. He then joined Hitachi, Ltd., where 
he has been engaged in the development of 
electronic analog computers in the Hitachi 
Central Research Laboratory. In particular, 
he has been working on the development of 
dc amplifiers. 

Mr. Nagata is a member of the Institute 
of Electrical Communication Engineers of 
Japan, and the Institute of Electrical Engi- 
neers of Japan. 


0 


Noah S. Prywes (M’55) was born on 
November 28, 1925. He received the B.Sc. 
degree in electrical engineering from The 
Technion, Israel Institute of Technology, 
Haifa, Israel, in 1949, From 1948 to 1950, 
he served in the Israeli Navy. In 1951, he 
received the M.Sc. degree in electrical en- 
gineering from Carnegie Institute of Tech- 
nology, Pittsburgh, Pa., and in 1954, re- 
ceived the Ph.D. degree in applied physics 
from Harvard University. 

Dr. Prywes worked from 1954 to 1956 at 
RCA, Camden, N. J. as a development 
engineer engaged in computing systems 
utilizing magnetics and transistors. From 
January 1956 to September 1958, he has 
been a department manager on the LARC 
project with Remington-Rand UNIVAC, 
Philadelphia, Pa. In this capacity he was 
responsible for all high-speed circuit design, 
interunit transmission, power supplies and 
control consoles for the LARC system. Later 
he was responsible for the design, develop- 
ment, construction and test of the comput- 
ing unit of the LARC system. Since last 
September, Dr. Prywes has been assistant 
professor of electrical engineering at the 
University of Pennsylvania, retaining a 
consulting position with Remington-Rand 
UNIVAC in Philadelphia. 

He is a member of AIEE and of Sigma 
Xi. 
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Frank L. Ragonese (A’46) was born in 
Bridgeport, Conn., on July 7, 1915. He 
graduated from the Bridgeport Engineering 
Institute in 1938 with a certificate in electri- 
cal engineering. 

After working as an electrical engineer 
at the Remington Rand Shaver Division in 
Bridgeport, Mr. Ragonese entered the U. S. 
Army and served from 1941 to 1946 as an 
Army Signal Corps Officer. While in the 
service, he attended radar and ultra-high- 
frequency courses at Harvard University and 
Massachusetts Institute of Technology. 

Since 1946, Mr. Ragonese has been em- 
ployed as an electrical engineer by the 
Remington Rand Laboratory of Advanced 
Research, South Norwalk, Conn., where he 
has engaged in the design, development, 
production and application of both black 
and white, and color industrial television 
systems, the design and development of a 
small magnetic drum unit, special permanent 
magnet type and electromagnetic type read 
and write heads for computer use, the de- 
velopment of a helicopter blade dynamic 
tracking system and, most recently, in the 
development of a magnetic character read- 
ing system. 
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N. Rouche (M’58) was born in Huy, Bel- 
gium, on June 2, 1925. He was graduated Ing. 
Elec. from the University of Liége, Belgium, 
in 1950. 

After a year as lecturer on the Faculty 
of Sciences, University of Liége, he did re- 
search on nonlinear oscillations at the Insti- 
tute for Mathematics and Mechanics, New 
York University, New York, N. Y. 

During the years 1953-1957, Mr. Rouche 
was an engineer in the Department of Auto- 
mation at the Bell Telephone Manufactur- 
ing Company, Antwerp, Belgium. In 1955, he 
was graduated Agrégé de 1’Enseignement 
Supérieur from the University of Liége. 

Mr. Rouche has been professor of ap- 
plied mathematics and electronics at 
Lovanium University, Leopoldville, Belgian 
Congo, Africa, since 1957. 


Significant contributions to the transistor circuit art. 


function. 
Solid-state 


microwave 

mechanisms. 
Solid-state memory, storage, and logic devices. 
Low-temperature digital and linear applications. 
Novel types of solid-state devices. 


Solid-state devices performing an integrated circuit 


devices and amplifying 


For further information, write to: Arthur P. Stern, 
Program Committee Chairman, 1959 Solid-State Cir- 
cuits Conference, General Electric Company, Elec- 
tronics Park, Building 3—Room 113, Syracuse, N. Y. 
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SCIENCE 


Vol. 1, No. 3 


PURPOSE 


It seems worthwhile to reiterate the purpose of 
SENEWS and to stress the need for your contributions 
if the program is to succeed. SENEWS is a newsletter 
addressed primarily to computer-oriented members of 
the IRE, ACM, and AIEE to aid them in promoting 
interest and knowledge in computing among _ high- 
school age students. It wants also to communicate di- 
rectly with secondary school teachers and other science- 
education groups, to report on their projects and prob- 
lems and offer help in bringing together all those inter- 
ested in such projects. 

The response to date has been quite gratifying, with 
the majority of items coming directly from teachers or 
leaders of other groups. We do, however, need and con- 
tinue to solicit your contributions and comments. 

Write to: 

Michael Warshaw, Chairman 

JCC Science Education Subcommittee 
The RAND Corporation 

1700 Main Street 

Santa Monica, Calif. 


CHANGING OF THE GUARD 


With the publication of this issue, we announce a new 
JCC Science Education Subcommittee Chairman and a 
new AIEE representative to the Editorial Board. 

Replacing Clare Farr as Chairman and Editor-in- 
Chief of SENEWS will be Michael Warshaw of The 
RAND Corporation. Robert Goldman of Bendix Com- 
puter Division is taking over from G. N. Hollander as 
the AIEE Editorial Board Member. 

As before, R. W. Melville of Stanford Research In- 
stitute and G. E. Forsythe of Stanford University will 
remain as representatives of the IRE and ACM, re- 
spectively. 

A vote of thanks to Messrs. Farr and Hollander for 
their efforts in getting SENEWS through its nascent 
phase. 


IRE TRANSACTIONS ON ELECTRONIC COMPUTERS 


December 


JOINT COMPUTER COMMITTEE 


SENEWS 


EDUCATION SUBCOMMITTEE NEWSLETTER 


December, 1958 


DISTRIBUTION AND PREPRINTS 


Volume 1, Nos. 1 and 2 of SENEWS were carried in 
the Communications of the ACM and IRE Transac- 
TIONS ON ELECTRONIC COMPUTERS. We are pleased to 
announce that the AIEE is going to distribute SENEWS 
separately to interested members and will make pre- 
prints available of No. 2 and all succeeding issues at 
their standard charge. Interested parties may obtain 
these by writing: 

R. S. Gardener 

American Institute of Electrical Engineers 
33 West Thirty-Ninth Street 

New York 18, N. Y. 


FUTURE ENGINEERS OF AMERICA—PHILADELPHIA 
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Daniel Ashler, from the Math Department of Abra- — 


ham Lincoln High School in Philadelphia, has reported 


on the activities of the group he sponsors: two decimal- — 


to-binary converters have been built, and an adder and 
multiplier are nearly finished. Ultimately, some sort of 
computer may be assembled. 

The designs for the decimal-to-binary converters are 
original: the two solutions to the problem are quite dif- 


ferent, and one of the boys has devised a circuit nota- _ 


tion for which he claims superiority over the conven- 
tional schematic diagrams. One of the converters won 
an outstanding award from the Ford Motor Company’s 
Industrial Arts Award competition. Components were 
gleaned from worn out pinball machines contributed by 


a local dealer and from scrap piles in junk yards except | 


for a few switches that had to be bought at a “war sur- 
plus” store. 

The boys are studying Boolean algebra and symbolic 
logic, and they have been taught to code for the IBM 
650. Last spring they worked out programs involving 
three-stage looping and ran them on the 650 at the local 
IBM office. 


Time has not permitted us to get all the details on 


the conversion schemes or on the new notation. They 


958 


hould both prove to be interesting, and we hope to 
ave more information for the next issue. 

| On the editorial side, we would like to commend Mr. 
ishler and his boys for showing what we feel is remark- 
ble restraint in not plunging right in with a full- 
lown computer project, but rather taking it piece by 
iece until a sufficient backlog of experience has been 
ained. Projects such as the design and construction of 
converter or an adder are fascinating tasks in them- 
elves and have the added feature of reaching comple- 
ion in a finite length of time. 


FUTURE ENGINEERS OF AMERICA—VAN Nuys 


|The FEA in Van Nuys, Calif., has recently started a 
omputer laboratory where each year ten students will 
ceive theoretical and practical training in digital com- 
juting. The technical training will be supplied by Don 
inderson and Jim Russell from the local IBM office. 
“wo evenings a week they will give instruction in com- 
juter logic, circuit theory, and programming. Sufficient 
fomponents will also be made available to give the 
itudents realistic laboratory experience. It is planned 
i) make this a permanent program, accepting ten new 
tudents each year, and Gerald Speen of the Valley 
council FEA estimates that a fully operating computer 
nay be realized in five years. 

There is, however, still a great need for professional 
celp at this section of the FEA. According to Mr. Speen, 
hirty-five students want to work in some phase of 
jomputing, but only ten could be selected. So any com- 
juter engineers who live in the valley and would like 
pb spend one night a week helping out a few interested 
iigh-school students should contact: 

Gerald Speen 

Future Engineers of America 

Box 2664 

Sepulveda, Calif. 

Phone—EM 72211 


Hi1GH-ScHOOL COMPUTER 


Aaron Buchman and his mathematics students at 
dutchinson Central Technical High School, Buffalo, 
N. Y., are shown in the accompanying photograph with 
he relay computer they designed and built (see SEN- 
£WS Vol. 1, No. 1). The relays are arranged on hinged 
»anels so that the student can actually see each step 
the computation as it is performed by the machine. 
Note the punched “paper” tape (9-inch wide press- 
soard) which the student feeds into the tape reader by 
and. In this manner, he controls the speed of compu- 
ation or stops it at any desired point in the program. 
“he conceptual advantages of this computer for ele- 
nentary instruction can hardly be overstated. Mr. 
3uchman staged a most convincing demonstration for 
four SENEWS editor last August, and is extremely 
nergetic and cooperative in personally answering in- 
juiries from others who wish to build a school computer. 
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F. W. Goro, © 1957 Time Inc. 


Pato ALTO COMMUNITY SCIENCE SEMINAR 


Although the following report from George Forsythe 
is not devoted primarily to the computing art, we feel 
it is of sufficient general interest and importance to 
print without deletion. 

“The Joe Berg Foundation (1712 South Michigan 
Ave., Chicago 16, Illinois) is described in the Reader’s 
Digest of August, 1958. It is subsidized to help communi- 
ties set up and run their own science seminars for gifted 
high school students with a real interest in and talent 
for scientific study. Things worked in Palo Alto some- 
thing like this: at the initiative of someone, probably 
Henry Martin, science teacher at Palo Alto High School, 
the school district contacted the Joe Berg Foundation. 
The foundation sent a number of standard tests and 
advice on how to use them to select appropriate stu- 
dents. In the spring of 1958, certain students of the 9th, 
10th, and 11th grades, recommended by their science 
or mathematics teachers, took the test. Of these stu- 
dents, 35 were selected. Meanwhile, Martin rounded up 
a comparable number of scientists in the area from in- 
dustry and university circles. 

“One night Jacob Shapiro of the Berg Foundation 
met with the scientists, the students, and their parents. 
He explained the general nature of the program, the 
need for regular meetings, and so on. He envisioned 
individual project work between youth and scientist, 
plus general lectures and outside reading. It was made 
clear that the Berg Foundation has no control over the 
program, which is run by the group itself. 

“Since last April the group has met from 7:30 to 10 
P.M. each Wednesday evening at Palo Alto High School, 
with no break for vacation. From 7:30 to about 8:30 
there is an expository lecture, usually by an adult scien- 
tist or engineer. After questions and discussion, the 
group breaks up into sections from 9:00 to 10:00. Sec- 
tions have been formed for biology, physics, chemistry, 
mathematics, electronics, geology, astronomy and pos- 
sibly others. 

“Here are some of the titles of the general talks: 
‘Radiological Safety Procedure,’ ‘The Electromagnetic 
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Structure of the Proton and the Neutron,’ ‘Deforma- 
tion of the Earth’s Crust,’ ‘Immunity,’ ‘Genetics,’ 
‘Instrumentation Chemical Analysis,’ ‘Nuclear React- 
ors,’ ‘Microwave Tubes,’ ‘Magnetic Resonance,’ ‘Prob- 
lems in Mathematics,’ etc. The speakers have included 
at least one Nobel prize winner, and numerous persons 
with an international reputation. 

“For the most part, the activities of the groups con- 
sist of individual projects worked on by one student and 
one adult—often in the adult’s laboratory. The follow- 
ing projects are typical of those proposed: study of 
variable stars, study of the effects of a nonlinear circuit 
element, study of an optimal automobile turn indicator, 
fluid flow visualization, a magnetic survey, and seasonal 
variation in quality of the Palo Alto water supply. 

“In the mathematics section we have so far empha- 
sized group study more than individual research, and 
are alternating two activities. Benjamin Epstein is 
leading the group in statistics using the excellent text by 
the College Entrance Examination Board, ‘Introduc- 
tory Probability and Statistical Inference for Secondary 
Schools.’ George Forsythe is leading the group in learn- 
ing to code and operate the IBM computer type 650. 
The IBM manual is the text, to be supplemented, when 
it arrives, by R. V. Andree’s ‘Programming the IBM 
650 Magnetic Drum Computer and Data Processing 
Machine’ (Holt, 1958). At the moment, the group is 
coding the Monte Carlo evaluation of the area of a 
circle, combining coding and experimental statistics. 
Computer engineers will be interested to know that 
Bill Kautz of the Stanford Research Institute discussed 
logical circuits and demonstrated his electronic com- 
puter components. 

“We are all learning how much bright youngsters can 
do, given some intellectual leadership by mature scien- 
tists. And the students are learning that membership 
in the group gives them real prestige in their schools.’’ 


UNIVERSITY OF OKLAHOMA TEACHER- 
STUDENT COMPUTER COURSE 


Because of the enthusiastic response given the 1957 


IRE TRANSACTIONS ON ELECTRONIC COMPUTERS 


December 


University of Oklahoma’s lecture series in modern 
mathematics for high school students and teachers, a 
full semester, noncredit course in computer program- 
ming and an introduction to related mathematics has 


been started for the coming school year. The history — 


behind this program has been communicated to us by 
Dr. Richard Andree from the University of Oklahoma. 

In September, a one-page announcement of the class 
was sent to high schools in the area, stating the goals 
and establishing October 4th as the starting date. The 
requirements were that no teacher would be admitted 
without an accompanying student and conversely. 
The announcement further made it clear that a great 
deal of homework would be required, the course would 
be free of charge, but the student would have to furnish 
his own transportation and text. Not everybody could 
be accepted. Classes would meet for a full day on alter- 
nate Saturdays, with the mornings used for laboratory 
work and the afternoons for lectures. 

Because of the competing activities at school, such as 
service teacher training, football, scheduled trips, etc., 
plans were made to accept no more than 30 students. 
Almost 300 applications were received! Hurried con- 
sultations and changes in plans allowed the program to 
be enlarged first to 90 and finally to 106 persons, as it 
was difficult to turn down students who were willing to 
commute 300 miles one way to take the course! The 
teachers were chosen first, then each teacher was per- 
mitted to select one of his students. 

A grant from the Oklahoma Frontiers of Science 
Foundation is providing the reimbursement for the 
instructor and two laboratory assistants, and the Uni- 
versity of Oklahoma furnishes classroom space and 
machine time. The text is R. Andree’s “Programming 
the IBM 650” (Henry Holt Co.). 

Dr. Andree reports that enthusiasm is running high, 
and the resulting programs are well written and often 
ingenious. He said that it is proving to be a rich, heart- 
warming experience to teach a group like this, and he 
sincerely recommends it to other computer installa- 
tions, 
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